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We review the available atomic data used for interpreting and modeling X-ray observations. The 
applications for these data can be divided into several levels of detail, ranging from compilations 
which can be used with direct inspection of raw data, such as line finding lists, to synthetic spectra 
which attempt to fit to an entire observed dataset simultaneously. This review covers cosmic 
sources driven by both electron ionization and photoionization and touches briefly on planetary 
surfaces and atmospheres. We review all of this, the applications to X-ray astronomy, the available 
data, recommendations for astronomical users, and attempt to point out the applications where 
the shortcomings are greatest. 
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Essentially every observation in astronomy which in- 
volves spectra requires some aspect of atomic data for 
its interpretation. This can range from stellar colors as 
probes of age or abundance to high resolution spectra 
of the Sun and the synthetic models which are used in 
their interpretation. The importance of atomic data to 
the field of X-ray astronomy comes from the fact that 
there are many atomic features in the X-ray band and 
the physical conditions in the sources are usually far 
from thermodynamic equilibrium, so that interpretation 
of line fluxes or ratios requires knowledge of collisional 
rate coefficients and radiative rates. Even the instru- 
ments with the greatest spectral resolving power are not 
capable of fully resolving many line blends, so that mod- 
eling is of particular importance. Furthermore, the field 
has evolved rapidly, and many of the needs for new or ac- 
curate atomic data have only recently become apparent. 
In this paper we point out the motivation for accumu- 
lation of atomic data, and review available atomic data 
used for interpreting and modeling X-ray observations. 

A review of atomic data for an application such as 
astronomy spans an audience including specialists with 
tools capable of providing new data, astronomers seeking 
to understand the accuracy and range of available data, 
and those who can use the data to synthesize spectra or 
other observables for comparison with observations. The 
subject differs from some other review topics since it is 
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challenge than by a single physical process. In this paper 
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we can. 

For the physicist reader, we attempt to describe the ca- 
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pabilities of current X-ray astronomy instruments, some 
of the astrophysical science challenges, and the scope of 
existing work on the measurement and calculation of rel- 
evant atomic constants. It appears that much of the re- 
cent calculation and experimental work which is most 
useful to X-ray astrophysics has been done when there 
is a clear understanding of the astrophysical application. 
We also attempt to draw attention to those who have 
made important or extensive contributions in the past, 
even in cases where their work has been superseded. In 
doing this, it is not feasible to be comprehensive. A great 
deal of work has been done in this area, and we refer the 
reader to the many extensive bibliographies which pro- 
vide a more complete set of references. For some physical 
processes we do provide relatively long lists of references 
to work which spans the many ions or atomic species of 
interest. These are cases where we think such lists are 
both relevant and useful, since they convey the volume or 
comprehensiveness of past work and also are of current 
interest. 

For the astronomer reader, we attempt to describe the 
technical challenges of calculating or measuring relevant 
quantities. The goal is to provide an appreciation of the 
uncertainties imposed on astrophysics by an incomplete 
knowledge of atomic quantities. In addition, in many 
cases, the data cited or described here can be directly 
applied to observations. 

For the reader interested in the synthesis of atomic 
data into model spectra or diagnostics which can be di- 
rectly fit to astronomical observations, we attempt to 
provide pointers to data which is appropriate for a given 
application, the limitations in accuracy and comprehen- 
siveness of existing data, and the prospects for future 
measurements or calculations. We also hope to foster 
cooperation between modelers or astronomers and the 
atomic physics community. This has proven to be an ef- 
fective way to stimulate work on problems which are rel- 
evant to astrophysics, while also educating astronomers 
and modelers about the limitations and applicability of 
atomic structure and cross section data. 

We feel that the time is right for a review such as this, 
because the accuracy and quantity of atomic constants 
describing many processes has progressed significantly in 
the last 10-15 years. This is driven both by advances in 
observational data with the launch of the Chandra and 
XMM — Newton X-ray astronomy satellites, both sta- 
tistical quality and spectral resolution, and by advances 
in the tools available to the physics community. There 
now exist values for rate coefficients and cross sections 
which will be accurate enough to allow interpretation of 
much of the available observational data for the near fu- 
ture, and current techniques will continue to improve the 
situation. In principle, the project of calculating or mea- 
suring atomic cross sections for any specific application 
has finite scope, and can at some point be considered 
complete. Although this is not imminent, we can esti- 
mate the quantity and accuracy of atomic data that will 
be needed for most applications to X-ray astronomy in 



the near future. 

In addition, the advent of on-line databases has rev- 
olutionized the dissemination of atomic data. There is 
no longer a need for publication of rate coefficients or 
cross sections in obscure laboratory internal publications 
or conference proceedings, and results based on propri- 
etary or fee-based databases or computational tools are 
receiving limited attention. The use of the web, faster 
computers and essentially unlimited storage have made 
it practical to save and disseminate great quantities of 
data. So there is no technical limitation to the shar- 
ing of data. In this review we therefore intentionally 
omit reference to unpublished data, data which is part of 
a fee-based or proprietary database or code package, or 
which is published in obscure laboratory internal reports 
or conference proceedings which are not widely available. 
We feel strongly that scientific information of all types, 
to be truly useful, must be freely and easily available, 
and its provenance must be transparent to all. 

Atomic data can be used in different ways according to 
the level of detail needed for a specific application. The 
levels of detail range from compilations used in direct 
inspection of raw data, such as line finding lists, to syn- 
thetic spectra which attempt to fit to an entire observed 
dataset simultaneously. The cosmic sources also differ, 
ranging from electron-ionized and photoionized plasmas 
to planetary surfaces and atmospheres. Although tradi- 
tional X-ray astronomy is generally agreed to span the 
energy range 0.1 - 10 keV, modeling of X-ray data in- 
evitably requires knowledge of cross sections which ex- 
tend to the ionization potential of hydrogen. We adopt 
this broader definition, and include some discussion of 
cross sections and rate coefficients associated with pho- 
tons emitted in the EUV and UV bands. X-ray spectra 
thus far are only capable of detecting the 15 or so most 
cosmically abundant elements, including C, N, O, Ne, Na, 
Mg, Si, Al, S, Ar, Ca, Fe, and Ni, and we generally re- 
strict our consideration to calculations and experiments 
relevant to these elements. 

Much of the information relevant to this topic has 
already been reviewed by others. Recent advances 
in astrop hysical X-ray spe c troscopy have be e n re- 
viewed bv iKahn et aD l)2004 . iPaerels and KahnI l)2003|) 
and some experimental results in atomic data by 
iBeiersdorfed |^03). The subject of atomic data for 
X-ray astronomy has been discusse d in specialized 
confer ence proceedings, including ISilver and KahnI 
lll993l). iBautista. Kallman and PradhanI 1)20001) . 
iFerland and SavinI l|200l|) and ISmithI l)2005l) . More 
specialized reviews pertaining to many of the subtopics 
discussed here will be mentioned in what follows. 

This review is organized as follows: We begin with 
a discussion of some examples of the astronomical X- 
ray spectra, choosing from recent observations using 
Chandra and XMM — Newton, followed by a review 
of the techniques used to derive atomic data, both the- 
oretical and experimental. We then discuss the atomic 
data currently available, classified according to the level 
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of modeling needed in order to apply the data. These 
range from line finding lists which can be applied by in- 
spection, to diagnostics of physical quantities involving 
a small number of measurable parameters such as line 
ratios or equivalent widths, to models which attempt to 
synthesize a large portion of the spectrum. These we di- 
vide into the likely excitation mechanism: coronal, for 
sources where the primary energy input is mechanical; 
photoionized, where photon excitation, heating and ion- 
ization plays a dominant role; and planetary, associated 
with neutral, molecular or solid material. Then we fur- 
ther subdivide according to physical process. 

We divide the discussion of each physical process or 
subtopic into sections: In the first, denoted 'background', 
we provide an overview of the work on the topic, dat- 
ing from the early suggestions of the importance of the 
process and simple approaches to estimating rate coeffi- 
cients and cross sections. In the second section, denoted 
'Recent Developments', we change tactics slightly, and 
provide an overview of the work which we consider to be 
the current state of the art, or which point the way to 
future developments or improvements in the state of the 
art. Inherent in this is the omission of some work in favor 
of that which we consider to be either of sufficient accu- 
racy and comprehensiveness to be adequate for the needs 
of X-ray astronomy today, or which we think points the 
way toward such data. 

The study of X-ray spectra associated with planets, 
comets, and other primarily neutral or solid objects is 
newer than the study of other emission mechanisms. As 
a result, there is less work so far on interpretive tools 
for such spectra, and we do not give a separate discus- 
sion of line finding lists or discrete diagnostics for them. 
We also include a glossary in order to help keep track 
of the many acronyms and terms describing computa- 
tional techniques, physical effects, and approximations. 
Throughout this paper we adopt chemical notation for 
ionic charge states, i.e. with a superscript denoting the 
ionic charge. This is in contrast to the spectroscopic 
notation used in much of the astrophysics literature, in 
which the charge state is denoted with a roman numeral: 
I=ncutral, Uncharge -1-1, etc. In the final section of the 
paper we discuss remaining needs for atomic data and 
prospects for future work. 



II. ON THE NEED FOR ATOMIC DATA 

The recent impetus for atomic data applications to as- 
tronomy comes from developments in instrumentation. 
Evolution of instruments has been rapid in the field of 
space astronomy, which has arisen and grown to a level 
of sensitivity and spatial and spectral resolution which is 
comparable to that of ground-based optical instrumen- 
tation in approximately the last 40 years. The most 
recent developments are summarized in the review by 
iPaerels and Kahnl l)2003 ^ for the X-ray band. Much of 
the current work is motivated by the launch of the grat- 



ing spectrographs on the Chandra and X MM — Newton 
satelli tes in the year 1999. These are described in de- 
tail in lBrinkman et aTl l|200Cll ) for the Low Energy Tr ans- 
mission Grating (LETG), in 'Caniza res et all pOOSj) for 
the High Energy Transmiss ion Grating (HETG) on the 
Chandra satellite, and by iRasmussen et all l)200lj) for 
the reflection grating spectrograph (RGS) on XMM — 
Newton. The gratings provide spectra with resolving 
power e/Ae ^300 - 1000 in the energy range between 
0.2 keV and 8 keV, with effective areas ~ 10 - 100 cm^. In 
terms of wavelength, the resolution is AA = 0.01 — G.OSA 
depending on the instrument, and the wavelength range 
extends from ~1.8 A for the Chandra HETG to 170 
A for the Chandra LETGS. This is in contrast with 
the previous genera tion of instrum ents, exemplified by 
the ASCA satellite l|TanakaL ITgsll . which carried mod- 
erate resolution CCD (Charge Coupled Device) instru- 
ments and had a spectral resolving power ~20. The grat- 
ing instruments have discovered line features in spectra 
previously thought to be featureless, and require a sig- 
nificantly more comprehensive and accurate database of 
atomic constants than had previously been in use. How- 
ever, both Chandra and XMM — Newton also carry 
CCDs whose use is necessitated by cosmic sources which 
are too faint for adequate signal in the gratings, or which 
arc spatially extended and hence cannot make use of the 
full grating spectral resolution. 

The Chandra and XMM — Newton gratings repre- 
sent a tremendous advance in resolution, but are not ca- 
pable of unambiguously resolving thermal or natural line 
shapes from many sources. The nominal resolution (a) of 
the HETG at the energy of the 0^+ Lyg line corresponds 
to a Doppler velocity of 170 km s~^ ijCanizares et al.l 
120051) . In a thermal plasma this would be obtained at 
a temperature of 10^'^ K, which is significantly greater 
than the temperature where this ion is predicted to be 
abundant under most plausible conditions. The natural 
width of the line corresponds to 0.008 eV or an equivalent 
Doppler velocity of 3.7 km s^^. Natural widths increase 
rapidly with the nuclear charge, and so for iron the cor- 
responding velocity is 550 km s~^. A typical dielectronic 
recombination (DR) satellite line to a hydrogenic reso- 
nance line (these terms are described in detail later in thi s 
review) is shifted by Ae/e ~ 0.002 l|Boiko et alill977|) . 
This is at the nominal resolution of the Chandra HETG. 
With good statistics, parameters describing line ccntroids 
or widths can be determined to precision better than the 
nominal resolution by factors of several. Nonetheless, 
it is clear that the Chandra spectra are not capable of 
measuring thermal line shapes unambiguously, and the 
wavelength precision is not comparable to that obtain- 
able in the laboratory or in solar observations with Bragg 
crystal spectrometers, for example. On the other hand, 
the HETG can resolve the important diagnostic lines dis- 
cussed in this review, and it can detect the effects of su- 
personic turbulence or bulk motions at Mach numbers of 
'--^ 2 - 3 or more. This is of great astrophysical interest, 
since many X-ray plasmas are located in regions of strong 



4 



gravity, or near shocks or other energetic sources. 

As an illustration of the capabilities and limitations of 
the data obtained with Chandra and XMM — Newton 
we present spectra obtained using these instruments 
which test the available atomic data and synthetic spec- 
tra. These are not chosen because they are typical, but 
rather because they best show the capabilities of the in- 
struments and the available atomic data. Figure ^ shows 
the spectrum of the Seyfert galaxy NGC 3783 obtained 
with the Chandra HETG. This object is an active galaxy, 
and the spectrum shows absorption from a large num- 
ber of resonance lines from the K shells of highly ionized 
ions of medium-Z elements, O - S, and also from the L 
shell of partially ionized iron. The solid curve is a model 
fit, described bv lKrongold et alJ l|2003j) . showing general 
consistency but not perfect agreement for the majority 
of the strong lines. In this object the gas is likely to 
be photoionized, and all the lines are shifted to the blue 
from their rest wavelength by an amount corresponding 
to a velocity of ^^500 - 1000 km s^^. This is interpreted 
as being a Doppler shift in an outflowing wind, and sim- 
ilar spectra have been observed from most of the other 
objects of this class. Discovery of such outflowing gas is 
surprising, since active galaxies are strong sources of con- 
tinuum X-rays, and are thought to contain black holes 
which radiate as gas is pulled inward. This spectrum 
also illustrates the fact that the data from the gratings 
on Chandra and XMM — Newton can be used to con- 
strain observed line wavelengths to within ^ 0.25%. This 
exceeds the precision of most ab initio calculations and 
therefore requires laboratory measurements to provide 
wavelengths precise enough to allow reliable line identi- 
fications and to infer line shifts. 

Figure |21 shows the spectrum of the nearby ac- 
tive star Cap cUa obtained with the Chandra HETG 
ijCanizares et al .. 2000). This shows rich line emission, 
similar to the Sun, and which is typical for a coronal 
source. For the most part the lines in the 10 - 18 A 
range are due to emission from L-shell iron ions. The 
observed iron L spectrum can be reproduced almost en- 
tirely b y assuming a single e lectron temperature of kT^ = 
600 eV ijBehar et al.U2001b|) . This temperature is consis- 
tent with both the measured fractional ion abundances 
of iron and with the temperature derived from ratios of 
Fe^^+ lines. However, there are some remaining discrep- 
ancies between single temperature models and the data 
for the lines of Fe^^+ near 16 A and also Fe^^"*" and Fe^^+ 
at 15.015 and 14.206 A which are all overestimated by 
current models. The statistical quality of this spectrum 
is comparable to that obtained from the Sun, illustrating 
the power of X-ray spectroscopy to study nearby objects 
outside the solar system. 

An illustration of the power of high resolution X- 
ray spectroscopy to study distant objects is shown in 
Figure |21 which shows the spectrum of a composite of 

J :alaxy clusters taken w ith the RGS on XMM — Newton 
Peterson et all l2003^ . The blue curve shows the data 
coadded for 13 clusters, which clearly exhibit coronal line 



emission from Fe^^"*", Fe^'^"'", and 0''+ Lya. This is consis- 
tent with emission from gas at an electron kinetic tem- 
perature of kTe ~ 1.5 keV. This temperature has also 
been derived from lower resolution observations of these 
and similar clusters. The gas densities can also be de- 
rived from the X-ray luminosities and sizes. Based on 
these quantities, the time for the gas to cool radiatively 
is predicted to be short compared with the probable clus- 
ter age, hence the lines from lower ionization species are 
expected as the gas cools to ~ 1 keV and below. The 
vertical dashed curves in the Figure denote the positions 
of the lines which are predicted from such cooling flow 
models, such as 0^+ and Fe^^+. Their absence requires 
a major reexamination of the prevalent model for cluster 
X-ray gas. One possibility is that the gas, rather than 
cooling steadily over cosmological time, is reheated by 
active galaxies in the cluster. If so, it would provide an 
additional example of the influence of black holes on their 
environments and on their likely influence on the growth 
of cosmic structures l)Ruszkowski et al.l[200i|) . 

Figure 0] shows the spectrum of a Seyfert galaxy 
NGC 1068 taken w ith the LETG on Chandra 
l|Brinkman et all l2002() . The extended long wavelength 
coverage is a unique feature of this instrument. The black 
curve shows the data taken of the nuclear region of the 
galaxy, and the red curve shows a model. The model 
consists of emission due to gas which is emitting primar- 
ily as a result of photoionization, photoexcitation and 
scattering, as evidenced by the radiative recombination 
continuum (RRC) features due to O^"'" near 16.5 A and 
C'^^ near 31.5 A and their strengths relative to the reso- 
nance lines from these ions. This is interpreted as being 
due to the fact that the galaxy contains a strong source 
of EUV and X-ray continuum radiation, likely a black 
hole, but that this source is hidden from our direct view 
by cold obscuring gas. This appears to be typical of 
type 2 Seyfert galaxies, which are likely intrinsically sim- 
ilar to the type 1 Seyfert galaxy NGC 3783 (Figure [ij. 
Both types contain a black hole partially surrounded by 
torus of nearly opaque cold gas, but their spectra reveal 
the fact that they are view ed at different orientations 
l|Antonucci and Milleil IT985^ . Thus X-ray spectra such 
as these provide evidence for the presence of hidden black 
holes via their scattered radiation. 

These examples illustrate the statistical quality of re- 
cent high resolution X-ray spectra, and, crudely, the ex- 
tent to which current synthetic spectra are able to match 
the observations. They also provide a glimpse of the as- 
trophysical issues which can be addressed using X-ray 
spectra and the potential importance of atomic data to 
address these issues. They also illustrate the unique na- 
ture of the information provided by X-ray spectra; the 
scientific insights provided by these observations are not 
accessible to any other waveband, or to observations us- 
ing lower resolution X-ray instruments. 
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FIG. 1 Chandra HETG spectrum of NGC 3783, from iKroneold et alJ ^ 
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FIG. 2 HETG spectrum of Capella, from ICanizares et al.l ll200Cf) . Insets show the regions of the spectrum in the vicinity of 
the lines of He-like Ne*^ and O^^. In the left inset panel the two curves are the data from the two arms of the dispersed 
image. The thin line is the data from the higher energy arm (HEG) and the thick line is the data from the medium energy arm 
(MEG). 
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FIG. 3 Composite spectrum of galaxy clusters where cooling gas is expected, from iPeterson et alJ i2003l) . The solid curve 
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FIG. 4 Spectrum of the Seyfert 2 galaxy NGC 1068, taken with the LETG on the Chandra satellite, f rom iBrinkman et alJ 
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A. Atomic Data Accuracy Requirements 

Spectra such as the ones shown so far in this section 
can be used as a practical guide for X-ray atomic data ac- 
curacy requirements. Simple estimates can be obtained 
by examination of the spectra. The strongest emission 
line (the Lya line of 0''+) in the spectrum of Capella, 
as shown in Figure [3 for example, has flux 0.02 photons 
cm^^ s^^, corresponding to a total of ~1300 counts in a 
3x 10^ s observation. The dominant error for such an ob- 
servation is counting statistics, and assuming Gaussian 
errors this corresponds to a statistical error on the line 
intensity of ~3%. This may be viewed as the implied 
accuracy requirement for a single atomic rate coefficient 
which affects the line flux linearly, if all other quantities 
which can affect the line flux are known or not of inter- 
est. Such rate coefficients, when they describe collisional 
processes, are derived by averaging the energy depen- 
dent cross section for the process over the Maxwellian 
distribution of electron velocities. So an accuracy re- 
quirement on a rate coefficient translates approximately 
linearly into the requirement on the energy dependent 
cross section if interpreted as applying uniformly over all 
energies. On the other hand, a greater uncertainty in the 
cross section can be tolerated if it applies over a narrow 
energy range. An observed line flux can be used to infer 
the total number of emitting ions if the excitation rate 
for the line is known. The accuracy of such a determina- 
tion will be limited by the atomic data if the uncertainty 
in the excitation rate exceeds the statistical error in the 
observed spectrum. 

On the other hand, inferring physical quantities such 
as gas temperature, density or elemental abundance from 
a line or a spectrum in many cases requires compar- 
ison with models. Model results depend on various 
atomic rate coefficients: those affecting the ionization 
balance, those affecting excitation and those affecting 
photon emission. The dependence of observables such 
as line strengths on rate coefficients affecting ionization 
balance or excitation can be non-linear, and are not well 
suited to analytic estimates for atomic data accuracy re- 
quirements. Numerical experiments can illustrate the 
sensitivity of line strengths to the rate coefficients af- 
fecting ionization and excitation. As an illustration. Fig- 
ure shows the ionization balance in coronal equilib- 
rium for iron. Ion fractions are plotted as a function of 
temperature T in the range 4<log(T)<8 and the dom- 
inant ionization stage ranges from Fe^ at low tempera- 
ture to Fe^^+ at the highest temperature. The left panel 
show s these quantities calcu l ated using rate coefficients 
from lArnaud and Ravmondl l)l992|) for iron. We denote 
these as our baseline ionization balance (but note that 
the work by,Brvans ct al. ( 2006) represents a significant 
update). The right panel has been calculated using the 
baseline rate coefficients, except that the rate coefficients 
for recombination for each ion have been perturbed, i.e. 
multiplied by a randomly chosen factor ranging between 
0.79 and 1.28 (i.e. A log[ratecoeflicient] = ±0.1). This 



range is only slightly greater than the magnitude of the 
effect which is likely to be associ ated with pl asma mi- 
crofields, cf.. Section IVlTX2l and iBadnell et a l. (2003), 
and so reflects a realistic limit to the level of accuracy 
needed from zero-external-field rate coefficients for this 
process. Within this range, the multiplicative factor af- 
fecting recombination is randomly chosen and evenly dis- 
tributed. The recombination rate coefficients are stored 
as analytic functions of temperature and, for each ion, the 
same multiplicative factor is applied at all temperatures. 
Comparison of the two panels of Figure [S] indicates the 
sensitivity of the ionization distribution to changes in the 
recombination rate of this magnitude: at a given temper- 
ature ion fractions can differ by factors ~2-4, although 
closed shell ions such as Ne-like and He-like ions are rela- 
tively unaffected. The temperature where the maximum 
fraction occurs can also change by ~0.1 - 0.2 in log(T). 

The effects of the ionization balance change on the 
synthetic spectrum, and on the inferred distribution of 
gas temperatures in a real astrophysical source, are il- 
lustrated in Figures and [7| Figure shows an ob- 
served spectrum from the active star Capella taken with 
the HETG grating on the Chandra satellite, the same 
source shown in Figure [3 The observed spectrum is 
shown as the data points with error bars. The solid 
(red) curve in Figure El is a synthetic spectrum calcu- 
lated using the baseline ionization balance shown in the 
left panel of Figure [S] This spectrum is calculated using 
the ionization balance at two temperatures, log(T)=6.9 
and log(T)=7.1. The element abundances and the rela- 
tive normalizations of the two components were allowed 
to vary in order to achieve a fit which has x^=3267 for the 
1602 independent energy channels shown in this figure. 
This is not in a range of which permits application 
of traditional statistical measures of probability of ran- 
dom occurrence (Bevington, 1969), and is primarily of 
illustrative value. Similar values of are found when 
these data are fit to two temperature component mod- 
els calculated with other modeli ng packages which a re in 
widespr ead use, such as AP EC (jSmith et all l200l|) and 
MEKAL (iMewe et al.Ul995b|) . 

In the region of the spectrum shown most of the fea- 
tures are due to iron, the exception being the Lyman 
lines of Ne^+. We find a best fi t value of [Ne/Fe], th e 
abundance ratio relative to solar (jGrevesse et allll996l) . 
of 0.66. The strongest lines in this spectrum are indicated 
in Figure Fe^^+ between 12 and 17 A, and the Lya 
line of Ne^+ at 12.13 A, which is blended with the corre- 
sponding Fe^^+ line. Higher ionization lines include lines 
from Fe^^+ at 10.5 - 10.8 A, Fe"+ between 11 and 16 
A, Fei8+ at 13.5 - 13.8 A and Fe^^+ at 12.85 A. Lower 
ionization lines include 0*"+ at 16.00 A, Fe^^+ blended 
with Ne9+ L/3 at 10.24 A and Ne*+ at 11.03 A . The 
baseline model accounts for the strengths of many of the 
strongest lines, but under predicts the Fe^^"*" 15 and 17 
A lines; this is a manifestation of problems with rate co- 
efficients affecting this ion which are in widespread use 
and which are discussed in section FVl. A. 41 



8 



The emission line spectrum radiated by the plasma 
is calculated using rate coefhcients for electron im- 
pact coUisional excitation culled from the references dis- 
cussed later in this paper. These are enumerated in 
iBautista and KallmanI l)200l|) . with the addition of rate 
coeffici ents from th e recen t CHIANTI version 5 compi- 
lation l)Landi et all |200^. It is well known that a 
source such as Capella is best fit using a broad distri- 
bution of temperature s, ranging from log(T)~6 - 7.4 
l|Canizares et all 1200(1) . but we adopt the two compo- 
nent model shown here for illustrative purposes. 

The effects of perturbing the ionization balance on the 
fit to the spectrum can be seen in Figure |3 in which 
the analogous procedure is carried out using the ioniza- 
tion balance shown in the right panel of Figure In 
making this fit we have not attempted to reoptimize any 
of the parameters describing the fit, i.e. the normaliza- 
tion, temperatures or elemental abundances in the two 
temperature components. Rather, we have left them the 
same as in Figure in order to illustrate the sensitivity 
to the ionization balance. Comparison with the fit shown 
in Figure El shows that the fits to the strongest lines, i.e. 
those of Fe^^"*" and Ne®+ , are essentially identical. This is 
a reflection of the fact that the ionization balance of the 
closed-shell Fc^^+ is little affected by the perturbation 
in recombination rate coefficients. However, significant 
differences are apparent in the fit to the ions in adjacent 
ionization stages, notably Fe^^"*" and Fe^^+. These show 
up in the lines near 14.2-14.5 A for Fe^''+ and 12.8 A for 
Fe^^"*" , which are under predicted by the model shown in 
Figure [7| The lines of Fe^^"*" are not strongly affected by 
the perturbation in the recombination rate coefficients. 
The fit between data and model in FigureElhas x^=3610 
for 1602 energy channels and is clearly worse than for 
the baseline model. However, when the perturbed model 
is iteratively fitted to the Chandra data a better fit is 
found when the high temperature component increases 
to log(T)=7.2. The resulting fit is improved, although 
still inferior to the baseline model: x^=3522 for 1602 en- 
ergy channels as compared with x^=3267 for 1602 energy 
channels for the baseline model. 

These experiments illustrate some of the effects of 
changes in the rate coefficients affecting ionization bal- 
ance on fitting to observed astrophysical X-ray spectra 
and on the inferred temperature distribution. They in- 
dicate that changes of ~25% in these rate coefficients 
can affect the model strengths of the lines of the more 
delicate ions by factors ~3 - 10, while the strong lines 
from stable ions are more robust. Such changes in model 
line strengths lead to changes in the inferred physical 
conditions derived by iterative model fitting which are 
significant. In this case we find changes of 0.1 in log(T) 
from such a procedure, which is at the resolution of the 
grid of models used to calculate the spectra. 

In a coronal plasma there is an exact symmetry be- 
tween recombination and electron impact coUisional ion- 
ization, so that perturbations to recombination rate coef- 
ficients can be interpreted as corresponding inverse per- 



turbations to the ionization rate coefficients. The in- 
ferences from these numerical experiments are recipro- 
cal between the two types of rates. At the same time, 
it is clear that the above model fits have been applied 
to only one spectrum, and only one numerical experi- 
ment at perturbing the rate coefficients has been per- 
formed. In order to adequately characterize the effects 
of perturbing various rate coefficients, this experiment 
should be repeated in order to truly randomize the per- 
turbations, and experiments should be carried out with 
different amplitudes and with fits to other astrophysical 
X-ray spectra. More extensive experimentation of this 
type has been carried out with application to solar EUV 
and X-ray line emission b jjSavin^jn^Lamir^^ and 
iGianetti. Landi and Landiniri|200o|riMasa11^ 
ined the effects of uncertainties in rates affecting ioniza- 
tion balance on the results of fitting to moderate resolu- 
tion X-ray spectra obtained using CCD instruments. 



III. THEORETICAL TECHNIQUES 

Theoretical calculation of atomic processes begins with 
the solution to the multi-electronic Schrodinger or Dirac 
equation, 

j2K+j2v{fA^ = E^. (1) 



where hi are the one-electron Hamiltonians and Vf, the 
electron-electron interaction potentials. In the fully rela- 
tivistic case, hi are Dirac Hamiltonians, Vfj~'^ include the 
direct and exchange part of the electrostatic interaction 
(the interaction energy due to the positional correlation 
of parallel-spin electrons according to the Pauli exclu- 
sion principle) and the Breit interaction (the two-body 
magnetic interaction), and 5" is a four-component wave 
function. In the non-relativistic case, hi are Schrodinger 
Hamiltonians and V^~'' include only the electrostatic in- 
teraction. In the Breit-Pauli (BP) relativistic approx- 
imation, the non-relativistic multi-electronic Hamilto- 
nian is corrected by adding terms resulting from the re- 
duction of the one-electron Dirac Hamiltonian and the 
Breit interaction to the Pauli form ijBethe and Saloeteil 
I1972D . The relativistic correction terms are then the 
mass-correction, Darwin, and spin-orbit terms that are 
added to hi, and the spin-other-orbit, spin-spin, orbit- 
orbit, spin-contact, and two-body Darwin terms that are 
added to Vlj~'^- Higher-order relativistic interactions like 
QED have also to be treated in the case of highly-ionized 
heavy atoms. The angular part of the solution of the 
Schrodinger equation is known exactly and is found us- 
ing the Racah algebra l)Edmond^ll96(]l:lFano and RacahL 
119591: iRacabL Il942l Il943|) . With regards to the radial 
part, the exact solution cannot be found due to the 
V^^!~'^ terms. One of the simplest approximation is the 
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FIG. 5 Left Panel: Ionization balance for iron in a coronal gas, calculated using rate coefficients described in the text. Right 
panel: Ionization balance calculated using perturbed rate coefficients as described in the text. Curves correspond to various 
ions of iron: in the left panel Fe+ dominates at log(T)=4, Fe"""^ dominates at log(T)=5.2, Fe+* dominates at log(T)=5.9, Fe"""^^ 
dominates at log(T)=6.3, Fe^^'' dominates at log(T)=6.6, Fe^'^* dominates at log(T)=7.5. In the right panel, the temperature 
of maximum abundance of Fe"""® - Fe^^'' is displaced to lower temperature by ~0.1 in log(T). 
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FIG. 6 Fit to a 3x10" s Chandra HETG observation of the 
active star Capella in the 10 - 18 A region using the baseline 
ionization balance shown in the left panel of figure |S] 
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FIG. 7 Same as figure |^ except using ionization balance 
shown in figure |S| derived from perturbed rate coefficients as 
described in the text and displayed in the right panel of figure 
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independent-particle or central-field model where 



e — e 



(2) 



where Vi{r) are local central potentials in which each 
electron moves independently. The atomic state function 
is then a product of the atomic orbitals which are 
eigenstates of the mono-electronic Hamiltonians, 



{h, + Vi{r)}(P,^ ei<j)i, i = l,...,N 



(3) 



where N is the number of electrons. Different forms 
of the functions Vi{r) can be considered l)Cowanl IlljSl'l. 
Forms that depend upon the orbitals themselves lead to 
the u se of the self-c onsistent- field (SCF) iterative proce- 
dure l|CowanL[T98l|) . 

A more elaborate solution of Eq. |^ comes from the 
application of the variational principle. Using the trial 
atomic state function defined by 



^ = {N\)-^''^det (f>i{xj) 



(4) 



where the electron space and spin coordinates and 

the determinant insures the Pauli exclusion principle, the 
are obtained by requiring that the expectation value 
of the multi-electronic Hamiltonian is minimum 



(5) 



where the variation is taken with respect to the radial 
parts of th e orbitals (t>i. This leads to the Hartree-Fock 
equations l)Hartred . ll957|) 



+ I dr'Mx'), 
j 



i^l,...,N. (6) 
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These non-homogeneous coupled integro-difFerential 
equations are solved iteratively in a SCF procedure. An 
important property of the Hartree-Fock approximation 
follows from the Brillouin's theorem which implies that 
the diagonal matrix elements of single-particle operators 
are given correctly to first order by Hartree-Fock atomic 
state functions. This approach essentially takes into ac- 
count the part of the electron correlations within an elec- 
tronic configuration. 

Further improvement can be made using trial atomic 
state functions 4" that are expansions of configuration 
state functions $, 

^-^^Cfc^fe (7) 

k 

where are themselves expansions of Slater determi- 
nants (Eq. 10}) and are the mixing coefhcients. The 
summation in Eq. Q can be systematically extended 
in principle to yield results of arbitrarily high accuracy 
or to achieve convergence. This approach is referred 
to as the configuration interaction (CI) approximation. 
When the variational principle (Eq. ||SJ)) is applied to 
Eq. ((7J varying the Cfe and the orbitals (pi , one talks about 
the multi-configuration approach su ch as in the m ulti- 
configuration Hartree-Fock (MCHF) llFischeilll996D and 
multi -configuration Dirac-Fock fMCDF'l ijParpia et all 
[1996') methods. Otherwise one talks rather about the su- 
perposition of configurations approach where Eq. (|5Jl or 
Eq. © is used to determine the orbitals cf>i and the multi- 
electronic Hamiltonian matrix is diagonalized afterward 
to obtain the Ck and the eigen value E as, for example, 
in Cowan's code ijCowanl Il98l|) . In practice, the CI ex- 
pansion in Eq. ITj) is limited and in order to include more 
correlation in the multi-electronic Hamiltonian matrix 
elements are corrected to reproduce the few available ex- 
perimental energy levels. The approximation becomes 
then a semi-empirical one. 

A different approximation which is also capable of 
systematic improvement of the atomic state function is 
based on the man y-bo dy perturbation t heory (MBPT) of 
iBruecknei ^ (1955) and' Goldsto^ l)l957|) . In this method, 
the multi-electronic Hamiltonian is written as a sum of a 
zero-order Hamiltonian and a perturbation term 

H = Hq + Hpert (8) 

where 

i 

and 

^pert=E^r'-E^'W- (10) 

i<j i 

A complete set of orbitals (pi are then obtained by solv- 
ing Eq. (O and used in the order-by-order Rayleigh- 
Schrodinger perturbation expansion of the atomic state 



energies E and other properties such as oscillator 
strengths. An example of the reduction of the per- 
turbation expansion to MBPT formulas using Feynman 
diagrams and second-qu antization method is given in 
lAvgoustoglou et alJ l)l99 2') for closed-shell atoms. 

For processes involving the continuum, the continuum 
atomic state function ^'(e) for the (A^-l-l)-electron system 
has to be evaluated, 

a 

where e is the energy of the free electron, (pa are free- 
electron orbitals and Xq sue atomic state function of 
the A^-electron target Hamiltonian. The free-electron or- 
bitals are solutions of the following Schrodinger equations 

{h^ + e}c^a^^Uc.c'(pc' (12) 

a' 

where are one-electron kinetic Hamiltonians and the 
matrix potential Uaa' is defined by 

Uaa' = J dxi...dxN Xc^i.Xi...XN) /^g^ 
X U{xi...xnxn+i) Xa'{xi.-.XN) 

where U is the sum of the nuclear and electron-electron 
interaction potentials acting on all the {N + 1) electrons 
of the target-plus-free-electron system. 

The widely used distorted wave (DW) approximation 
consists in neglecting the coupling between the different 
channels or the interaction between continuum states, i.e. 

Uaa'^0,a^a', (14) 

by using the central- field model defined in Eq. . It be- 
comes a good approximation for both the non-resonant 
background and resonance contribution. However, many 
DW excitation rate coefficients omit the resonance con- 
tribution entirely. 

At sufficiently large ionization or free-electron energies, 
DW orbitals approximate closely to Coulomb orbitals cal- 
culated using 

U^^-^-i^ (15) 

r 

where Z is the atomic number. It is then called the 
Coulomb-Born (CB) approximation. The exchange part 
of the electrostatic interaction is neglected in CB approxi- 
mation but included in the Coulomb-Born-Oppenheimer 
(CBO) and Coulomb-Born-Exchange (CBE) approxima- 
tions. 

The simplest approximation consists in considering 
plane waves to describe the continuum orbitals. It is 
equivalent to switching off the matrix potential, i.e. 
Uaa ~ 0. This technique called the Born approxima- 
tion is most likely to be accurate for energies far above 
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threshold, and for an electron colliding a neutral atom. It 
is a good approximation for the background cross section 
over a wide range of energies for positive ions, because 
the infinite energy limit for scattering of a positive ion is 
exactly the Born cross section and this limit controls the 
finite energy cross section over a wide range of energies 
due to the slowly varying behavior of the cross section 
with energy. 

The three above-mentioned approximations make as- 
sumptions about the wave function of the continuum elec- 
trons which are only appropriate when the energy of the 
projectile electron (or free electron) is much greater than 
the binding energy of the target electrons (or bound elec- 
trons) of interest. In contrast, approximations such as 
the more elaborate an d wide ly used i?-matrix method 
l)Burke and BerringtonL Il993|) can be invoked for free 
electrons at low energies and for low ionization charges 
up to negative ions where the coupling between channels 
is generally strong. 

Moreover, the importance of resonant enhancement 
of rate coefficients for collisional excitation (see e.g. 
Section IVII A.ll means that the R-matrix method is 
widely used for highly-charged ions as well, as it is a 
very efficient way of determining the resonance contri- 
bution. Distorted-wave cross sections are traditionally 
non-resonant. There are two ways in which resonances 
can be taken account of subsequently. Firstly, (as ex- 
emplified by HULLAC) autoionizing states can be in- 
cluded in the coUisional-radiative equations when calcu- 
lating level populations and effective recombination and 
ionization rate coefficients. This can become very de- 
manding computationally but it does allow for collisional 
redistribution of autoionizing states in dense plasmas. 
Secondly, the low-to-medium density approach, appro- 
priate to the solar atmosphere, magnetic fusion plasmas, 
is to include the resonances in the excitation calculation, 
either via R-matrix (always present), or perturbatively 
(as has been done for some ions using FAC) and so omit 
autoionizing states from the coUisional-radiative popu- 
lation rate equations (otherwise they would be double- 
counted). This hugely simplifies the coUisional-radiative 
problem but means there is a density range below that 
for collisional equilibrium, where partial collisional redis- 
tribution of autoionizing states takes place, and where 

non-autoionizing level populations are not well described 
1 

The R-matrix method consists of dividing the space 
into an internal and external regions. In the internal 
region where the interaction between the free electron 
and the target electrons is strong, the {N + l)-system 
atomic state function ^'(e) is expanded in terms of target 



^ The authors thank Dr. N. Badnell for pointing this out and 
providing the text for this paragraph. 



eigenfunctions X: 

V&(£)^^^Xa^^+E^^*' (16) 

a i 

where A is the anti symmetrization operator insuring the 
Pauli exclusion principle, the bar sign over x indicates 
the coupling with the angular and spin parts of the free- 
electron orbital 0(e), Ue{r)/r represents the radial part of 
the latter, and <I>j are bound states of the (iV-|- l)-system 
that are constructed with target orbitals to ensure com- 
pleteness of 4'(e) and improve short-range correlations. 
The functions Ue{r) and the coefficients Ci are obtained 
applying the Kohn variational principle to Eq. (|16|) giving 
rise to a set of coupled integro-differential equations. 

These are solved by determining the i?-matrix which 
is a real, symmetric matrix whose rows and columns are 
labeled according to the possible channels of the scatter- 
ing problem. At this point, one should notice by looking 
at Eq. (|16|l that resonances (quasi-bound states) are au- 
tomatically included either in terms of closed channels 
in the first expansion or as linear combinations of bound 
states in the second expansion. In the external region, 
^[e) is represented by Eq. l(TT|l and Eq. ((T^ is solved 
with the following long-range matrix potential, 

C/a.'(r)=^citlA'+^ (17) 

A 

The MBPT method described above is also a power- 
ful technique to compute continuum processes. Here, the 
perturbation expansion is reduced using Feynman dia- 
grams involving, for instance, the photoelectron of the 
photoionization process. 

We now describe some specific codes for computations 
of atomic structure and scattering. Codes employing 
similar assumptions and computational algorithms are 
grouped together. 

-HULLAC: 

The Hebrew University Lawrence Livermo re Atomic 
Code (Bar-Shalom. Klapisch and Ored. l200ll) is a full 
package for structure and scattering calculations. It 
solves the multi-electronic Dirac equations and includes 
QED corrections. An analytical parameterized central- 
field potential is used to obtain the analytical Slater-type 
spin-orbitals. The atomic states 'I' are constructed using 
CI expansions. The parameters appearing in both the 
central-field potential and the spin-orbitals are varied in 
a MC approach in order to minimize the energies of any 
set of levels or, in its semi-empirical mode, to minimize 
the rms deviation between a set of experimental energy 
levels and their corresponding theoretical values. The 
continuum is treated in a DW approximation. Using 
interpolation techniques for the radial part of collision 
strengths, it can calculate large quantities of data. 

-FAC: 

The Flexible Atomic Code (|g3, l2003a|) is also a fully 
relativistic program computing both structure and scat- 
tering data. It uses a modified Dirac-Fock-Slater central- 
field potential which includes an approximate treatment 
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of the exchange interaction. The orbitals are optimized 
in a SCF iterative procedure during which the average 
energy of a fictitious mean configuration with fractional 
orbital occupation numbers is minimized. This mean 
configuration represents the average electron cloud of the 
configurations retained in the CI expansion. This is an 
intermediate approach between MC and SC. The DW 
method is used to calculate continuum atomic state func- 
tions '^{e). 

-MCDF /grasp /GRASP92 : 

The Multi - Confi guration Dirac-Fock structure code 
( Grant et'all llOSCj) and its more re cent versions, GRAS P 
(D^iO^U Eilll and GRASP92 l)ParDia et all Il996|) . 

use the elaborated fully relativistic Hartree-Fock MC 
variational principle to obtain the Ci and four-component 
(pi in the CI expansion defined by Eq. One energy 
level or the average of a group or all energy levels is min- 
imized in the orbital optimization. The QED interaction 
is treated as corrections to the energy levels. 

-SUPERSTRUCTURE / AUTOSTRUCTURE / AU TOLS J: 

SUPERSTRUCTURE (jEissner et al.l Il974^ solves the 
Breit-Pauli Schrodinger equation. CI atomic state func- 
tions are built using orbitals generated from the scaled 
Thomas-Fermi-Dirac central-field potential that is based 
on the free-electron gas with exchange approximation. 
The scaling factors for each orbital are optimized in a 
MC variational procedure minimizing a LS term energy 
or weighted average of LS term energies. Semi-empirical 
corrections can be applied to the multi-electronic Hamil- 
tonian in which theoretical LS term energies are cor- 
rected in order to reproduce the centers of gravity of 
the available experimental multiplet s. Oth er versions are 
AUTOSTRUCTURE LBadnell, 1986a, Il997|) and autolsj 
ijPubau et all Il98l|) . These include the treatment of the 
continuum in a DW approach. In AUTOSTRUCTURE, an- 
alytical Slater-type orbitals can be also used, and non- 
orthogonal orbital basis sets can be considered for the 
calculations of inner-shell processes. 

-Civ3: 

The Config uration Interaction Version 3 code of 
iHibbertI l|1975|) uses analytical Slater-type orbitals in CI 
expansions to represent atomic state functions and solves 
either the non-relativistic or the Breit-Pauli Schrodinger 
equations. The orbital parameters are varied in a MC 
optimization procedure to minimize one or several en- 
ergy levels using an analytical central-field potential. 
The atomic state function can be improved in a semi- 
empirical procedure called fine tuning in which the multi- 
electronic Hamiltonian matrix elements are corrected in 
order to reproduce the available experimental energy lev- 
els. 

-R- MATRIX / BPRM / DARC : 

The BPRM scat te ring package described in 
iBurke and BerringtonI l)l993f) implements the elabo- 
rated i?-matrix method in the Breit-Pauli relativistic 
approximation. In the latter, only the one-body Darwin, 
mass-correction and spin-orbit terms are considered. 
The scattering package can use either superstructure. 



AUTOSTRUCTURE Or Civ3 target orbitals. The radiation 
and Auger dampings which are important in processes 
involving i nner-shells have been included using an optical 
potential feorczvc a and BadnellL Il996l l200Clf) . A fully 
relativistic implementation of the i?-matrix approach 
using GRASP target sp in-orbitals is also available (darc; 
lAit-Tahar eFall |l996j) ). 
-MCHF: 

In the Multi-Configura tion Hartree-Fock structure 
code (jFischer et al.l Il997j) . the non-relativistic multi- 
electronic Schrodinger equation is solved using the elab- 
orated Hartree-Fock MC variational principle to obtain 
the orbitals and the mixing coefficients utilized in the CI 
expansion of the atomic state function. The relativistic 
corrections are included by diagonalizing afterward the 
Breit-Pauli multi-electronic Hamiltonian. 

-HFR: 

The pseudo-relativistic Hartree-Fock code of ICowanI 
l)l98lj) solves the Hartree-Fock equations for the spher- 
ically averaged atom for each electronic configurations. 
These equations are the result of the application of the 
variational principle to the configuration average en- 
ergy. Relativistic corrections are also included in this set 
of equations, i.e., the Blume-Watson spin-orbit, mass- 
variation and one-body Darwin terms. The Blume- 
Watson spin-orbit term comprises the part of the Breit 
interaction that can be reduced to a one-body operator. 
CI is taken into account in an SC approach. The ra- 
dial parts of the multi-electronic Hamiltonian can be ad- 
justed to reproduce the available energy levels in a least- 
squares fit procedure. These semi-empirical corrections 
are used to allow the inclusion of higher-order correla- 
tions in the atomic state functions. Continuum orbitals 
are calculated in a DW approach using HXR central- 
field potential which uses an approximate exchange part 
of the Hartree-Fock non-local potential. 

-MZ: 

This code is base on the Z expansion technique 
l|Safronova and Urnovl [1980) in which a MBPT is ap- 
plied using screened hydrogenic orbitals. The perturba- 
tion expansion then contains powers of Z, the atomic 
number. The Breit interaction and QED are also treated 
in this expansion. Convergence is best with this method 
for highly ionized ions, and all members of an isoelec- 
tronic sequence are treated simultaneously. 



IV. EXPERIMENTAL TECHNIQUES 

Experiments have the potential to directly simulate in 
the laboratory the same processes observed remotely us- 
ing telescopes, and so can provide the most accurate and 
appropriate data for use in interpreting astronomical ob- 
servations. This cannot always be achieved, owing to 
practical challenges of achieving the densities or radia- 
tion environments occurring in some sources, and also 
because the scope of the data needs exceeds the avail- 
able experimental resources. So in many cases we rely 



13 



on experimental results as crucial checks on computa- 
tions, while for many other quantities of interest exper- 
iments represent the only means to determine atomic 
quantities with sufficient accuracy for comparison with 
observations. This accuracy is needed in the case of line 
wavelengths or energy levels, where inaccurate values can 
lead to line misidentification or misleading redshifts or ra- 
dial velocities. Experiments are also crucial for checks of 
the accuracy of DR rate coefficients at low temperature 
(~ lO''^ K), owing to the fact that this process requires 
energy levels which are accurate to ~1 cV. 

Experimental techniques have progressed at a rate 
comparable with the development of theoretical tools 
and computer platforms. Notable in their applica- 
tion to X-ray processes are the electron beam ion 
trap (EBIT), wh i ch ha s been described in detail by 
iBeiersdorfer et all lIlQQOl ) , and t he storage rings d e scribe d 
by lAbrahamsson et al.' fl993), ISchippers et all l)l998(l . 
and [G winner et a l. ( 2001). These dev e lopme nts have 
been reviewed recently by iBeiersdorfed l)2003|) and we 
will not repeat them here. In what follows we will de- 
scribe, for each process, which experiments have con- 
tributed to the available database and what future de- 
velopments might include. 



V. LINE FINDING LISTS AND RELATED TABULATIONS 

In order to organize the many different types and appli- 
cations for atomic data, we divide them according to the 
level of computation required for their use. Atomic data 
which is simplest to use, in principle, includes things like 
line finding lists, which can be utilized with a minimum 
of detailed fitting or model calculation. Such tabulations 
are crucial to interpretation of datasets taken at high res- 
olution. They allow a crude determination of the likely 
conditions in the gas under study, and may be able to ad- 
dress simple questions about element abundances and the 
likely range of plasma conditions. In addition, they are 
needed for successful application of the tools described in 
the later sections. These differ according to the physical 
excitation mechanism, and we discuss them separately. 



A. Coronal Plasmas 

Background 

Coronal plasmas are those in which the ionization 
balance is determined primarily by collisions between 
ions and thermal electrons. They are common to situ- 
ations involving mechanical heating, such as supernova 
remnants, stellar coronae and clusters of galaxies. For 
our purposes this includes objects such as supernova 
remnants which may not be in ionization equilibrium 
l|Shklovskvl Il968f) . A feature of coronal plasmas which 
are in or near equilibrium is that the ionization balance 
adjusts so that the ionization potentials of the most abun- 
dant ions are comparable to the gas temperature. As a 



result, the electrons have sufficient energy to excite many 
atomic bound states, and so the spectra are generally rich 
in line emission. For coronal plasmas, the most useful low 
level tools are likely to be finding lists for emission lines. 
This typically includes resonance lines, lines which have 
the ionic ground term as the lowest level, since these are 
often the strongest and therefore most easily measured. 
It also includes satellite lines, which resemble resonance 
lines but with the addition of a spectator electron, since 
these can be used to infer temperature. Finally, it in- 
cludes electric-dipole forbidden lines and lines connecting 
excited levels (subordinate lines) which are potential di- 
agnostics of various plasma parameters such as electron 
density and the effects of cascades. 

Prior to the launch of Chandra and XMM — Newton 
the need for accurate X-ray line wavelengths in astron- 
omy was driven by the study of solar X-rays. This 
led to the development of line lists in close cooperation 
with laboratory spectroscopists. Early work originated 
from data taken by rocket fiights and the OSO satel- 
lites. These include rocket observations of first solar de- 
te ction of 0^"'' lines and oth e rs in the 14 - 22 A range 
bv lFreeman and Jones! l|l970j) . IWalker and Ruggd l)l97(i|) 
presented observation of the forbidden lines in the He 
isosequence. iMcekins et al. ( 1968) presented fiare spec- 
tra from crystal spectrographs on OSO-4 in the wave- 
length ranges 0.5 - 3.9 A and 1.0 - 8.5 A. An analysis 
of the sol ar spectrum in the range 33 - 110 A was pre- 
sented bv IWiding and SandlinI l)l968tl . and a rocket UV 
spectrum of the Suri from 60 to 385 A was presented by 
iBehring et al.l l)l972j) . Thes e and other solar X-ray ob - 
servations were reviewed bv ICulhane and ActorJ l)l974f) . 
The launch of the second generation of solar instruments, 
the N RL spectroheliograph on Skylab and SMM led to 
lists bv lWiding et al .' fl98^ of solar fiare fines in the 170- 
345 A range, and by Phillios et aL (1982) of the flare lines 
in the 5.7 - 19 A range. iLemen et aITl|T983l published 
a compilation of observations of the iron K line from so- 
lar flares. A compilation of sol ar spectra of the K line of 
Fe"+ - Fe23+ was performed bv lSeelv et al.l l|l986lp , along 
with a comparison the oretical values calcula ted with the 
Z expansion technique. iDoschek and CowanI (|^84) com- 
piled a line list based on observations of the Sun for the 
10-200 A range. 

Motivated at least in part by solar observations, 
laboratory measurements of spectra using plasma ma- 
chines or sparks provide spectra in the EUV and X- 
ray for many io ns. Examples include s t udy o f S, A r 
and Ca ions b y Deutschman and House (1966, '1967^. 
iFeldman and Cohen (J,968) published a list of lines from 
iron m the 10-18 A range from both laboratory and so- 
lar flare spectra. iKunze et alJ l)l968l) measured He-like 
line intens ities from a theta pinch machine for C^+ and 
0^+ ions. IConnerade et al.l (^T^) compared laboratory 
plasma measurements with Hartree-Foc k structure cal- 
culations for Fe between 10 and 17 A. iFawcettI l|l970() 
classified line s betwee n 240 - 750 A produced in a laser 
plasma, and iFawcetti l)l97 j) classified lines of Fe^^"*" - 
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Fe^^ from a theta pinch machine. iGoldsmith et al.l 
lll972D classified Li - like sp ectra of K - Mn in the EUV. 
iFawcett and Haves! (119721 ) identified - and 

Co^"*" - 00^^+ 3 s — 3p and 3p — 3d lines from labo- 
ratory spectra. iFawcett et all l)l974j) made ab initio 
calculations of 2p — 3d line wavelengths and oscillator 
strengths in Fe^^+ - Fe^'^+, and compared the results 
with experimental values. Goldsmith et al. (197^ clas- 
sified 2s — 2 p transitions in l abora tory spectra of Ca^^+ 
and Tii8+. iFeldman etldl l|l97.'^ classified transitions 
of Fe"'^'^^ and Fe^^"*" obser ved in laser produced plasmas. 
IFawcett and Have^ l)l973l) classified n=2 - 3 transitions 
from 8^+ - S^^+ and Ar "+ - At^^+ . Othe r exam- 
ples include study of Fe bv IFawcett et alJ l)l972l) and by 
iBoiko et al.l (^23) ■ Comparison of theoretical and ex- 
perimental n=2 energy levels for the C, N, Li Be, O and 
F isoelectro nic sequences was carried out bv lEdlenI lll983[ 
ElsllllHl- Identifications of EUV lines of Fe^+ - Fe^^+ 
based on solar and laboratory data were performed by 
I.TTroenI lll993D . 

Recent Developments 

A commonly accepted standard set of line lists 
for all applications have been compiled from 
laboratory spectra and energy levels for vari- 
ous elements by NIST. These inclu d e com pila- 
tions for: Mg llKaufman and MartlnL Il991a^ Al 
llKaufman and Martin'. 'l991 b|; iMartin and Zalubal 
119791). Si "( Martin an d Zaluba4 119831). a nd S 
(|Kaiifman aTd Martini Il993f iMartin et all 1199(1) . 
Collections of energy levels and transition probabil i- 
ties for iron were made by ISugar and Corlis^ lll979ll . 
ICorhss and Sueail 11198211 ISugar and CorhssI ()l985j) 
and bv IShirai et al.1 ()199(I) . The data are defined as 
NIST Standard Reference Data, implying known or 
estimated accurac y. An exception is the energy level 
data compiled by iKellvl |l987) for some ions, which 
also includes some inner shell energies. Together 
with o lder d ata on spectra of light ions (|Wiese et all 
Il966l Il996ll this has been incorp orated into the 
NIST spectroscopic database (|Fuhr et al.L Il999^ at: 
^ttp://physics. nist.gov/PhysRefData/ASD/index. html 
The accuracies of the transition strengths are evaluated 
according to the NIST convention: A=3%, B=10%, 
C=25%,D=50%,E=50%. Many of the lines in the X-ray 
band have ratings of C or below. A consequence of the 
restriction to critically evaluated data is that the data 
included in this database lag behind the production of 
data available for such a compilation. This does not 
greatly affect the selection or accuracy of data for strong 
lines in well studied wavelength bands, but it does result 
in incompleteness of some line lists for use in X-ray 
astronomy. 

Recent work using the EBIT instrument at Lawrence 
Livermore National Laboratory (LLNL) represent the be- 
ginning of an effort to fill the need for line wavelengths 
which are both comprehensive and sufficiently accurate 
to be used in fitting high resolution space astronomy 
data. iBrown et al.l l)2001a|) have measured the spectrum 



of ions of Fe^^+ - F e ^^+ in the 10 - 20 A wavelength 
range. iLenson et al.1 (12203) have measure d the spectra 
of ions o f Ar in the range 20 - 50 A, and iLenson et alJ 
(2005a''b') summarize measurements for Ne-like through 
Li-like Ar, S and Si. 

Usin g many-body perturbation theory approaches. lGul 
()2005j) has calculated the level energies of excited L shell 
complexes in iron and nickel ions. Comparison with ex- 
perimental results shows that the line wavelengths are 
generally in agreement to within <10 mA. Energy lev- 
els with comparable ac curacy have been obtained by 
iKotochigova et alJ ^|200fi^ also using a many-body pertur- 
bation method. This is a significant improvement over 
that obtainable with ab initio calculations using stan- 
dard configuration interaction methods, and is adequate 
for many purposes for fitting to observed spectra using 
Chandra and XMM — Newton. 

Attempts to update and extend the NIST database 
for use in analyzing Chandra data include the lists 
for Ne, Mg, Si, S in t he 10 - 30 A wavelength range 
by iPodobedova et all l)2002|) . Lists specific to the 
X-ray band are given as part of the CHIANTI database 
(iDere et al. L ll997HLandi" and Bhat ia i2005ilWii et all 
119991) . Other compilations widely used in analyzing 
X-ray and E UV spectra a r e APEC (jSmith et al.. 200 1| 
and MEKAL (|Mewe et all ll995lJ) . These generally 
represent compilations from diverse sources, including 
both experimental and theoretical work, in order to 
be sufficiently comprehensive to be used in analyzing 
Chandra and XMM — Newton data. Effort has been 
devoted to updating these to the accuracy required for 
interpreting high resolution observations, since most 
of these compilations were begun before the launch 
of these instruments. These updates are available as 
part of the corresponding data analysis packages in 
widespread use for analyzing X-ray spectral data (cf. 
http : / /heasarc . gsfc . nasa. gov /docs/ software/lheasoft / 
ihttp:/ /cxc. harvard. edu/ciao/J . 



B. Photoionized Plasmas 

Background 

In astrophysical sources containing a strong source of 
ionizing continuum radiation and where mechanical heat- 
ing is negligible, the ionization balance in the gas will be 
determined by photoionization. Electron-ion collisions 
play a role in recombination and cooling the gas, and 
in equilibrium the temperature is determined by balanc- 
ing the heating by pho toelectrons with coUisional cool- 
ing (jTarter et alJ . ll969(l . A key distinction between pho- 
toionized and coronal plasmas is that the electron tem- 
perature in a photoionized plasma is ^ 5 — 10% of the 
ionization potential of the dominant ions. Line emission 
is due to recombination, inner shell fluorescence, or coUi- 
sional excitation of low-lying levels. Emission line equiv- 
alent widths are defined as the ratio of the line flux to the 
average continuum level in the vicinity of the line. These 
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are generally smaller in photoionized plasmas than for 
coronal plasmas, owing partly to the contribution from 
the continuum source and also because of the reduced 
energy available from electron collisions. 

The interpretation of photoionized plasmas shares with 
coronal plasmas the need for line lists. Photoionization 
by a radiation field which is non-thermal can lead to a 
situation where there is a substantial flux of ionizing pho- 
tons with energies far above the first ionization potential 
of the abundant ions. This is in contrast to the typical 
situation in coronal plasmas. It greatly increases the im- 
portance of inner shells in opacity and ionization. As a 
result, spectra can have strong fluorescence line emission 
and inner shell lines seen in absorption as the source of 
photoionization is viewed through the transmitted spec- 
trum of intervening gas. Fluorescence emission and in- 
ner shell absorption are unique X-ray signatures which 
can be observed from gas under a wide range of con- 
ditions, including neutral or near-neutral gas. Few in- 
ner shell transitions, other than satellite lines excited by 
DR, are included in the compilations of lines for coro- 
nal plasmas. Theor e tical w avelengths for K lines were 
calculated bv iHous j lll969l) using a single c onfigu ration 
Hartree-Fock method. iKaastra and Mewd l)l993tl com- 
piled wavelengths and branching r atios for m any inner 
shell lines based on energy levels o f lLotjlll96/D . A com- 
pilation of solar s pectra of the K lines of Fe"'^^+ - Fe^'^^ 
was performed bv lSeelv et alJ (|l986), along with a com- 
parison with theore tical values c alculated with the Z ex- 
pansion technique. IShirai et al.i (,,200Q.') compiled experi- 
mental energy levels for Fe ions. 
Recent Developments 

Wavelengths and oscillator strengths of inner shell 
n=2 to n=3 transi tions in iron i ons (Fe ° - Fe^''^+) have 
been calculated by iBehar et al.l l)2001al) using hullac. 
These features form a broad unresolved transition ar- 
ray (UTA) near 16 - 17 A. Calculations of transition 
probabilities for a subset of these lines, using Civ3 and 
an extensive selection of config urations in the C I ex- 
pansion was carried out bv iKisielius et al.l l)2003(l . Ex- 
tensive calculations of the UTA transitions, including 
transitions t o n=4 states, have been carried out by 
iDubau et"al] l)2003l) using AUTOLSJ. Wavelengths, A- 
values and Auger rate coefficients for K lines of all 
the ions of iron have been p ublished bv iPalmeri et all 
<2002VPalmeri ct al.' ('2003aD. lPalmeri et al.l ll2003b') and 
Bautista et al. (2003). These represent a compilation 
of experimental results where available, supplemented 
by calculations using AUTOSTRUCTURE and hfr for the 
lower ionization states. 

An additional observable which is associated primar- 
ily with photoionized gases is the photon spectrum due 
to radiative recombination (RR). The continuum radi- 
ation associated with radiative recombination (RRC) is 
emitted as a feature at the threshold energy which has 
an exponential shape, i.e. j{e) oc exp(— (e — eTh)/kTf.), 
above threshold ej/j, with a width proportional to the 
electron temperature Tg. RRC features in coronal plas- 



mas are broadened by the higher characteristic tempera- 
ture and therefore are not easily observable. There is no 
published database of experimental RRC threshold ener- 
gies as such, although energy level databases such as the 
NIST database can be used to derive them. 

Lists of resonance lines, which can be applied to spec- 
tra seen in absorption in fr ont of a strong continuum 
source , were compiled by IVerner. Barthel and Tvtleil 
1I1994D . This includes lines from all wavelengths long- 
ward of the He"*" Lyman limit at 227.838 A and all 
the ion states of all elements from hydrogen to Bis- 
muth (2'=83) and includes experimental and criti- 
cally evaluated wav elengths . At shorter wavelengths, 
IVerner. Verner and Ferland il996fl compiled experimen- 
tal and critically evaluated energy level and derived wave- 
lengths for lines originating from ground-term multiplets 
in the spectral region 1-200 A. 

Ah initio wavelengths for some K lines of medium- Z el- 
ements have been calculated bv .Behar and Netzeil 1)20021) 
using HULLAC. Comparison with Chandra spectra sug- 
gests a likely accuracy of > 20 mA. Wavelengths for the K 
line s of oxygen hav e been calculated bv IPradhan (200j) 
and iPradhan et aP l|2003tl using BPRM. Using a multi- 
platform ap proach (autostructure, hfr and BPRM), 
iGarcia e't'a l. (2005) calculated the atomic data relevant 
to the photoabsorption near the K-edge of all the oxy- 
gen ions. The accuracy of the K line wavelengths was 
estimated to be better than 20 mA by comparison to 
the available experimental data. Experimental measure- 
ments of inner shell lines include K lines_ofEV£J__^Fe^f+ 
measu red from tokamak plasmas by IBeiersdorfer et all 
()l993l), measurements of K lines of Fe^+ -Fe^^"*" by 
iDecaux et al. (,1997) using EBIT, and the EBIT measure- 
ments o f the K lin e s from 0*+ and 0^+ bv lSchmidt et all 
1)2004 . iGu et alJ l)2005|) combined EBIT measurements 
of the K lines of oxygen ions with FAC calculations of 
population kinetics to provide line wavelengths and iden- 
tifications with accuracies ranging from 5 to 20 mA . 
Wavelengths of n = 1 — 2 inner-shell transitions in the 
Be-like up to F-like ions of magnesium, aluminum, sil- 
icon, phosphorus and sulfur were measured in a laser- 
produced plasma a nd were interpreted using the Z ex- 
pansion technique l)Faenov et al.L 1199^ . iBiemont et alJ 
(,20(10,) observed the Ka lines of Ar^+-Ar-'^^+ in a plasma 
focus discharge and modeled the spectra using hfr and 
MCDF. K shell binding energies and line wavelengths can 
also be derived from Aug er electron spectroscopy. Exam- 
ples of this ar e given by lB7uch et alJ l)l979l Il987l Il992|) 
and iLin et alJ l)2001allJ) . Ex perimental techniq ues and 
results have been reviewed byHtoEerfoh] (Hiil). 



VI. DISCRETE DIAGNOSTICS 

Beyond the detection and identification of spectral fea- 
tures, physical parameters of the gas can be derived from 
the shapes and strengths of the features. Line ratios can 
be used to infer density or temperature. These require 
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the accurate measurement of a subset of hues or spec- 
tral features, and can only be applied in situations where 
the measurement is clean enough to make these features 
unambiguous. If so, they serve to allow accurate deter- 
mination of some physical quantities describing the gas 
under study often without requiring extensive numerical 
calculation or use of spectral deconvolution. 



A. Coronal Plasmas 

Use of line ratios has been extensively applied to coro- 
nal plasmas, where they have been used to infer density, 
temperature, elemental abundance, departures from ion- 
ization equilibrium, and non-Maxwellian electron veloc- 
ity distribution. Density diagnostics rely on the com- 
petition between collisions and radiative transitions for 
two or more lines with greatly differing radiative decay 
rates. Ratios of lines with differing excitation energies 
provide temperature diagnostics. Examples of the use 
of these diagnostics in the context of classical nebulae 
such as H II regions and planetary nebulae are discussed 
in detail bv lOsterbrock and Ferlandl l)2006j) . In the Sun 
the conditions in X-ray gas vary both in time and with 
position, and the use of discrete diagnostics provides a 
means to study a limited spatial region, corresponding to 
the atmospheric zone where a given set of lines is emit- 
ted most efhciently. This may provide insight which is 
easier to interpret for objects which are clearly inhomoge- 
neous, than global modeling , which involves mode l ing th e 
entire atmosphere at once llMariskai ITooS IZirini ITosl . 
Use of such diagnostics must be applied with care, how- 
ever, since there may be neighboring regions which may 
have lower emissivity but larger emission measure and 
which can thereby contaminate the diagnostic line mea- 
surements. 



1. He-like Diagnostics 

Background 

He-like ions are stable against ionization over a com- 
paratively wide range of conditions, and they emit strong 
n=l - 2 lines which are relatively free from confusion 
with other lines. The three characteristic lines are the 
resonance r (Is^ ^Sq — \s2p ^Pi), intercombination i 
(ls2 - ls2p 3p2,i), and forbidden f (Is^ ^Sq - 
ls2s ^Si), forming the "triplet" lines. The i line actually 
consists of two components corresponding to the two dif- 
ferent total angular momentum states in the upper level, 
but these are blended at the resolution of all available 
astronomical X-ray observations. The relative strengths 
of these lines are studied through conventional line ra- 
tios R=f/i and G=(f-}-i)/r. The ratio R is sensitive to 
gas density via collisions or the effects of a strong UV ra- 
diation field driving the transition ls2s ^Si — ls2p ^P2,i- 
The ratio G is sensitive to temperature owing to the fact 
that the excitation energies for the and levels are 



lower that for the ^P levels. The G ratio is also affected 
strongly by recombination, which increases the popula- 
tions of the and ^P upper levels relative to the ^P 
level. Radiative excitation of the ^P from ground can 
have the opposite effect on G. 

The import ance of the He-like lines as diagnostics was 
pointed out bv lGabriel and Jordan! l)l969tl . who discussed 
the density dependence of R and made a densi t y dete rmi- 
nation for the solar corona. iBlumenthal et af] l)l972f) dis- 
cussed the processes populating the upper levels of these 
lines and showed that R is also temperature dependent. 
The dependence of G and R on temperature and density 
for many elements, including the effects of satellite lines 
from DR into the Li-like ion, resonances and inner shell 
ionizatio n for corona l equilib rium c onditions were calcu- 
lated by IPradhan and Shull l|l98l|) . iMewe and Schriiveil 
calculated population densities of all levels with 
principal quantum number n = 2 in several He-like ions 
with Z ranging from 6 to 20. Calculation of level popula- 
tions including the effects of photoe xcitation within the 
n = 2 manifold were performed by ' Mewe and Schriiverl 
l)l978a|) . ^ ewe and Sch riivcr (1978bj) examined the ef- 
fects of time-depende nt ionization appropriate to solar 
flares. iPradhaiil l)l98,'j) showed that recombination drives 
up the value of G in a plasma which is recombination 
dominated, i.e. photoionized, owing to the fact that RR 
preferentially populates the triplet states. 
Recent Developments 

Diagnostic use of He-like line ratios depends on rate 
coefficients and cross sections populating the n=2 lev- 
els, and on the line wavelengths needed for accurate 
identifications. Other rate coefficients needed for com- 
plete modeling of He-like diagnostics include radia- 
tive rate coefficients for the non-dipo le allowed transi- 
tions, which have been calculated by iLin et al 

] II1977D . 

and DR sat ellite line intensities, w hich ha v e bee n cal- 
culated by ISafronova. Vasilvev a nd Smith' (2001) and 
others as discussed in section IVII.A.2I Theoretical 
wavelengths for He-like ions up to Z = 100 us- 
ing non-r elativis t ic var iational calculations were calcu- 
lated by iDrakd l)l988(l . Rate coefficients for 2 pho- 
ton decay of metastable He-like ions were calculated by 
iDerevianko and Johnsorj l)l997|) . R-matrix calculations 
of transition probabilities for He-like ions were presented 
by iFernlev. Taylor and Seatoiil lll987^ and have been in - 
corporated into TOPbase ijCunto and MendozaL Il992|) . 
Wavelengths for the He-like lines in solar flares have been 
measured for the ls2s '^S'l — ls2p ^Pq.2 lines in He-like 
Ne, Mg, Na and Si to within 20-30 mA using the SOHO 
SUMER instrument ^Curdt ct al., 2000). Lifetime mea- 
surements have been carried out for the excited states 
of many ions, not just for the He-like isosequ ence, and 
many of these are reviewed bv lTraber3 l)2002f) . Experi- 
ments using ion traps can measure lifetimes for forbid- 
den and intercombination transitions, and these provide 
an important check on transition probabilities for use in 
modeling He-like lines. 

The rate coefficients for coUisional excitation of 
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He-like ions were review ed by iDubad l)l994) and 
iKato and Nakazakil l)l989j) . who showed the importance 
of resonances in the ls2s manifold for high Z ele- 
ments, i.e., Ca and Fe. These have been taken into 
account in the modified re lativistic DW calculations of 
IZhang and Sampsonl (^^. Adoption of these rate co- 
efficients is recommended for high Z elements, and R- 
matrix calculations for low Z elements. These rate coeffi- 
cient s, together with radiati ve recombination coefficients 
from lPeauignot et a l.l ('1991') , calculations of /- values and 
energy levels from|Fernley, Tavlor and ScatoiJj^^^LSiii^ 
DR satellite calculations of Bclv-Dubau et alJ l)l982albl) 
were used bv lPorauet and Dubau (..2000.') to model He-like 
line ratios, including a discussion of recombining plasma 
compared with coronal plasma. iBaiitista an d Kallman' 
12000) calculated the line emissivities and level popula- 
tions using similar coUisional rate coefficients plus a re- 
combination cascade treatment, including the effects of 3- 
body recombination and suppression of DR a t high den- 
sity. The calculations bvlPorauet and Dubaul l)2000() and 
iBautista and KallmanI l)20r)ri^ provide consistent phys- 
ically reasonable values for density when applied to 
Chandra spectra fro m both recombination and co Uision- 
dominated objects. iHarra-Murnion et alJ l)l996(l calcu- 
lated rate coefficients for the emission of the lines from 
S^'*+, bench marked these against observations of toka- 
maks, and then used them to set limits on the density in 
solar flares. 



2. Dielectronic Satellite Lines 

Lines emitted during dielectronic recombination (DR) 
are in principle sensitive temperature diagnostics because 
the emission process involves collisional excitation of a re- 
combining ion, so the rate depends on the fraction of elec- 
trons capable of surmounting the core excitation energy 
barrier. DR is a two-step recombination process that be- 
gins when a free electron approaches an ion, coUisionally 
excites a bound electron of the ion, and is simultane- 
ously captured. This process is discussed in more detail 
m Section lyiTOl here we limit ourselves to discussion 
of satellite emission. Satellite lines are emitted during the 
stabilization of the doubly excited state, when the core 
electrons relax to their ground level in the presence of 
the additional recombincd electron. This spectator elec- 
tron is responsible for shifting the line wavelength away 
from the wavelength of the resonance line in the parent 
ion. Detection of satellite lines requires spectral resolu- 
tion e/Ae > 500. Interpretation of satellite intensities is 
simplified by the lack of ambiguity about the excitation 
mechanism, although satellites can also be emitted fol- 
lowing inner shell collisional or radiative excitation. Even 
if they are not resolved, satellites can contribute signif- 
icantly to the intensity of the adjacent resonance line. 
Thus they can affect the use of discrete diagnostics, such 
as the H-like diagnostics, and must be taken into account 
for accurate treatment of discrete diagnostics. 



A basic description of the calculation of intensities 
of DR satell ites and their diagnost ic value was first 
presen ted by iGabrieJ l)l972() and by iGabriel and Paeel 
l)l972j) . Most DR proceeds through high n [n > 50) 
for ions in isoelectronic sequences beyond H and He, 
and satellites from this are indistinguishable from res- 
onance lines. For satellites, only states with n < A pro- 
duce lines which can be resolved from the adjacent reso- 
nance line. The low n satellites require different assump- 
tions to calculate compared with the high-n sta tes, and 
the im portance of low-n states increases with Z. iGabriell 
II1972D showed that satellites are temperature diagnos- 
tics, largely independent of ionization balance. The satel- 
lite intensity depends on Z^, and so they are relatively 
unimportant for Z< 10. The formalism for treating satel- 
lite intensities makes the assumption of LTE between 
satellite upper level and the relevant autoionizing con- 
tinuum. The satellite intensity then depends on energy 
separation and branching ratio for line emission. The 
resonance line intensity must also be corrected for unre- 
solved satellites, and satellite intensities must include the 
effects of direct excitation of inner shells in the adjacent 
ion. Enhancement of resonance line intensities by un- 
resolv ed satellites has been treated by lAnsari and AlamI 
1)19751) . IGabriel and Phillips! |)1979D showed that satel- 
lites can be a diagnostic of non-Maxwellian electron en- 
ergy distributions. This is because of the relative im- 
portance of inner shell excitation to DR. The resonance 
line can be excited by all electrons with energies greater 
than threshold, while DR is a resonant process, selec- 
tively excited by electrons at an energy corresponding to 
the particular resonance. Thus, with 2 or more measured 
satellites plus the resonance line it is possible to deter- 
mine if the excitation rate is consistent with a Maxwellian 
electron velocity distribution, and to constrain the form 
of the distribution if not. 

Calculations of satellite intensities using non- 
relativistic Hartree-Fock orbitals were re ported in a 
series of papers by Karim and Bhalla. iBhalla et alJ 
l)l975j) calculated intensity factors for DR satellite spec- 
tra for highly charged He- like ions. This includes more 
accurate calculations of autoionization rate coefficients 
and more accurate calculations of inner shell collisional 
excitation rate coefficients than the previous work by 
IGabriell This results in smaller emissivities for 

lines arising from ls2p{^P)2s^P levels. These calcula- 
tions used non-relativistic Hartree-F ock-Slater orbitals 
in LS coupl ing with CI and exch ange. iKarim and Bhallal 
11986') and iBhalla et all ljl986bfl used this technique to 
calculate i ntensity factors for the sate llite to the Lya line 
of Ne9+. iKarim and Bhallal l)l988bj) calculated DR for 
selected hydrogen fike ions (Z=10, 14, 18, 20, 22, 26, and 
28) using the Hartrec-Fock-Slatcr orbitals. The effects 
of configuration interaction and spin-orbit coupling 
were included for n < A. Configuration average rate 
coefficients were used to calculate DR rate coefficients 
for n = 5, and 1/n^ scaling was assumed for all higher 
states. This shows that the maximum DR rate decreases 
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with Z and the position of maximum shifts to higher 
temperature with increasing Z. At lower temperatures 
DR rate coefficients for light elements are greater than 
those for heavy elements while at higher temperature the 
trend is reversed. Satellite line intensities were compared 
with calculations using Thomas-Fe r mi JD^iba^i et all 
Il98l|) and Z-expansion fBitter et al.', '1984V and ~10% 
difference is found. Karim and Bhalla (1988b) found an 
^--^40% diffe rence in the maxim um DR rate for Fe com- 
pared with lDubau et al.l l)l98l|) . poss ibly due to neglect 
of som e high nl states {2lnl', n >4). iKarim and Bhallal 
lll988cl) calculated X-ray and Auger transition rate 
coefficients from doubly excited states of He-like ions for 
Z=10, 14, 18, 20, 22, 24, 26, and 28, along with intensity 
factors of satellites originating from DR of grovm d-state 
hydro genic ions via the autoionizing states. B halla et all 
l)l988|) calculated satellite intensity factors for H-like 
Fe^^+, showing good agreement with Thomas- Fermi 
calculations for n = 2, but discrepancies for n > 2. DR 
satellite spectra for high- lying resonance states o f H-lik e 
Fe and Ni were calculated by |K^i^ and Bhallal \199^ . 
and for H-like Si, Ca, and Fe l|Karim et all 11992(1 . DR 
rate coefficient for some selected ions in the He isoelec- 
tronic seque nce were calculated bv [ Karim and Bhallal 
l|l989b() and iKarim and Bhallal l|l989a|) . These authors 
presented a comp arison with other calculations, includ- 
ing that of lNilserJ |1986.) which employs relativistic wave 
functions but which uses rate coefRcients which are 
extrapolated for levels above n = 4. They argued that 
accurate treatment of the high-n levels is important for 
ions with Z< 20. Effects of radiative cascades on H-like 
PR sat ellite spectra were studied by iKarim and Bhallal 
lll988al). in an atte mpt to address conflict ing claims 
bv iGau et all l|l980|) and iDubau et all l|l98lD . Possible 
systematic discrepancies are discussed between results 
obtained using Hartree-Slater, Thomas-Fermi and 
Z-expansion techniques. Calculations of DR satellites 
arising from levels with n >4 for He- li ke Cr and Ni 
were presented by iKarim and Bhallal l|l990ll . The 
validity of the commonly used scaling for high-n 

satellites was discussed, and it was pointed out that 
high-n satellites are likely to be unresolvable, but can 
appreciably affect total line intensity and should be 
taken into account. The dependence of DR satellite 
intensity factors on n for Isel — > 2l'nr —>■ Isnl" in 
select ed hydrogen-like ions was studied by l|Karim et all 
Il99l|). and for Is'^el ls2 pnl Is^nl in He-like ions 
bv IKarim and Bhallal l|l99l|) . 

Calculations of satellite intensities using wave func- 
tions based on a Thomas-Fer mi potential were ca rried 
out by Dubau and coworkers. iDubau et all ||1980|) cal- 
culated the emission from DR satellites to the Mg^^+ 
resonance lines including CI, intermediate coupling and 
wave functions calculated with SUPERSTRUCTURE. This 
showed that the relative importance of satellites increases 
up to n — 4, and decreases for higher n. Also in- 
cluded were the cascade contribution to the lower n satel- 
lites. The results are consistent with previous work of 



IVainshtein and Safranoval l|l978(l for wavelengths and ra- 
diative transition probabilities, but differ systematically 
for th e autoionization probabilities. iDubau and Volontd 
(|198Q!) reviewed calculations and solar observations of 
satellite line intensities. iVolonte et alJ l|l987|) calculated 
DR satellites to Ca^^+ Lya. This resolved a prob- 
lem dating from the work of iDubau et alJ l|l98lf) and 
Blanchct et al. ( 1985) in which the extrapolation of the 
contribution of the low-n satellites does not properly 
converge to the resonance line for n > 4, and may ex- 
plain why the calculated intensit y for Lyct wa s lowe r 
than observed from solar flares in iDubau et al. I lll98lD . 

-Dubau et "all lll982all calculated rate coefficients for 
production of satellite lines by DR and inner-shell exci- 
tation in Ca, as well as the production of He-like spectra 
by excitation, rad i ative a nd DR and through cascades. 
iBelv-Dubau et alJ l|l982bD calculated n=l - 2 spectra of 
Fe23+ and Fe'^^+ due to inner shell direct excitation, cas- 
cade, RR an d DR. The results were compared with so- 
lar spectra. iBelv-Dubau et l|l983|) compared results 
of these calculations with spectra from the PLT toka- 
mak. The Fe^'^+ DR rate coefficient associated with the 
Is — 2p core excitation was measured for temperatures in 
the range 0.9 - 3 keV and good agreement was obtained 
with contemporaneous calculations. 

The use of DR satellites as temperature di- 
agnostics in__tlie__so]Bi_COTona w as discussed for 
iron by iBelv-Dubau et alJ lll979al). iBelv-Dubau et alJ 
lll979bl) and for Ca bv jBelv-Dubau et all l|l982a(l . 
iDoschek and FeldmanI l|l987(l discussed temperature de- 
terminations from solar flares. These authors presented 
model temperature distributions and studied their effect 
on the satellite line strengths and on the temperatures 
which would be inferred from them based on available 
calculations. A review of the X-ray emission from so- 
lar flares was provided by .Doschek U£72) , who also dis- 
cussed temperature determinations using H/He like ra- 
tios for S, Si, Mg and satellites due to DR onto He-like 
ions. The He-like rate coefficie nts for Mg^°+ were cal- 
culated by iKeenan et alJ l|l986() , who showed that the 
Is — 3p line, when compared with the n=l~2 lines, can 
be used as a temperature diagnostic in solar flares and 
active regions. 

IPhilliDS et all ||1983|) used hfr to calculate inner shell 
spectra of Fe^^"*" - Fe^^"*" under conditions of solar flares. 
The excitation was found to be prim a rily du e to DR 
rather than direct excitation. IChenI ^86^, iNilsenI 
II1987D and INilsenI l|l988(l also carried out calculations 
of satellite spectra of H- and He-like ions using relativis- 
tic multi-configuration wave functions. Calculations of 
DR coefRcients and satellite spectra for He-like ions was 
carried out by 'Vainshtein a nd Safranoval l|l978|) and by 
Safranova et al. (2000), using the Z-expansion technique 
and including Breit-Pauli operators. The systematic de- 
pendence with Z is examined and compared with previous 
work. An extensive discussion and comparison of differ- 
ent computation al methods for sat ellite spectra for iron 
was presented bv lKato et al ] |ll997ll . 
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Calculations of satellite line wavelengths and of 
the cross section for satellite emission by DR 
have been bench r narke d by the EBIT experiment. 
iBeiersdorfer et al ] lll992^ measured the satellite spec- 
trum of Fe^'*"'" and Fe^^+ in order to measure the cross 
section for DR capture and stabilization. They find rea- 
sonably close agreement with the theoretical calculations 
of Vainshtein and Safranova ( 19781. B cly-Dubau et ^l.1 
l|l979a|) . and lBelv-Dubau et alJlfl982aj) . IGu et alj ' i|2001 
measured the EBIT spectra for the other ions of iron 
Pg20+ _ pg23+ Q.^gj. wide range of incident elec- 
tron energies, and showed that unresolved DR satel- 
lites can contribute as much as 15% to the intensities 
of strong resonance lines. Satellite spectra of He-like Fe 
an d Ni obtained from a tokamak plasma were studied 
bv ISmith et al.l l)l993tl . Experimental measurements of 
satellite lines from He-like ions of Ne a nd heavier ele- 
ments using EBIT were carrie d out bv [ W argelin et al.^ 
((2001,) . .Smith et al.. (,1996,1 and lSmith et all ((2000h) . 



Fe^°+l)Mason et all Il979l) . Fe^^+ llMason and Storev , 
Tml. Fe^^+ llBhatia and MasonL 11981 . 

1986al . Ni^^+ llBhatia and Doschek , 

19981 ). Ni^s+fBhatia and Doschel Il999|) . 



Ni^^+l Bhatia and Kastncr ll980D, and 

Ni^°+ |Bhatia et all l200l| . These include energy 
levels based on CI calculations using superstructure, 
transition probabilities and collision strengths in the 
DW approximation with the addition of a correc- 
tion calculated in the Coulomb-Bethe approximation 
l|Burgess and Sheore v. 1974) to take into account high 
partial waves. These references are only a representative 
sample of available work on such density diagnos- 
tics. Other references can be found in the extensive 
bibliographic tabulations contained in the AMBDAS 
( http: / /www-amdis. iaea.org/ AMBDAS/J and ORNL 



{http://www-cfadc.phy.ornl.gov/bibhography/search.htmll 
databases. Other relevant data^ Sr application 
to fusion plasmas is collected in the ALADDIN 



I http://www-amdis.iaea.org/ALADDIN/ 1 database. 



3. Other Diagnostics 

Ions in the Be - Ne isoelectronic sequences have ground 
terms with many levels which can be mixed collisionally 
at high density. When this occurs it opens channels for 
emission in X-ray lines whose strength is therefore a den- 
sity diagnostic. This occurs for densities greater than 
~ 10^^ cm~^ for lines in the 2s^2p^ — 2s2-p^^^ transition 
array of the ions of iron Fe"+ - Fe20+ in the in the 80 - 
140 A range (Stratto n et al.lll984() . iMauc he (200l) has 
modeled the effect of the ground term mixing on lines 
in the 10 - 20 A wavelength range from these ions. The 
lines from the ions Fe^^+, Fe^^+ and Fe^^"*" in th e 14 - 1 9 
A range have been modeled bv ICornille et a l.'' ('1994b'). 
They calculated the structure, radiative and coUisional 
rate coefficients using AUTOLSJ. These can be used to 
identify lines in the spectrum of solar active regions. The 
intensities of the Fe^''^+ DR satellites to Fe^^+ at 15 A and 
of the lines of Fe^^^ at 14.2 and 16 A are shown to be 
sensitive to temperature. 

Calculations of level populations affecting emis- 
sion in strong UV, EUV and X-ray lines demon- 
strating density dependence due to the coUisional 
mixing of levels in the ground configuration, have 
been carried out in a s eries of papers b y Bha - 
tia and coworkers: Si^ (iBhatia and LandiL l2003a^ 
Si^+(lBhatia and LandiL l2003d^■ S^+ llBhatia and Land! 
I2003cl). 3 ^+ (IBhatia and Landi' l2003b^ 

Si°+ (IBhatia et al.L [l987) . Ca^-^+ ( Bhatia and MasonL 
\\9i?,m . Ca^s-Hfehatia and MasonL Il983l) 

Fe^+ iB hatia and Doschek', jiggSbT . 
Fe^°+(|Bhat ia and Doschck, 1996; Bh atia et all [2002) . 
Fe"+( Bhatia and Dos chell il993at iBhatia et al., 199J), 
Fe^'^^ (jBhatia a nd Kastne^ 'l98(]f; 'B hatia and MasonL 

Il997|) Fe^+( Bhatia and Dosclie^^^^i]200J, 

Fe^^+ (IBhatia et all Il989bl). Fe^Q+ (IBhatia and MasonL 
ll980aL iMason and BhatiaL Il983() . 



4. Fe^«+ 

The strong lines of Fe^^+ near 15 and 17 A are among 
the most prominent in the X-ray spectra of many coro- 
nal sources. The relative strengths of the 2p^ — 2p^3s 
lines near 17 A the 2p^ — 2p^3(i lines near 15 A and the 
2p6 - 2p54d lines near 12 A are temperature sensitive ow- 
ing to the differi ng rate coef ficients fo r elect ron impact 
ionization (EH) (jRavmond a nd Smithl. Il986l) . The ratio 
of the lines 2p^ — 2p^3(i lines near 15 A and the 2p^ — 2p^3s 
lines at 17 A can als o be an indicator of recombination 
llLiedahl et allll995l) . 

Calculation of transition probabilities in intermediate 
coupli ng with CI have been carried out by lLoulergu3 
l|l97lj) . The wavelengths of the 3s-3p subordinate 
lines for severa l Ne-like ions have b een calculated by 
iKastnej ()1983|) . IBhatia and Kast'ner (1,985) discussed 
the use of the the 2/ - 2p53s forbidden fine 
as a density diagnostic. They also presented observ- 
ability diagrams which serve as a convenient overview 
of the known and unobserved lines, both for Fe^^+ 
and for other ions in the neon isoelectronic sequence. 
Cross sections for excitat i on of t he Fe^^"*" lines have 
been ca lculated by 'Mann' (1983), H aeelstein and Jund 
lll987l). IZhang and Sampson (1989,) and for Fei^+ by 
[Zhang et alj l| l989(). The l i nes have bee n modeled by 
jLpulergue and Nussbaumer (197?, '1975*1, Smith et aLj 
Jl985), Bha tia and Doschck (19 92) , Cornille et_al., (1992 
Il994al) and lPhillips et al] (|l997j) . Calcul ations of excita - 
tion for Fe^'*^ - Fe^^"*" were reviewed byHadneH (Em. 
Cross sections fo r excitation of the Fe^^"*" lines have 
been calculated bylCuuta. Deb and Msezanel ()2000|) and 
IChen. Pradhan and Eissner I ()2003() 

It has long been known that the ratios o f the 2iP^ — 
2iP^3d comp onents observed from the sun (jParkinsonL 
I1973L I1975I) differ from calculations. This was empha- 
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sized by iBeiersdorfer et alJ l)2004(l . who pointed out that 
calculations of Fe-'^^+ excitation cross sections differ sys- 
tematically from the ensemble of laboratory and astro- 
physical data. Earlier suggested explanatio ns for this dis- 
crepa ncy include the effects of casca des fiGoldstein et all 
inner shell ionization of Fe^^"*" llBautistaLl200(l l and 
optical depth effects (jBhatia an d Kastner*. '1985 1). The 
effects of cascades were further explored by L och et al.l 
l)2006j) . who point out the importance of c orrecting ex- 
perim ental measurements for polarization. ISmith et al.l 
l)l985l) have shown that the Fe"'^^+ DR satellites, such as 
the 15.226 A satellite to the Fe^^+ 2p^ ^Sq - 2p^M ^Di 
line at 15.261 A can be comparable to the intensity of 
the Feis+ A 15.01 A lines. 

The strongest lines from Fe^^+ under coronal con- 
ditions are the 2p^ "^Sq — 2p^3d^Pi resonance and 
2p^ ^S'o — 2p^3c? ^Di intercombination lines at 15.01 and 
15.26 A. The relative stren gths are not highly sensitive 
to temperature or density ijLoulergue and Nussbaumerl 
ITOTSI , and they have been observed and modeled exten- 
sively l)Fawcett et all Il979j) . The ratio of the strengths 
of these two features, as measured in the laboratory 
|Brown et al., 199^ is found to be independent of con- 
ditions and excitation process, including radiative cas- 
cades, resonance excitation, and blends with unresolved 
DR satellites. The laboratory ratio is measured to be 
3.04, and is greater than that measured in many astro- 
physical sources. KBIT measureme nts and modeling of 
the lines of ijGu et alJ. |2004 confirm the rcsuhs 

for Fei^+. 



5. Optical Depth Diagnostics 

Multiple resonance scatterings will affect conveniently 
measured ratios of lines with differing thermalization 
properties. Thermalization depends on the line optical 
depth and on the probabilities of destruction or conver- 
sion per scattering. Such ratios, whose excitation is rela- 
tively insensitive to other conditions, include the relative 
strengths of the Fe"'^^+ 15 A lines. The effects of opti- 
cal depth on the ratios of X-Ray and EUV lines o f this 
ion have been studied bv lBhatia and Kastned l)l999|) . Al- 
though optical depths have been suggested as the origin 
of ratios of the 15 A lines which are discrepant with cal- 
culations, the discrepancy now appears to be affected by 
blending and possibly with o mission of important physi - 
cal effects in the calculations ijBeiersdorfer et all l2004|) . 
Other line ratios which are affected by optical depth in- 
clude the Lyman series of hydrogenic ions, and the G 
ratio of He-like n = 1 — 2 lines. In both of these cases the 
ratio in question also depends on other factors such as 
density or temperature, and this complicates diagnostic 
measurement of optical depth using these lines. 



6. Abundance Diagnostics 

Elemental abundances cannot be directly derived from 
line intensity ratios, owing to the dependence of emis- 
sivities on temperature and density. Line equivalent 
widths have been used to derive abundances in coro- 
nal plasmas, since the continuum in many spectral re- 
gions is dominated by electron-proton bremsstrahlung, 
and so is only weakly dep endent on abundances of met- 
als l)Svlwester et all ll998(l . This requires that the tem- 
perature be derived from measurement of the continuum 
shape, for example, since both the line and continuum 
emissivities depend on temperature. Abundances can 
also be derived through full spectral fitting, as described 
in the next section. 



7. Non-equilibrium Diagnostics 

Departure from ionization equilibrium shifts the ion- 
ization balance, at a given electron temperature, away 
from the equilibrium value. This has the effect of in- 
creasing the importance of processes associated with 
adjacent ions in emitting lines from a given species, 
such as inner shell ionization (excitation-autoionization 
or direct ionization) when the ionization balance shifts 
to the low side, or recombination (DR or RR) when 
the balance s hifts to the high s i de. A n example 
was studied bv iMewe and Schriiveil l)l9781j) for He-like 
lines, who showed that non-equilibrium effects can en- 
hance the inner shell contribution from ionization of 
the Li-like species, if the temperature is high and the 
ionization is low, or the recombination from the H- 
like species if th e temperature is low a nd th e ion- 
ization is high. lOelgoetz and PradhanI l)2004|) stud- 
ied time-dependent re c ombin ation dominated plasmas. 
lOelgoetz and PradhanI l)200ll) showed that the DR satel- 
lites dominate the emission of the Fe^''+ at temperatures 
below that of maximum abundance in coUisional ioniza- 
tion equilibrium. Owing to their extreme temperature 
sensitivity, the satellites are excellent spectral diagnos- 
tics for such temperatures in photoionized, coUisional or 
hybrid plasmas, whereas the forbidden, intercombination 
and resonan ce lines of Fe^^^ are not. Similar effects occur 
with Ye^^+ l)BautistaLl2000|) . and other ions. The impor- 
tance of iron K lines as diagnostics of non-equilibrium, 
and detailed modeling of th e excitation and emis sion of 
these lines, was discussed bv lDecaux et all l)2003|) . 



8. Non-Maxwellian Diagnostics 

iGabriel and Phillip! |)1979|) first suggested the use of 
DR satellites as indicators of non-thermal electron ve- 
locities. They showed that the intensities of two satel- 
lite lines Is^nl — \s2pnl with n= 2,3 relative to the 
Pg24+ resonance line Is^ — \s2p can be described by a 
single temperature only if the emitting plasma has a 
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Maxwellian distribution of ele ctrons. This effe c t was 
detected in the labor atory bv [Bartirom o et al.l l)l985f) 
and in solar flares by ISeelv et'an ~ i^8^ The effects 
of non-Maxwellian velocity distributions on the coronal 
ioniza tion bal ance have been ex amined by Dzifcakova 
lll992ft and bvlDzifcakoval Jl99j^ . following on work by 
lOwocki and Scudded l)l 9831 . who showed that the oxygen 
ionization balance in the solar corona is more sensitive 
to non-thermal effects than is iron, and that these effects 
can cause apparent differences in temperatures inferred 
from the two elements. The effects of non-Maxwellian 
velocity distributi ons on the He-like lin e ratio s have been 
examined bv .Dzifcakova and Kulinoval l)2002|) . They pa- 
rameterize the ele ctron velocity distribution a ccording to 
the formulation of lOwocki and Sc^idd cr (1983) , and show 
that at large departures from Maxwellian the value of G 
can greatly exceed that of an equilibrium coronal plasma 
at the same value of R. This corresponds to the appear- 
ance of greater temperature in the G ratio, for a given 
value of R. 



B. Photoionized 

RRCs have an exponential shape with characteris- 
tic width equal to the electron temperature, and only 
in photoionized plasmas are they likely to be narrow 
enough to be clearly detected. This was pointed out 
by iLiedahl et alJ 1)1995(1 . and has b een used to infer 
the temperature in X-ray binaries (^ Woidowski et all 
20011). and in Seyfert 2 galaxies ( Brin kman et al L I2002t 
Kinkhabwala et all 120021 ISako et all l2000|) . This pro- 
cedure relies on the determination of the background 
continuum level, and so is subject in principle to blend- 
ing and confusion with other features. The influence of 
bremsstrahlung con tinuum emission cannot be ignored 
llMewe et all ll98(J) . Blending with lines of iron com- 
plicates use of the features from Ne, for example, but 
is less important for the RRCs from Si and S. Use of 
relative strengths of RRCs to determine abundances re- 
quires atomic data for the photoionization cross section. 
The atomic data needed for this is discussed in the next 
section. 

Metastable levels in ions of the Be isoelectronic se- 
quence can provide diagnostics of the combined effects 
of gas density and UV photoexcitation. These can be 
applied to gases which are photoionized via their effect 
on the absorption spectrum. As shown bv lKaastra et al.l 
l|2004|) . the absorption from metastable 2s2p^P in the ion 
O^^ provides a constraint on density in the range 10^ - 
10^^ cm~^. They show a marginal detection of lines aris- 
ing from this level from an active galaxy observed with 
the LETG on Chandra. 

Owing to the presence of a strong continuum source, 
spectra of photoionized sources in the X-ray band can 
be viewed in transmission, and if so exhibit primarily 
absorption due to bound-free continuum and resonance 
lines. Line strengths, as measured by the line equivalent 



width, are therefore diagnostic of the gas column density 
and the line intrinsic width. Curve of gro wth analy sis, fa- 
miliar from the analysis of stellar spectra l)Mihalasl ll978'). 
has been applied to a nalysis of Chandr a spectra of ac- 
tive galaxies by, e. g. iLee et all l|200l|) . An important 
difference between the optical/UV and X-ray case is that 
the damping parameter a, which is typically < 10~® for 
a line such as Lya of hydrogen, can be as great as ^ 1 
for X-ray lines such as Fe Ka lines owing to the typically 
very short lifetim es of line upper levels l|Liedahl l2003t 
iMasai and Ishidal |1}04) . 

Broadening of bound-free absorption features is dom- 
inated by the intrinsic properties of the transition, i.e., 
level or resonance structure near the threshold and the 
phase space above the threshold. So such features serve 
as diagnostics of the presence of certain ions or atomic 
species. Examples include the resonance structure and 
edge position in ions of oxygen. These have been calcu- 
lated foi_Of_^_^Pradh^j200m, neutral oxygen, 0+ and 
llPradhan et all 1200 3|) and for all ions of oxygen by 
iGarcia eT^Lrij2005f) . The position and shape of the neu- 
tral oxygen e dge, and the reso nance structure have been 
calculated bv iGorczvca et all |2003lj) using the R-matrix 
code package. Comparison w ith experimental measure- 
ments ( Schmi dt et all |2004() using the KBIT device at 
Lawrence L ivermore Nat i onal L ab shows that the cal- 
culations of iGarcia et alJ l|2005() are accurate to within 
~ 10 mA for 0''+ and 0^+. 0^+ - 0^+ K line wave- 
lengths were measured with an accuracy ranging from 
5 to 20 mA in an KBIT experiment and were inter- 
preted in_a^_collisional-radiative model using the FAC 
code (iGu et all 120051) . These lines were observed by 
Paerels et all l|200ll) . and have been used bv lJuett et alJ 
112664}) to infer the mean ionization balance in the in- 
terstellar gas along the lines of sight to several X-ray 
binaries. 



VII. SPECTRUM SYNTHESIS AND GLOBAL FITTING 

In many situations the most detailed information 
about the conditions in astrophysical plasmas can be 
gained by modeling the microphysical processes affect- 
ing the excitation/decay, and ionization/recombination 
in sufficient detail to synthesize the effects on the spec- 
trum from many ions and elements simultaneously. This 
entails calculating the ion fractions and level populations 
and also the emitted spectrum or opacity or both. This 
procedure, sometimes described as spectrum synthesis or 
global modeling, is necessary in cases where the observed 
spectral resolution or signal-to-noise is not sufficient to 
cleanly resolve individual spectral features to be used as 
discrete diagnostics. This is the realm of data analysis in 
X-ray astronomy prior to the gratings on Chandra and 
XMM — Newton, when the instrumental resolution was 
capable of resolving only the strongest lines in uncrowded 
spectral regions, such as the Lyman-a line from 0^+ in 
coronal sources. Much of the data obtained by Chandra 
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and XMM — Newton still is obtained with the CCD in- 
struments alone, for sources which are too faint for grat- 
ings or which are spatially extended, and even grating 
spectra or bright sources are unable to resolve or unam- 
biguously detect weak lines in crowded spectral regions. 
Global modeling is also useful when there is strong cou- 
pling expected between the ions responsible for spectral 
formation. An example is in a photoionized gas, where 
in equilibrium the temperature will couple the emission 
and absorption properties of all the ions in the gas. 



A. Coronal Plasmas 

Earl y work on coronal io nization balance was carried 
out by 'Jordan! lll96!ll Il97(t . who first included a treat- 
ment of DR and autoionization and the suppression of 
these processes at high densities. Early work on the 
cooling function of coronal plasma was presen ted b y 
ICox a nd Tucker (1969) and ICox and Daltabuit l)l97l|) . 
Other calcu lations of corona l ioni z ation and emis- 
sion i n clude iLandini and Fossil |T973). I Jain and Naraic 
1 1978 [) . lAlTen and Dupred l)l969() and lTucker and Koren 
( 197lfl . The work of Mewe and coworkers served 
to better characterize many of the rate coefficients 
needed for coronal emissivity calcu lations . These 
began wi th calculations of X-ray fMew^ Il972albl : 
iMewe et al. . X985) and EUV ( Mewe, 1975J lines from 
the so lar corona. He-like lines (jMewe and Schriiveil 

Il978a|) . exploration of abundance dependences on 

lines ([IVfe we and Groncnschild, 1981,) and continuum 
(|Gronenschild an d Mew d 119781). a n d non -equilibrium 
ionization ijCronenschild and Mewel Il982|) . Many of 
these rate coefficients and cross sections were incorpo- 
rated into the mekal code (iKaastral Il992t iMewe et all 
119851 119861 Il995alb^■ and updated to the spex code 
()Kaastra. Mewe and Nieuwenhuiizenl Il996j) . a unified 
plasma model/analysis environment. 

Current ioniz ation balance calculations in w idespread 
use are those of lArnaud and Rothenflu3 l)l985fl who cal- 
culated the ionization balance for most elements of in- 
terest to X-ray astrophysics, and Arnaud and Ravmond 
II1992D who evaluated and compiled recombination and 
collisional ionization rate coeffic ients for all ions of iron , 
and the ionization balance of iMazzotta et afl lll998l). 
Both Arnaud and Rothenflud ()l985j) and lMazzotta et al.1 
1995 ) make use of many of th e rate coefficients of 



IShuU and van Steenberd il982alhD. which in turn use 
the DR rate coefficients of 'J acobs et al.l ()l977aj) . As 
po inted out in Section I VII. A. 21 the rate coefficients 
of iJacobs et all ^ll977a^ have been shown to overes- 
timate the effect of cascades to autoionizing levels, 
^rvans et al. ( 200QD have published ionization balance 
calculations which make use of the most recent rate co- 
efficients (summarized later in this section) and which 
have now been extensively bench ma rked and fitted to 
exper iments. The CHIANTI database ()Dere et all Il997l 
contains evaluated collisional excitation and ra- 



diative data , augmented to in clude wavelengths shorter 
than 50 A I^ Landi etall 119991) and excitation by pro- 
tons l|Young et all l2003j) . and including much of the 
data appropriate to coronal plasmas reported in this 
re view. The most r e cent up date to CHIANTI, version 
5 ijLandi and Bhati 1 I2005ali also includes transitions 
to the n > 4 levels of iron ions based on calcula- 
tions using FAC, and new UV lines from C-, N-, and 
0-like isosequence s . and updated level energies from 
iLandi and Phillipgl l|2005|) . 

Other calculations of coronal ionization balance 
and spectra include studies of time-dependent ioniza- 
tion ()Shapiro and Moorel Il976t [Sutherland and Dopital 
Il993|) . exploration o f the effects of ionization balance o n 
plasma diagnostics (|Gianetti. Landi and LandiniL 12000V 
calculatio n of non-equilibrium ionization in supernova 
remnants l)Hamilton. Chevalier and Sarazirj . E^83), and 
spect r a of s upernova remnants in the adiabatic phase 
()ltohl . Il979f) . Comprehensive calculations of coro- 
nal emission spectra for the X-ray band, including 
both ionization balance and spectrum synthesis by 
iRavmond and Smith ( 1977) remain in widespread use. 
These hav e recently been updated and g reatly ex- 
panded bv iBrickhouse. Ravmond and SmithI ()1995() and 
ISmith et all ()200ll) to form the core of the APEC code, 
which is in use in analyzing high resolution astrophysical 
X-ray spectra. 



1. Electron Impact Ionization (Ell) 

Background 

Electron impact ionization (EII) is of fundamental im- 
portance to coronal plasmas, in which electrons are ener- 
gized by some mechanical agent such as shocks, acoustic 
waves or MHD dissipation. It can be divided into direct 
ionization (denoted DI), which we discuss in this sub- 
section, and other processes involving an intermediate 
state which are addressed in the succeeding subsections. 
It is challenging to compute owing to the fact that it re- 
quires a treatment of two continuum electrons in the final 
state. Therefore experimental cross sections play a key 
role in the determination of accurate rates. Measurement 
techniques include crossed beams, merged beams, traps, 
plasmas, or indirect methods. Of these, crossed beam 
measurements have been most widely applied. Storage 
ring methods have the potential to eliminate the great- 
est systematic uncertainty, namely the metastable states 
in the target beam. 

A convenient point for comparison comes from th e clas- 
sical treatment of electron scattering (ISeatonLll96l . In a 
collision between two free electrons, one initially at rest, 
the energy transferred is e = E/{1 + [RE / e^Y) , where E 
is the kinetic energy of the incident electron and R is the 
impact parameter. If this energy transfer is instead inter- 
preted as the energy available for ionization of a bound 
electron, then the cross section is determined by the max- 
imum impact parameter such that the transferred energy 
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is the ionization potential, /: 



^classical 



Ih 



1 I ira, 

E 



(18) 



This formula is not accurate enough for quantitative 
work; at low energies it overestimates the cross section 
and at large energies measured cross sections decrease 
odog{E)/E. The functional behavior derived from this 
formula, with modifications, has been used to parameter- 
ize the cross section in many tabulations of cross sections 
and collision strengths derived from both experiment and 

theory. 

' ■ ■ semi-empirical formula of iLot j l)l967|) was de- 
veloped at a time when few accurate experimental 
measurements were available. Selected experimental 
rate coefficients and isoelectronic interpolation, based 
on the compilation of e xperimental rate coefficients by 
iKieffer and Duniii l)l966f) . were fitted by hand for each 
subshell to a formula 



a(E) = C 



a HE/x) 
XE 



(19) 



where x is the ionization potential of the subshell. This 
formula captures the correct behavior in the asymptotic 
(Born) limit, while allowing semi-empirical adjustment of 
the cross section near threshold. It was found that a good 
fit to the available experimental data was obtained if a 
is approximately constant and C, is the average number 
of electrons per subshell. This was done for the lowest 
two charge states of H, He, Li, N, Ne, Na, K an d for 
neutra l Ar, Kr, Rb, Xe, Cs and Hg. Based on this, iLotd 
l|l968l) derived the energy dependent cross section and 
ionization rate coefficients for all ions of elements up to 
C a, estimat e d to b e accurate to ^30%. 

ISummer'3 l)l974l) developed calculations for direct ion- 
ization using the semi-classical Exchange Classical Im- 
pact Parameter (ECIP) method. This combines a 
classical binary treatment of close collisions, which 
gives accurate cross sections at low energies, with a 
treatment which has the pr oper asymptotic behavior 
llBurgess and Su mmerslll97(j) . This method was applied 
by Summers ( 19 7i Il974l^ to calculations of ionization 
rate coefficients. iBurgess et al.l l|l977j) evaluated these 
rates by comparison with experiment, and found better 
agreement than other calculations available at the time 
for most ions. Calculatio ns using the CBO app r oxima - 
tion, were ca r ried o u t bv iGolden and SampsonI lll97^ , 
Golden et al. lll978l) . iGolden and SampsonI l)l98(]() and 
Moores et al I ()1980|) . Born approximation cross sec- 



tions were computed for ions of Al and Na (McGu ire. 
ll977llT98^ . The importance of excitation-autoionization 
and direct ioniza ti on of Na-like ions was pointed out by 
ISampson I |)1982() . IShevelko et a l. (1983) used the CBE 
approximation to calculate the rate coefficients for EII for 
ions belonging to the isoelectronic sequences from H to 



Ca. DW calculations in LS cou pling were carried out 
for the isosequen ces: H and Li llYoungeii ir980a.). He 
llYoungeil Il980b^■ Ne llYoungeil. ll981aD. Na llYoungeil 
Il981b^ . CI ( |Youngei ^. '1982a'), and Ar ("YoungeiJ '1982b'). 
These remain in widespread use for astrophysical model- 
ing. 

DW calculations are likely to be accurate for ions 
with charge greater than a few, but in order to pro- 
vide calculations for less ionized species techniques such 
close-coupling must be used. Close-coupling calcula- 
tions for the electron-impact ionization include the time- 
de pendent close-coupling ni e thod, whic h was applied 
bvlPindzola and Robicheaiu] l)200Cl() . and lPindzola et all 



(|2000D to calculations of ionization of He, C and Ne. In 
this technique the time-dependent Schrodinger equation 
is solved for the radial wave functions. These calcula- 
tions used a configuration-averaged potential due to the 
core electrons and resulted in poor agreement with ex- 
periment for Ne, suggesting the need for a full Hartree- 
Fock treatment of the interaction with the core electrons. 
Comparison of time-independent and time-dependent 
close-c oupling methods was carried out by Badnell ct aj 
||1998|) for Na-like Mg, Al, and Si. The time-independent 
methods were R-matrix and convergent close-coupling so- 
lutions, based on a total wavefunction constructed using 
antisymmetrized products of Laguerre pseudo-orbitals 
and physical bound orbitals. General agreement was 
found betwee n the results of the methods and with the 
experiment of Peart et alJ l)l99l|) at the 10% level. Cal- 
culatio n of ionization from metastable Ne was carried 
out bylBallance et all l)2004 . and for C2+ bv lLoch et J] 
l)2005j) . using the R-matrix with pseudostates method. 
Relativistic methods have not been extensively applied to 
calculations of coUisional ionizatio n. The use of MCDF 
methods has been demo nstrat e d bv | Moores and Pindzolal 
11990) and by Pindzola et alJ \l'dS^ . 
Experimental Measurements 

Experimental measurements of coUisional ionization 
have been carried out for many ions of interest to astro- 
physics, and these are summarized in Tabled Measure- 
ments using crossed beam s include those of Broo k et al.l 
(J978) for He, C, N and O. ICrandaU et al l (1982 1 for Na- 
l ike ions of Mg, Al, an d Si. lDiserens et alJ (Il984) for Ne+ 
.iGregorv etall ()1983|) for Ne'^ + and Ar'^+.lGregorv et alJ 
(jl986) for Fe5+, Fe6+, ¥e^+ . IGregorv et all lll987D for 
Fe"+, Fei3+, and Fe^^+.lMontague et al.' (1984') for Fe+. 
Measurements of single and multiple ionization of sulfur 
atoms by EII were made bv lZiegler et al ] (I1982D . 
IStenkd 



Experim ents by 



and by 



iLinkemann et al ] II1995D illustrate the complications 
inherent in the measurement of coUisional ionization. 
These include the effects of metastable levels in the 
initial state ions, which can have greater cross section 
for coUisional ionization than the ground state, and 
excitation-autoionization (EA) in which an ion is col- 
lisionally excited to a level which autoionizes. IStenkd 
()1999() made crossed beam measurements of ionization 
of iron ions, Fe"*" - Fe^"*", showing the influence of 
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metastables as evidenced by the ionization at energies 
below the ground state threshold, and also excitation- 
autoionization on the total cross section. Metastables 
are not likely to be present in the lowest density astro- 
physical plasmas, and so must be separated from the 
total cross section for use in astrophysics. However, their 
inclusion in a self-consistent way is a goal for the real- 
istic simul a tion of finite-density astrophysical plasmas. 
iFalk et all lll988aD used crossed beams to measure the 
influence of metastables on the EII for Be- like ions of B, 
C, N, and O by varying the metastable fraction in the 
target beam, demonstrating that metastable levels can 
dominate the coUisional ionization in this isoelectronic 
sequence. 

Experimental and theoretical data have been reviewed 

by several groups^ Early reviews include those of 

Rcicffcr and Dunn l|l966f) and and van Reeemorteil 
\l97(t\. A review by the Belfast group includes 
iBeU et all lll983l) f or ele ments lighter than fluorine, 
and iLennon et all (^^3) for elements from fluorine 
to Ni. These result in a se t of r ecomm ended data 
utilizing the scaling formula of iLotd lll968l). normalized 
to the d i storted w ave calculations of lYoungeil l)l980albL 
ll981albL ll982albD and, where available, experimen- 
tal measurements. These rate coefficients ha v e bee n 
further evaluated by iKato. Masai and" ArnaudI jl99l|) . 
who find the rate coefficients too low for neutral and 
near-neutral species, and also for Na -like ions due to 
neglect of excitation-autoionization. IVoronov 
has compiled rate coefficients based on the Belfast com- 
pilation, but which attempts to correc t the se prob lems 
by adopting the rate coefficients of iLot j l)l968j) for 
many of the species in question. An extensive review 
of the physical princip les and many experimental and 
theoretical results is bv lMiilleill[T99lh . In addition, rate 
coefficients for coUisional ionization have been reviewed 
as part of all the pr eviously mentioned io n izatio n 
balance calculations: 'Arnau d and Raymond! lll992t): 
^rnaud and Rothcnflue (198 5^) and iMazzotta et al.l 
lfl998|) . Sources include lKato. Masai and ArnaudI l)l99lj) 
an d correctio ns to errors in the Belfast collections 
in |http: / / dpc.nifs.ac.jp/aladdin/ Bibliographies of 
measurements and calculatio ns of coUisional ioniza- 
tion are given by Jjurgcss and Chidichimd l)l983|) 
and by lltikawa et all l)l984|l . Iltikawal l|l99lh . and 
lltikawal i|1996() . Both theoretical cross sections 
and extensive bibliographic data have been col- 
lected at the ORNL coUisional database website, 
jhttp://cfadc.phy.ornl.gov/astro/ps/data/home.html 

The extent of the available experimental data is il- 
lustrated in Table which list the ion stage and ref- 
erence for various experimental papers appropriate to 
astrophysically abundant elements, along with the re- 
ferring compilation in cases where these have been 
adopted by a compilation such as Belfast, ORNL 
( http: / / cfadc.phy.o rnl.gov/ astro /ps /data/home. html) or 
iBdv" and van Reeemorterl ljl97Cl(l . It is apparent that 
there are multiple measurements for many ionic species. 



TABLE I Experimental measurements of coUisional ioniza- 
tion cross sections. 



Ion Reference Compilation 

H" ^Mgggwa^^m^^lgrke (1^^ Bely 

H" £gthg_gt_aL (1^^ Boly 

H" ^^^n^Brackmann (J^^ Boly 

^^^^^^^^^^^^on (1^^ Boly 

^loku^n^am Boly 

^anc^^^ (1^^ Boly 

^igon^^^^^^msco^^ ^^^^ Boly 

Ho" ^ontegu^^^ Belfast 

Ho" ^oldere^^ ^gSi) Boly 

Ho" ^^^^^^^9W) Belfast 

C" ^^^^^^ ^9W) Belfast 

C+ ^amdai^^^ (1^^ Belfast 

C+ ^Ucei^^L Burgess 

"FMken^rT lQSSal ORNL 

J^^^^^Tal^ Burgess 

^regor^e^^ ^98^ ORNL 

Crandall et al. (1979bl ORNL 

C*+ Crandall et al. (1979bl ORNL 

C^+ ^^^^^^^^998) 

N" ^^^^^^r^2Z^ Belfast 

N^+ ^^^^^^^ ^985) ORNL 

N^+ ^^er^^L (1^^ Burgess 

N^+ 'FMken!rT l983a^ ORNL 

N^+ ^^^^^^1^ ^98^ ORNL 

N*+ Crandall et al. ri979b'l ORNL 

N*+ £efrance_e^^ (1^^ 

N'^+ Crandall et al. a979al ORNL 

0~ ^ison^^^^^^msco^^ (^^^^ Boly 

0° ^^^^^^^^^^^^ Boly 

0° ^rook^^^ ^2Z^ Belfast 

0° Tite^T^T T 1968') Boly 

O^ ^^ke^^^^^^^^so^ ^^^L) Burgess 

0+ ~^Toch^r^r 72003l 

0^+ ^^^^^^1^ ^98^ ORNL 

O^"^ ^^ke^^^^^^^^so^ ^^^L) Burgess 

0^+ ~^Toch^r^r 72003l 

0^+ C^^^^^^ ^979a) ORNL 

0^+ ^^^^^^^(2^0^ 

0"^+ Falk et al. (1983a) ORNL 

0''+ ^^^^^^ ^200^ 

o^+ <2SSSi 

0^+ " "^6^7 71963) Boly 

0^+ ^ra^^^^^^98^ ORNL 

0^+ ^^^^^^^r^2£ga) ORNL 

0^+ Crandall et al. (1979bl ORNL 

0^+ Rinn et al. (19871 ORNL 

0^+ ^^^^^^^ 1,1998) 

No" ^^^^^^"(1280) Belfast 

No" ^aimiste^ ^99^ ORNL 

No+ ^^^^^^ Boly 

Ne+ ^i^eren^^^ (1^^ Belfast 

Ne^+ ^^teumo^^^L ^9^ 

Ne^+ ^^^^^^^^2g4_) Belfast 

Ne^+ ^anniste^ ^996) ORNL 

Ne^+ ^^^^^^1^ ^98^ ORNL 

Ne''+ ^armiste^ ^^2^ ORNL 

Ne''+ ^uponchol^^^ ^997) 

No^+ ^uponchoU^^L (jm^ 

No''+ " BannisteTT Tggel ORNL 



and that there is little overlap between the adopted 
datasets used by various compilations. This is a manifes- 
tation of the fact that the ionization database needs to 
be thoroughly reexamined bec ause di f ferent data bases 
give different rate coefficients lISavinL 120051) . There is 
no consensus among the widely used calculations of 
ionization balance as to which is most accurate. In 
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TABLE I continued: Experimental measurements of coUi- 
sional ionization cross sections. 



Ion 



Reference 



Ne"+ 
Ne^+ 
Ne'^+ 
Ne«+ 
Mg" 
Mg+ 
Mg+ 
Mg+ 
A1+ 
Al=+ 
Al='+ 
Al*+ 
Al=+ 
Al''+ 
Al^+ 
Si+ 
Si=+ 
Si='+ 
Si*+ 
Si^+ 

sr+ 

g4+ 

Ar° 
Ar° 
Ar+ 
Ar+ 
Ar+ 
Ar+ 
Ar=+ 
At^+ 
Ar=+ 
Ai^+ 
Ar*+ 
Ar*+ 
Ar=+ 
Ar6+ 
Ar^+ 
Ar^+ 
Ar»+ 
Fe=+ 
Fe' 



Fe^+ 
Fe^+ 
Fe"+ 



13+ 
15+ 



Fei=+ 

Nl*+ 
Ni=+ 
Nl=+ 
Ni«+ 

Ni**+ 



Ni 
Nil 



12+ 



Duponchelle e t al. 

^efranec et al, (1990,) 
Dupo nchelle et al. ('1997^ 
^ugonchcU^^^ 
JC^^^nse^^^^^c^eide^ ^^7^ 



^ecke^^^ 
Peart et al._ (1^^ 



^Uchel^e^^L 
^Uchel^e^^L 
^Uchel^e^^L 
^Uchel^e^^L 
^Uchel^e^^L 

I^ ^ric et al. 



, f 199,3'! 



(1968') 
(19941 
_ (1982) 
"2001a) 
(200la1 
(200la1 
(200la1 
^2001a) 
(19931 
(19931 
^randal^^L ^82) 
^homgsoi^^^ (139^ 
^homgsor^^^ (19941 
^e^^^i^^^^^ 
^e^hi^i^^^^^ 

'. ■ ■ ■ ■ - • 

^teghei^^^ 
^fcodruffe^^ (197^ 
^iseren^^L ^88) 
^ai^^^ QSS^ 
MWe^^^ ^85^ 
^ar^^^ ^93) 
^ateuHMt^^^ ^9^ 
^lelle^^L ^85a) 
^rego^^^L ^83) 
^randal^^^ (J^^^ 
^indzol^^^ (W8^ 
^rego^^^^ ^8^ 
^owal^^L ^8^ 
^^^^^^^^^^^son ^^9^ 
^arfia^^^ (13^ 
Zhang et al. 



(19911 
(1985al 
, (19861 
I (19861 
I (19861 
(1995bl 
._ (19871 
(19871 
(19871 

Wong et al. (19931 
^toik^^h ^95a) 
^toik^^h ^95a) 
^^^go^^^^^low^^ ^^8^ 

Stenke et al. (1995al 



^^eUe^^^^h 
^^^go^^^^^ 
^^^go^^^^^ 
^^^go^^^^^ 
Ste^c^^^^L 
^^^go^^^^^ 
^^^go^^^^^L 
^^go^^^^^ 



Ste^c^^^^L 

Wang et al. 
Wang et al. 



(19881 

(1995al 

(19881 

(19881 

(19881 

(1988) 

(1988) 



Compilation 



Belfast 

Bely 
ORNL 



ORNL 
ORNL 
ORNL 
ORNL 
ORNL 
ORNL 
ORNL 
ORNL 
Belfast 
Belfast 
Burgess 
Belfast 
Belfast 



ORNL 
ORNL 
ORNL 
ORNL 
ORNL 
ORNL 
ORNL 

ORNL 
ORNL 
ORNL 
ORNL 
ORNL 

ORNL 
ORNL 
ORNL 



ORNL 

ORNL 

ORNL 
ORNL 
ORNL 
ORNL 
ORNL 



fact. lArnaud and Rothenfliiel l)l98,'j) and lMazzotta et all 
l)l998|) primarily make use of distor ted wave calculations, 
such as those by lYoungeJ l)l980a|) , rather than experi- 
mental results directly. This is due to the remaining un- 
certainties with regard to the applicability of experimen- 
tal results to low density environments, as well as con- 



venience. Experimental resu lts such as those of IStenk3 
|)1999|) and lBannister 1 1)1996() . in which the beam compo- 
sition can be thoroughly characterized provide a hopeful 
step in this direction. 
Excitation Autoionization 

Excitation-autoionization (EA) is a process in which 
an ion is coUisionally excited to a multiply excited level 
which then autoionizes. This process can dominate the 
total coUisional ionization cross section for many ions at 
energies above the threshold for direct ionization. Its 
import ance w"as po inted out by lGoldber~ et all lll96,'i^ 
and bv iBelvl l)l968(l . An experimental de monstration of 
the importance of EA was performed by iMartin et alJ 
(11968.) . 

Calculations of th e effec t in Fe"'^^^ were carried out by 
ICowan and Mam] l)l979l) . iBurgess and Chidichimd 
ifTgS^ performed functional fitting similar to 
that of Lotz for EA. DW calculations include 
rate coefficients from ground and excit ed lev- 
els of Ar L shell ions bv ICohen et alJ lll9 981 , 
and by iGrifBn. B ottchcr and Pind zolal |l982a|) and 
iGrifBn. Pindzola and Bottcher (T987(). Also in this cat- 
egory are the calcul ations for Na-like and Mg-like ions 
and Fe^+ - F e |^+ bvlPi^dzola et all lll986alblll 9981. and 
iMitnik et all il9m . fo r Fe°+ bv iPindzola et al. (19951) 
and for all ions of Ni b v lGriffin gmd^indzolal ()1988[) and 
IPindzola et"an ljl99l|) . CBE calculations of EII with 
inclusion of E A in the H - C is oclcctronic sequences were 
performed bv ISampson andljoldeni l(l979L Il98lj) . and 
ISampsonI l)l982[l made similar calculations for Na-like 
ions 

Close-coupling calculations including inner shell ex- 
citation an d EA were made for Li-like ions of C, 
N, and O llHenrvl |1979D. Na-like ions of Al and Si 
l)Henrv and Msezane', '1982*) , and for Fe'^'^+ using R- 
matrix ( Butler and Moores^ il985j) . R-matrix with pseu- 
dostates calculations include pseudo-orbitals in order to 
allow for an accurate treatment of the bound and con- 
tinuum wavefunctions. Ions w hich have b een treated 
in this way in c lude H tBartsch at and Brav . lil996L 
([Hudson e t all Il996l) Li+ (iBrown et al.l 
(Bartschat and Bray|^ 11997^ IPindzola et al. 
iMarchalant and Br avl Il997tl . B^+ ('Marchalant et all 
ll9'97[iWoitke et al.. . lT998D. (iMitnik et al.i..l999ll and 
Na-li ke Mg, Al and SnB^neUenOllii^r/^l^+lSni 

l2nn(]D . 

An additional process which affects EII is resonance 
excitation double autoionization (REDA), which is the 
first stage of dielectronic recombination, dielectronic cap- 
ture (cL Section rVII.A.2|) . followed by double autoion- 
ization. A related process is resonant excitation auto- 
double-ionization (READI), in which the autoionization 
occurs as a single event resulting in the ejection of two 
Auger electrons. The net result is an effective ionization 
event. iLinkemann et aL I (I1995^ have studied the effects 
of EA and REDA in Fe^^"*" using the Heidelberg storage 
ring, and co mpared their r esults with distorted wave ca l- 
culations bv lBadnell et al ] (11993.) and lChen et al.l (Il99d) . 



He 

1999iL Be+ 
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EA dominates the cross section by a factor ^-^5 above 800 
eV for this ion, and REDA can contribute ^20-30% to 
the total ionization rate. Theoretical cross sections re- 
produce the magnitude of the experimental cross section, 
but they do not accurately reproduce the complex reso- 
nance structure. Possible reasons for this are the isolated 
resonance approximation which omits interacting reso- 
nance effects. Experimental meas urements of multiple 
ionization have been carried out bvlMiiller et alJ lll985bD 
for Ar + and Ar'^+, iTrair^OJ^^S^i^iEsoE^^D 
l)200lD for and Aic hele et airi|2001b|) for Ne^+. 

Semiempirical formulae for electron imp act double ion- 
ization cross sections were presented by IShevelko et al.l 
ll200,'l I200(tD . R-matrix calculations can reproduce the 
fine structure in the EA/REDA/READI cross section for 
some ions, as evidenced by the comparison shown in Fig- 
ure |S1 This shows the R-matrix calculation and crossed- 
beam measurements of cxcitation-autoionization for the 
K shell of 0^+ (,Miiller et al.-.2000,) . 



2. Dielectronic Recombination (DR) 

Background 

In coronal equilibrium ionization balance is determined 
by the relative rates of coUisional ionization and recom- 
bination, where the most important recombination pro- 
cesses at low density are radiative recombination (RR), 
discussed in the next section, and dielectronic recombi- 
nation (DR), which we discuss in this section. We offer 
first a note about terminology. It is most common to refer 
to recombination of ion j when discussing recombination 
from ion with charge state j to charge state (j — 1), and 
we will adopt this convention in most of our discussion 
of DR. In the less common case where DR is discussed 
in terms of the final state ion, (j — 1), we will refer to it 
as recombination into ion {j — I). 

DR dominates over RR in coronal equilibrium for many 
ions and at a wide range of temperature, but not uni- 
versally. The H- and He-like ions are an exception, in 
which RR dominates at the temperatures characteristic 
of equilibrium. Accurate ionization balance calculations 
require that both processes b e included. A des cription 
of the process is provided bv I Savin et al.l l)2003f) . which 
we quote with minor modification: " DR is a two-step 
recombination process that begins when a free electron 
approaches an ion, coUisionally excites a bound electron 
of the ion, and is simultaneously captured. The electron 



excitation can be labeled nL 



n'lj, , where n is the prin- 



cipal quantum number of the core electron, I is its orbital 
angular momentum, and j is its total angular momen- 
tum. This intermediate state, formed by the simulta- 
neous excitation and capture, may autoionize. The DR 
process is complete when the intermediate state emits a 
photon which reduces the total energy of the recombined 
ion to below its ionization limit. Conservation of energy 
requires that for DR to go forward i?^ = AE — Ei,. Here 
Ek is the kinetic energy of the incident electron, AE is 



the excitation energy of the initially bound electron, and 
El, is the binding energy released when the incident elec- 
tron is captured onto the excited ion. Because AE and 
Eb are quantized, DR is a resonant process.." 

For the purposes of computation, DR is generally 
treated as an independent process from RR (although 
a unified approach is discussed later in this section). 
The calculation of rate coefficients divides into two ba- 
sic parts: determination of the energy structure of the 
doubly excited levels which mediate the process, and the 
rate coefficients or branching ratios for the stabilizing 
decays. The firs t discussion of the impor tance of this 
process is that of iMassev and Bate^ ()1942|) . A more de- 
tailed sum mary of the history of the importance of DR 
is given by ISeaton and StorevI l)l976() . Historically, the 
understanding of this process has progressed at a rate 
determined primarily by the atomic structure calcula- 
tions and measurements. This is particularly true for 
DR at low temperatures, where the excitation energy of 
the initially bound core electron in the presence of the 
captured electron (i.e. the quantity AE defined above) 
is very small for each of the resonances which contribute 
to the rate coefficient. This increases the need for ac- 
curate calculations of these resonance energies, since the 
Maxwellian distribution of the recombining electron ki- 
netic energies will be sharply peaked at low energy if the 
temperature is low. An accurate calculation of the extent 
to which a given resonance overlaps with this distribution 
requires values for AE which are accurate to ^ fcTg, or 
<1 eV at Te - 10* K. 

An outline of the rate calculation was provided by 
iBates and Dakarnd l|l962|) . which we reproduce here. If 
the autoionization of the capture state is much more 
probable than stabilization, then the capture state can 
be regarded as being in LTE with the continuum to a 
good approximation. Then the DR rate is simply the 
equilibrium abundance of the capture state, given by the 
Saha-Boltzmann equation, times a (small) branching ra- 
tio expressing the fraction of ions which stabilize. That 
is, the capture and autoionization reactions will maintain 
a quasi-equilibrium between the forward and reverse re- 
actions of 



X 



(9-1) + 



(20) 



State d is an excited state of the recombined ion 

that lies above the ionization potential of ion 

and i is the ground state of the recombining ion X''^. 

The population of state d can be derived from the local 

balance between dielectronic capture and autoionization. 

The DR rate coefficient is then the density of the capture 

state times the stabilization rate: 



afiT) 



n(Xr^)+) 



(21) 



where n(X^^ ^^^): "-(Xf^) are respectively the number 
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FIG. 8 Cross section for electron-impact ionization of O^^ in the vici nity of the ls2s^ ^5 '° EA threshold. The DI contribution 
has been subtracted off. Data points are crossed beam measurements iMiiller et alll200Cfl . solid curve is R-matrix calculation. 



densities of state d and state i, rie is the electronic den- 
sity, and Fg is the damping constant for radiative stabi- 
lization. The factor 6(x|^'' ^^^) is a departure coefficient 
which reflects the fact that the true number in state d 
will be reduced from the Saha value by the branching 
ratio for autoionization, which is very close to unity: 



(22) 



In this equation Fq is the damping constant for the au- 
toionizing transition and n(K^^ ^■'^) is given by the Saha 
equation: 



(9-1)+^ 



-Si/kT 



(23) 



where and uji are respectively the statistical weights 
of states d and i, and £i is the energy difference between 
states d and i. So the rate coefficient may be written 



where Tj = Tj^ and = F"-'^ are respectively the 
lifetimes for radiative stabilization and autoionization. 
Since the lifetime for stabilization is generally much 
longer, this can be written: 



af(r)=F. 



2uj, (27rmfcT)3/2 



-Ci/kT 



(26) 



(note that modern calculations do not need to make use 
of this approximation). Fg can be expressed in terms 
of the absorption oscillator strength ftd for transition 
6 — > d, & is an excited state (or the ground state) of the 



recombined ion X 



(9-1) + 



that cannot further autoionize 



and to which the autoionizing state d can decay. 
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where 
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(28) 



and vm is the wave number of the emitted radiation. 

The DR rate for an ion is obtaine d by summing over all 
the levels h and d. iBurgessI ()l964al) pointed out the fact 
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that the summation over statistical weights and Boltz- 
mann factors can diverge formally, or can be very large 
if many states participate in the DR process. This is re- 
duced by the fact that the approximation used to derive 
Eq. H26() breaks down for large n, i.e. that the autoion- 
ization no longer dominates over stabilization. iBurgesj 
calculated the rate coefficient for recombination 
onto He^ using values for Tg calculated by extrapolating 
the coUisional excitation cross section to below threshold. 
A prescription for calculating DR rate coefficients was 
presented in the classic work of Burgess (196d), in which 
he demonstrated by numerical experiment that when the 
sum over b and d is performed on equation (|27|l the de- 
pendence on the nuclear charge and level energy separate, 
and the rate can be written as a sum over the oscillator 
strengths of the stabilizing radiative transitions and a 
polynomial in the energies of the same stabilizing transi- 
tions. These polynomials, known together as the Burgess 
General Formula (GF), allow the rate coefficient for DR 
to be calculated using an analytic function of nuclear 
charge, atomic number, and temperature. This remains 
the standard set of rate coefhcients against which others 
are compared, and which are in widespread use in various 
calculations of ionization balance. 

The GF is meant to be used for DR where the chan- 
nels which dominate the total rate are into high n lev- 
els near the DR series limit. This allows all these high 
n levels to be treated approximately as a single reso- 
nance at the energy of the series limit. It is most ac- 
curate when the temperature is comparable to the dom- 
inant core excitation, i.e. Si/kT ~ 1. It is not meant 
to be used where the dominant recombination channels 
are those at low energies far away from the series limit. 
The GF models DR as a dipole core excitation followed 
by the inverse decay. It does not allow for autoioniza- 
tion into excited states, non-dipole core excitation, or 
alternative radiative stabilization pathways by either the 
core or Rydb erg electron. Comparison with experiment 
l|SavinL ll999l) shows the GF, and modifications by Merts 
and bv lBurgess and Tworkowskil l(l976() . to be rehable to 
approximately a factor of 2. 

The importance of using accurate term energies rather 
than configuration-average energies in the calcul a tion o f 
DR rate coefficients was pointed out by IShord (^2^. 
This is because of the Boltzmann factor in Eq. H27(l . and 
also beca use of t he eff ect on the overlap of the resonance 
energies. IShord l)l969|) calculated DR rate coefficients for 
various ions of the H-like, Li-like, Na-like isoelectronic 
sequences, in addition to C+ and Ca"*" ions using mono- 
configurational screened hydrogenic wave functions in- 
cluding the effect of finite stabilization. These results 
demonstrate a qualitative difference between recombina- 
tion onto ions where An = transitions are allowed in 
the core excitation, such as Li- or Be-like ions, and those 
in which they are not, such as H-like ions. In the for- 
mer case, capture states are those with n > 100, while 
in the latter ca se capture occu rs primarily to states with 
n < 50. The iBurgesd l)l965j) formula provides results 



which are most accurate for ions dominated by the high-n 
capture states. Larger errors can result for ions where the 
low-n states dominate. For high Z the dielectronic cap- 
ture and stabilization can occur predominantly through 
states with n < A. The stabilization transition therefore 
can correspond to a spectroscopically resolvable emission 
line, although with an energy which is shifted from the 
corresponding line in the parent ion due to the partial 
screening provided by the recombining electron. These 
satellite lines have diagnostic value and are discussed in 
Section rvTX2l 

The calculations of IShord II1969D neglected exchange 
effects and channel coupling, both of which are ex- 
pected to reduce the autoionization rates. The ef- 
fects of stabilization and resonance overlap, in addi- 
tio n to exchange and chann e l coup ling, were examined 
by iBurgess and Tworkowskil l)l97(Jl using the Coulomb- 
Born approximation for recombination onto H-like ions 
with 1 < Z < 40. Comparison with the GF shows 
agreement to within 30% or better. Compilations of 
DR recombination rate cocfhcients calculated using the 
GF h ave been made by Aldrovandi a nd Peauianot (,1973l 

and remain in widespread use. 
Experiment 

Experimental treatments of DR divide into those in- 
volving plasma measurement and those involving direct 
measurement of reaction yield using a beam or trap. 
In plasma measurements ion fractions are measured or 
inferred from spectra and then the recombination rate 
is derived under assumptions abou t the ionization rate . 
This procedure was carried out bv lBreton et al.l (|1978|). 
who observed time variability of spectra from a toka- 
mak in which the time variability comes from sawtooth 
heating due to an MHD instability. By measuring the 
spectra of two ions simultaneously (Mg-like and Na-like 
Mo'^*'+ and Mo^^+), whose abundances peak near 2 keV 
in equilibrium, the io nization and DR r ate coefficients 
could be disentangled. iBrook et al.l (^2^ measured DR 
at kT = 100 eV for Fe*+-Fei°+ by adding Fe to a theta- 
pinch plasma. The measured rate c oefficients are approx- 
imately 50% of those calculated byO^HE^OD^^iZ^L, 
and a re mu ch less than those predicted by t he iBurges j 
l|l965j) GF. Ilsler. Grume and Arnuriud l)l982j) measured 
relative abundances of iron ions in a tokamak where coro- 
nal equilibrium is achieved. DR rate coefficients were 
inferred by assuming coUisional ionization rates. The re- 
sults lead to inferred DR rate coefficients whic h are ap- 
proxim ately ~10 % of those predicted by the iBurgesd 
iUHl GF. 

A related technique is the observation of DR satel- 
lites from plasma experiments . Suc h a measurement 
was carried out bv iBitter et for He-like Fe in 

toka mak. These were found to a gree w it h calculations 
of IJordaiJ l)l969() and ISummer^ J^^. iDecaux et alJ 
(|l991al) measured satellite spectra of Fe^^+ and com- 
pared with various available calculations, showing agree- 
ment in the satellite intensity factors to within ~ 10%. 
iDecaux et al.l l)l991bj) compared experiment with HF cal- 
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culations using the technique of lKarim and Bhallallll99l|) 

for n = 3 — 8 DR satelhtes of the Fe XXV Ka resonance 
hne, and discussed the diagnostic use of the DR satelhtes 
for deriving the plasma temperature. 

Measureme nts using trap s and related techniques in- 
clude those of iBriand et alJ l)l984D . who studied the en- 
ergies of the K shell resonance in Ar^^+ - Ar^^+ using 
an electron beam ion source (E BIS). Rat e s for PR onto 
He-like argoi i were measured bvlAli et alJ l)l99Cl|) also us- 
ing an EBIS.l KnapD et all l|l989l) measured DR satellites 
from recombination into Ni^^"*" using an electron beam 
ion trap (EBIT) , and inferred cros s sectio ns at the reso- 
nance energies. iBeiersdorfer et aL I (I1992^ measured DR 
satellites for He-like iron using an EBIT, along with ex- 
tensive comparison with theoretical predictions for the 
strengths of these features. Good agreement was found 
between experiment and theory for the strongest lines, 
within the 20% experimental error, although for weak 
lines much greater discrepancies were found. 

Beam experiments divide into those performed at high 
energy, using highly stripped ions, and those performed 
at energies similar to those expected in thermal equi- 
lib rium. An example of the former is the measurement 
of iTanis et al ] lll98lD . who measured the Ka radiation 
produced following electron capture by highly stripped 
§13+ _ gi6-(- ^^^^ neutral argon at 13 MeV. This fluores- 
cence process is the high energy a nalog of DR. Similar 
techniques were used by IClark et al . (1985) in a study 
of K shell excitatio n in Si"+ + He at 95 MeV, and by 
iSchultz et al ] ||1987D who measured cross sections for the 
correlated emission of two K X-rays following the colli- 
sion of S"'^^"^ ions with H2 in the energy range between 70 
and 160 MeV. 

Merged beam experiments allow reactions to be mea- 
sured at low energies, close to what is expected for ther- 
mal equilibrium. Key to merged beam measurements is 
reduction of background ions in the interaction region, 
which can be an important contaminant at low beam 
energies. Experimental measurements of DR were car- 
rie d out for C"*" recomb ining to C using merged beams 
bv iMitchell et all ljl983(l . This resulted in a lower limit 
which exceeds that calculated by Lagattuta and H ahn 
l)l982a|) . Also experimental measurement of DR for 
Mg"*" u sing merged beams was carried out bv lBelic et al.l 
|)1983D . with measured cr oss sections found to slightl y 
exceed those calcula ted bv lLaeattuta and HahiJ l)l982a|) . 
iDittner et al ] ^I1983^ carried out measurement of DR in 
Li-like C and B using merged beams, using a high en- 
ergy ion beam from a tandem VandeGraf accelerator. 
This has the advantage of low background due to the 
high beam energy, and sho wed good agreement with 
iLaeattuta and HahnI l)l982aj) in the region of the 2s — 2p 
resonance for B. In C^"*" a discrepancy with theory was 
found near threshold, in which the measurement is lower, 
possibly due to errors in the theoretical calculation of the 
2pAd and 2p4f levels. Merged beam experimental mea- 
su rement of DR i n Bor on-like N, O, F were carried out 
bv IDittner et all l)l988j) which resolved the "^D, "^P, '^S 



states of the 2s2p^ electrons. 
Theory 

As mentioned above, the GF provides DR rate coef- 
ficients which are easily applied and accurate to within 
a factor ~2 for temperatures kT ~ Si. Improvements 
to these rates divide into several categories: improved 
term energies and bound-state wavefunctions, associated 
with more accurate structure calculations; inclusion of 
other types of transitions such as autoionization into ex- 
cited states, non-dipole core excitation and stabilization 
by decay to excited levels; and examination of the com- 
putational formalism, such as the isolated resonance ap- 
proximation and the use of perturbation theory to cal- 
culate matrix elements. Improvements to the structure 
initially made use of single configuration non-relativistic 
Hartree-Fock wavefunctions in LS coupling, in order to 
study the systematic behavior of DR and related pro- 
cesses for many ions. 

Study of processes which affect DR rate coefficients, 
including Auger and radiative transition probabilities, 
and inner shell excitation, were carried out by Hahn 
and coworkers. For example, inner shell excitation in 
electron-ion collisions in which the incident electron en- 
ergy exceeds the first ionization threshold is the high en- 
ergy analog of DR. The relation bet ween this process 
and DR was explored bv iHahnl l|l977j) . who pointed out 
the importance of inner shell excitation followed by au- 
toionization. The importance of Auger ionization to the 
total ionization cros s section for highly charged ions was 
also pointed out bv iHahrJ l)l978|) . Excitation probabili- 
ties, both to discrete and continuum states, were calcu- 
lated for the inner- and outer-shell electrons, using an 
improvement to the Bethe approximation. The result 
showed a decrease in the relative transition strength to 
the continuum as the degree of ionization Zj increases. 
The branching ratios for the Auger ionization and flu- 
orescence decay were fitted as functions of Zj for ion- 
ized targets. The Auger ionization and electron fluores- 
cence cross sections were compared with the correspond- 
ing direct processes. Calculations of a related process, 
excitation followed by double autoionization of the ion 
Fe^^"*", were carried out by^jagattuta and Hahn (1981^ 
using single configura tion bound state and dis t orted wave 
continuum orbitals. iMcLauehlin and Hah^ llT98l cal- 
culated cross sections for the resonant-excitation of the 
Is electrons accompanied by the capture of an incident 
electron for the target ions Si^^+ and S"'^^+. A system- 
atic study of the dependence of the Auger and radiative 
transition probabilities of high Rydberg states on their 
pr incipal and orbi t al qu antum numbers was carried out 
bv lGau and HahiJ l)l978|) , leading to a simple empirical l- 
dependence. The scaling behavior of the transition prob- 
abilities associated with DR were studied by Rctter et a^ 
|)1978D . showing that F^ cx Z"*, F^ - const. These were 
tested against numerical computations using single con- 
figuration Hartree-Fock wavefunctions, for the Be and 
Ne isosequences. This simple scaling was found to break 
down at high Z, and a polynomial expression in Z was 
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given which is better approximation. Relativistic effects 
were shown to be important for Fc. 

Calculations of DR using single configuration non- 
relativistic Hartree-Fock wavefunctions in LS coupling 
include those of Hahn fl980j, who studies the scaling 
of a DB. with principal quantu m numbe r and with nu- 
clear charge. A correction to iBurgessI lll96,'^^ general 
formula was suggested based on these results, which 
were carried out only for the Be and Ne isosequences. 
The importance of the contributions of high Rydberg 
sta tes to DR of Ar'''+, Fe^^+ a nd Mo^^+ were studied 
by iLagattuta and HahnI l)l981a|) . leading to derivation 
of approximate Auger rate coefiicients and fluoresc ence 
yields for such states. ILagattuta and HahrJ l)l981bt) cal- 
culated the DR rate for Na-like Mo'^^+ and pointed out 
the importance of cascades to all excited levels, not just 
the ground level, which reduces the total rate, and of 
includin g doubly excited stat e s whic h increases the to- 
tal rate. ILagattuta and Hahiil l)l982bf) demonstrated the 
importance of 2p — 3d core excitations and Auger decays 
from excited states using a calculation of DR for 01^+ 
(Ne-like). These effects cause depa rtures from scal- 
ing for the transition probabilities. ILagattuta and HahrJ 
l)l982a|) calculated DR of Mg"*" and pointed out the im- 
portance of cascades to autoionizi ng levels, which leads 
to a re duction of the net rate. ILagattuta and HahnI 
l|l983a|) calculated DR for C+, for which the dominant 
excitation is 2s — 2p. The capture states are at high 
n > 100, and so are densely packed in e nergy close 
to threshold. ILagattuta and HahnI l)l983b(l calculated 
DR for Ar^''+, examining the influence of scaling 
and p ointing out the effects of Is ex citation at high en- 
ergy. iMcLaughlin and HahiJ ^ll983a^ ca lculated DR for 
C3+, and IMcLaughlin and HahnI l)l983b^ calculated DR 
for B^"*" with a n impro ved tre atment of high n states. 
iMcLaughlinand Hahn ljl983c^) calculated DR for 0^+ 
and discussed scaling of F g and F., vs. Z for vari- 
ous ions in the isosequence. ILagattuta et alJ l)l98 6^ dis- 
cussed the effect on the DR cross section of electric-field- 
induced mixing of high Rydberg state levels, fo r both 
Mg+ and Ca+ target ions. lOmar andjlalm[ j[l987t) calcu- 
lated DR for Cai2+, Ca"+. Ca^°+ . iMoussa et aLl ill 9881) 
calculated DR into Ne-like Mg^+, , Cl^+, in which 
DR involves An ^0 core excitat ion, and resul t s wer e 
compared with the expe riment of iDittner et alJ l)l983(l . 
iNasser and HahnI l)l989j) calculate DR for N'^+, 0^+ in 
order to compare with the ex periments of Dittner et al.l 
||1988|) . iRamadan and HahnI (0,989) calculated DR for 
the B-like ions of C, O, Ar and Fe. They pointed out 
the importance of accurate energies for An — 0. Owing 
to shortcomings of the single configuration Hartree-Fock 
treatment for this purpose, these authors use the HFR 
code. They estimate field effect enhancement based on 
state counting arguments, and derive an expression for 
the reduction: rp nf/{\ + Imax), and nj = 3.2 x 
10^/F-'^/'*, in which F is the field strength in V cm"^ and 
Imax is the maximum I which contributes to DR. Agree- 
ment is obtained with the experiment by iDittner et al.l 



l)l988|) using tf = 2 for C"*", but the experimental uncer- 
tainty is large. Metastables in the recombining ion can 
affect experi mental results from b eams, and these were 
calculated by 'Hahn| l)l989f) and bv lHahn and B ellantond 
(.igSQ.) . who calculated DR cross sections for metastable 
0^+ and found them to be large, ^--^10"^^ cm^ for reso- 
nances, leading to rate coefficients which are ~ 10"^ cm^ 
s~^. With some field enhancement and assuming that 
the initial ion beam is a mixture of metastable ls2s (^5 
and ^S) states , the overall fea ture of the experimental 
data of Anders en et al.l l)l992a|) was reproduced, includ- 
ing the broad peaks at the incident-electron kinetic en- 
ergies of 4.7 and 12.5 eV, and partially also at 2.5 and 
6.8 eV. Wit h a different mixture, the C "^"*" data were re- 
produced. iJaniusevic and HahnI l|l989fl calculated DR 
rate coefficients and cross sections for the 0'^+ ion where 
the 2s and 2p electrons of the initial state are excited 
to higher, n >3 states ( An ^ transitions). The 2p 
electron excitation dominates the An ^ process, while 
the 2s excitation is suppressed by the cascade corrections 
that strongly affect the intermediate states ls^2s2p'isnl 
and Is^2s2p3dnl. The 2s contribution is approximately 
15% of the An 7^ mode, while the total An 7^ con- 
tribution is roughly 10% o f the An = cross section. 
iBellantone and HahiJ |)1989D calculated DR cross sections 
and rate coefficients for the H-likc and He-like C and O 
ions. The DR cross sections for the initial metastable 
states of the He-like ions were also estimated for a few 
low-lying resonance states near the DR threshold. Com- 
parison between various electron coupling schemes was 
discussed and it was pointed out that configuration in- 
teraction (CI) w ill r educe the contrast between various 
coupling choices. iHa hn ( 1993) obtained rate formulas for 
DR by fitting all of the existing DR data for ions with 
core charges (Zc) less than 50 and the number of elec- 
trons in the target ions less than 13. A review of the 
physics of RR and DR was presented bv lHahnI ||1997|) . in 
addition pointing out more exotic modes of recombina- 
tion. These include off-shell dielectronic recombination 
(radiative DR = RDR), in which an electron capture is 
accompanied by simultaneous radiative emission and ex- 
citation of the target ion. 

Similar techniques were used by'Roszman' (1979'), who 
calculated the total rate of DR for Mo^^+ (Ne-like) us- 
ing radial orbitals obtained from central field solution 
and a Hartree-Fock calculation for ground states. The 
structure was LS term resolved, not configuration av- 
eraged. Agreement with the GF for the total rate is 
good in the temperature range 1.0-6.0 keV, but poor 
below 1.0 keV, which is attributed to the GF use of 
hydrogenic Coulomb wavefunctions for conti nuum and 
hydrogenic wavefunctions for bound states. iRoszmanl 
(il987a) calculated DR for members of the fluorine iso- 
electronic sequence, Ar^+, Fe^^+, and others using the 
single-configuration, LS-coupled, frozen-core approxima- 
tion. These calculations employed a different treatment 
for the ground state orbitals, in which the Hartree-Fock 
exchange is replaced with Cowan HXR exchange poten- 
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tial, and continuum exchange is treated using a semi- 
classical exchange potential. A significant discrepancy is 
found with the only other available calculation, that of 
iJacobs et all l|l977al) . This may be due to the overesti- 
mate o f the decays to autoionizing levels bv lJacobs et al.1 
as will be discussed below. The results also difi^er 
from the rate coefficients calculated using the GF. Simi- 
lar results were found by Roszman ( 1987b) from calcula- 
tions of PR for 0-like Ar ^°+ and Fe^*"*" using similar tech- 
niques. IKoszma ^ lll987cD calculated DR for ions of the Li 
isoelectronic sequence: Ne^"*", Ar^^"*", Fe^'^"'", and Kr^' ^+. 
When compared with iMcLauehlin and HahnI l)l983a|) a 
large discrepancy is fou nd, which is attributed to an 
incorrect assumption bv IMcLauehlin and HahnI l|l983aj) 
that all transitions from states ls3pn'l' and IsApn'V to 
ls2pn"l" are treated as stabilizing, while in fact many are 
not since, for n" > nn. they ar e not below the first ion- 
ization limit. iRoszmanI l)l989af) used the same technique 
to calcul ate rate coeffici ents for O^"*" and 0^+. C ompar- 
ison with iBadnelll l|l988|) and the lBurgessI l)l965j) general 
formula shows differences of ~40%. The effect of den- 
sity dependent correction formulas for DR as proposed 
by Jordan, Burgess, and Sho re (all unpublished) are ex- 
amined by iRoszma 3 lll989bD . along with the effects of 
metastables. Anoth er single- c onfigu ration Hartree-Fock 
calculation is that o f lYoungeil 111983^. who calculated DR 
for He-like ions of C, Al, Ar, and Fe using DW and in- 
cluding CI for a li mited set of states. The resu lts compare 
well with those of iBelv-Dubau et all l)l979aj) . 

iJacobs et al.l ()l977al) pointed out that autoionization 
can occur to another level besides the ground level of the 
recombining ion, and that this may be more probable 
than the inverse of the initial capture. This would lead 
to a greater autoionization probability and a smaller net 
recombination cross section. Also, the autoionizing level 
can be collisionally ionized as well, if the density is high, 
leading to a net reduction in the recombination cross sec- 
tion. These processes were incorpor ated into DR rate co- 
efficients f or: Fe^+ - Fe^ '^+ (Jacobs et al.lll977a^ and a ll 
ions of Si ll Jacobs et all Il9 77b'). S ( Jacobs et al].ll979|) . 
Ca, and Ni l) Jacobs et al. TT98Q) . These calculations were 
widely adopted and were later shown to be inaccurate 
owing to inclusion of autoionization into excited states 
which are energetic ally inaccessible (Badnell,, 1986a(l . For 
example, for Fe^^+ IJacobs et al.l ()l977al) included dipole 
autoionization of Is'^ 2 s3p{^P)nl ls'^2s3s{^ S)Eclc, but 
^adncU (1986a) showed that this is only energetically al- 
lowed for n > 19. At T=10^ K this translate s into a 
factor 2.5 difference in the total rate. IBadnelll lll986aD 
also showed that non-dipole autoionizing transitions are 
important for this case. 

The effects of fine structure were included in single 
configuration intermediate coupling calculations of DR 
rate coeffi c ients for excited configur ations of iron ions by 
iDasguptal l)l995|) . Comparison with ISavin et 

ID lll999afl 

shows adequate agreement, although their work leaves 
out important autoionizing channels by not explicitly 
calculating capture to states with n > 15. Scaled 



rate coefficients for O an d F-like ions were calculate d by 
iDasgupta and Whitiievl f|199Q,) and .Dasguptai ljl995j) . 

MCDF orbit a ls were used fo r DR calculations by Chen 
and coworkers. IChenI lll986c^ carried out such calcula- 
tions for the He isoelectronic sequence. Relativistic ef- 
fects drastically alter the satellite structure of Fe^''+ , but 
only affect the total rate by 20%. Relativistic effects can 
alter the rate coefficients by as much as a factor of 3 by 
altering the Auger energies for higher Z elements (e. g. 
Mo'*"+). In non-relativistic treatments and LS coupling 
the total rate is dominated by a few capture levels, while 
relativistic effects redistributes the rate to more levels. 
Calculations we re carried out for Ne-like ions {Z=18, 26 
and greater) by [Chel (|l9^i) and for F-like ions {Z=26 
and greater) bv lChenI l)l988a|) . Coster-Kronig channels 
are those in which an ion with an inner shell vacancy 
decays by autoionization of an electron with the same 
principal quantum number as the initial vacancy. The 
influence of these channels on DR rate coefficients were 
examined in the calculations of Chen (1988b). The total 
DR coefficients for , N5+, and F'^+ ions are reduced 
by 60%, 13%, and 4%, respectively, due to the inclu- 
sion of Coster-Kronig chan nels. These effects are found 
to be negligible for Z> 10. IChen and CrasemannI l)l988|) 
carried out DR calculations for Be-like ions with atomic 
numbers Z = 30, 34, 36, 42, 47, and 54. Effects of rel- 
ativity and configuration interaction on the DR satellite 
spectra and rate coefficients were studied by comparing 
the theoretical results from nonrelativistic and relativis- 
tic single-configuration and multiconfiguration Hartree- 
Fock calculations with and without Brcit interaction by 
IChenI l)l988c^) . Explicit calculations for H-hke Ne, Cr, 
Mo, and Xe show that relativistic and CI effects are im- 
portant in calculating satellite spectra. For the light ions, 
nonrelativistic calculations in intermediate coupling with 
configuration interaction may be sufficient. For medium 
heavy and heavy ions, however, ab initio relativistic cal- 
culations in intermediate coupling with CI including the 
Breit interaction are necessary. DW calculations using 
MCDF orbitals were carried out for Li- like ions bv lChenI 
II1991D . It was shown that the total rate coefficients for 
the fine-structure states for ions with Z>4 can differ by 
as much as one order of magnitude at low temperatures 
due to the effects of relativi t y and i ntermediate coupling. 
Comparison with iRoszmanI l)l987a|) . finds An=0 rate co- 
efficients 20% la rger for N e ''^ at l ow temperature, which 
is likely due to IRoszmanI l)l987a|) neglect of ls2pnl for 
I >8 and the use of non-relativistic 2s — 2p energies. For 
An >0, computed rate coefficients are lowe r at low T and 
higher at high T, than those of RoszmanI l)l987a|) . The 
discrepancy at low temperatures is likely due to discrep- 
ancy in Auger energies, and at high temperatures it is due 
to inclusion of K excitation channels in the work of Chen. 
DW calculations using M CDF orbitals have been car- 
ried o ut for B-like ions by l)Chen. Reed. Guo and Savlil 
Il998|) . Comparison with experiment ijSavin et al.lll999a|) 
shows that the MCDF rate coefficients agree with exper- 
iment to within ~ 30%. The discrepancy may be due to 
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the neglect of / > 8 Rydberg states, and possibly also 
due to overestimate of resonance energies. 

Badnell showed the importance of the intermediate 
coupling and CI effects on DR rate coefficients in a series 
of papers reporting calculations using his AUTOSTRUC- 
TURE code in Fe^^"*" and Fe^^+ . This enables the rel- 
atively rapid calculation of large numbers of radiative 
and autoionization transition rates for arbitrary atomic 
configuration s which are nee ded to calculate DR rate 
coefficients. iBadnelll l|l986al) calculated configuration- 
mixing LS coupling or intermediate coupling autoion- 
ization rate coefficients for Fe^''"'", and pointed out omis- 
sions in the s tate con t ribution in work on the same ion 
by iBelv-Dub au et alJ l)l979a|) . Calculations using the 
sa me techn i ques fo r Fe^^"*" and Fe^^"*" were carried out 
bvlBadne ll lll986bD . This was also done for Be-like ions 
llBadnellL 119871) ■ B-lik e ions llBadnell et all Il99ll) S+ - 
. ijBadndl Il99l[l and oxveen ions ijBadnell Il989l 
I1992D . Intermediate coupling calcul ations were carried 
out bv lBadnell and Pindzolal ^ll989a^ for the boron iso- 
electronic sequence, by Badnell and Pindzola flOSOb) for 
oxygen and bv Badnell and Pin dzola (X989c) for Na-like 
P and CI. The effects of coupling are expected to be less 
important for highly charged ions, i.e., charge greater 
than 20, since the change from LS coupling does not 
open significant additional radiative channels. DR for 
the ground and excite d sta t es of He-like C and O were 
calculate d by [B adnell et alJ l)l990(l . and for Li- like Al by 
iBadneH 1119901 ). These were compared with R- matrix 
calculations bv lTerao and Burkd l)l990() . and the impor- 
tance of high angular momentum states and stabiliza- 
tion were pointed out. The importance of fine struc- 
tur e and CI were exami ned in detail for He-like ions 
by iPindzola etlll l|l990(l and compared with storage 
ring experimental results. iGorczvca and Badnelll l)l996bj) 
demonstrated the importance of CI in Na-like ions, in- 
cluding Fe^^+. 

Calculations of DR cross sections are most often 
done perturbatively, using bound state wavefunctions 
which are calc u lated separately from the continuum. 
iGorczvca et all ()l996j) compared calculations using R- 
matrix with an optical potential against perturbative 
methods for calculating DR in the Ar^^+ ion, showing 
agreement between the methods. This provides valida- 
tion for the use of R-matrix methods for cases where 
the isolated resonance approximation is likely to be in- 
accurate. Rate coefficients which use continuum wave- 
functions calculated in the R-matrix method and which 
combines both RR and DR into a unified rate coefficient 
have b e en ca l culated by Nahar and Pradhan (1994) and 
iNahad (119951 Il996al) . Total recombination rate coeffi- 
cients were calculated using R-matrix for DR into Si, 
Si+, S+, S2+, C+, C, N+,_Of+ Ff+ Ne4+, Na5+, Mg6+, 
Al^+, Si^+, an d 8^°+, l|Nahail Il995(l and for iron ions 
llNahail 119961 . 

Important to any close-coupling treatment of DR is 
the effect of radiation damping. This corresponds phys- 
ically to the effect of radiative decays on the resonant 



capture state, and which formally determines the width 
of the resonance and therefore its effect on the rate co- 
efficient. The importance of this effect, and the limi- 
tat ions of R-matri x calculations of DR, were discussed 
bv IGorczvca et all |2002J, by comparing computational 
techniques applied to recombination into Fe^^+. They 
showed the importance of accurate treatment of reso- 
nance damping, adequate numerical resolution of reso- 
nances and inclusion of radiative decays to autoionizing 
states, and that these effects can be treated more effi- 
ciently and accurately using perturbative calculations. 
Recent Developments 

In a series of papers beginning with'Sav in et alJ l)l997|) . 
measurements were made of the resonance strengths and 
energies for several iron ions in the Li-Ne isoelectronic 
sequences using the heavy-ion Test Storage Ring in Hei- 
delberg, Germany. This apparatus has the combined 
advanta ges of low backgroun d and negligible metastable 
content. I Savin et all l)l999aj) measured the energy depen- 
dent cross section for Fe Fe-^^+ and Fe^^+ — > Fe^ 
An=0 DR and calculated DR rate coefficients. They 
found significant discrepancies between rate coefficients 
inferred from their measured cross sections and those 
of other published calculati ons. In a co mpariso n with 
the cal culations of Jacobs et alJ ^ll977a^ and .RoszmanI 
l)l987a|) . w ho only published M axwellian-averaged rate 
coefficients. ISavin et al.l ()l999a|) demonstrated that such 
comparison can be used to identify discrepant calcula- 
tions, but cannot clearly identify the approximation re- 
sponsible for the error, and that agreement between mea- 
sured and calculated rate coefficients is not a reliable test 
of the validity of a calculation. The accuracy of DR calcu- 
lations can best be evaluated by comparison of resonance 
strengt hs and energies in the cross section. ISavin etld] 
(|l999a') show that single configuration LS coupling cal- 
culations overestimate Fe^^+ Fe^^+ by a factor 1.6. 
This may be due to the use of LS coupling which leaves 
out important autoionizing channels, or it could be due 
to inaccurate resonance energies. 

ISavin et' al. ( 2002^ measured resonance strengths and 
energies for DR of Fe^®+ — > Fe^^"*" via n = 2 ^ n' = 
2 and n = 2 — > n' = 3 core excitations. They 
have also calculated these resonance strengths and en- 
ergies using two independent techniques: the perturba- 
tive multiconfiguration Breit-Pauli (MCBP) and multi- 
configuration Dirac-Fock (MCDF) methods, finding rea- 
sonable agreement between experimental results and the- 
oretical calculations. The left panel (a) of figure [^1 
sho ws the results of measureme nts made usi ng tokamaks 
bv llsler. Grume and Arnuriud llT982'l and IWang et a,lJ 
lll988bl). theoretic al calc ulations bylJacobs et al.l|l977a^ ■ 
'RoszmanI l)l987a|) . and iDasgupta an d Whit nevI lll994D 
and those adopted in the compilation of Ma zzotta et alJ 
(1998 ) and the RR rate coefffcients calculated by 
lArnaud and Ravmondl l)l992|) . This shows the dispersion 
in rate coefficients obtained by various workers is a fac- 
tor ~10 at low temperatures, and a factor ~3 at coronal 
temperatures, far greater than the estimated experimen- 
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tal uncertainties. The right panel (b) of Figure IHl shows 
the storage ring measurements of ISavin et al.. (,2002ai|) . 
along with c alculations using M CBP (autostructure) 
and MCDF. |Savin et al ] (|20'02b) measured the resonance 
strengths and energies for DR of Fe^^+ Fe^^+ An—Q 
core excitations. They have also calculated the DR res- 
onance strengths and energies using four different the- 
oretical techniques: autostructure, hullac, mcdf, 
and R-matrix methods. On average the theoretical reso- 
nance strengths agree to within < 10%. However, the 1 
a scatter in the comparison is 30%, so that calculations 
of individual emission lines due to DR will have a scatter 
of approximately this magnitude, although the ensemble 
will be more accurate. 

I Savin et"al] l)2003(l measured DR resonance strengths 



and energies for Fe 



,20-1- 



Fei9+ for Fe2i+ ^ Fe20+ via 



An = core excitations. They have also calculated 
these resonance strengths and energies using three tech- 
niques: multiconfiguration Breit-Pauli (MCBP) method 
using autostructure, and mcdf and FAC. Although 
there is general agreement between experiment and the- 
oretical calculations, discrepancies occur for collision en- 
ergies <3 eV. Nonetheless, the storage ring experiments 
and the careful comparison with various computational 
platforms provide crucial benchmarks. These make it 
clear which computational techniques and approxima- 
tions are most reliable, and the level of accuracy which 
can be expected. These techniques can then be applied 
on a large scale to provide DR cross sections and rate 
coefficients for astrophysical applications. 

Other stora ge rin g measurements include those of 
iDeWitt et ail lll99,'^^ for He+, IZong et alJ ^ll997^ for 



Ar^^+ IMannervik et all for 0^+, 

|1999^ for F^ +. lQlans et a,lJ ll200l fo r N^+. 



^Ign^^l 



iii.^ - -■ - Bohm et al.. 

1I2OO2I 12005^ for 0^+. iBohm et al .' (20011 for Ne^+ 
iFogle et al.l i2003D for Ki^'^+ . and iNikohc et^aP ll2004) 
for Na^+; these are summarized bv iGlansTT^dT l)2004) . 

Many of these include comparison with calculations 
which utilize relativistic many-body perturbation the- 
ory, and take into account QED effects in t he po sition 
and strengths of resonances. IKenntner et al.' f 1995) mea- 
sured DR for Li-like CI and Si and Fogle ct al. ( 2005) 
has measured DR for Be- like C, N and O using similar 
apparatus. Measurements of DR have been done using 
a single-pass merged-beam bv [Andersen et aTl l|l992alb| ) 
for He-like and Li-like ions of nitrogen, fluorine, and sil- 
icon. These authors point out the importance of reso- 
nances associated with the ls2s(^5') metastable state in 
the He-like case. They also point out the importance 
of intermediate coupling in determining the resonance 
structure and hence the DR rate coefficient in situations 
when the electron velocities sample resonances close to 
the ionization threshold for the metastable state. The 
use of merged beams for DR and for electron collisional 
excitation and for ch arge transfer has been reviewed by 
^hancuf et al.l (^2^) , and the use of storag e rings for 
measu rements of recombination is reviewed bv lSchippersI 
1I1999D. 



Building on the experimental results of Savin and 
coworkers, large scale calculations of state-selected cross 
sections and rate coefficients for DR in intermedia t e cou- 
pling using the AUTOSTRUCTURE (' Badnell et all l2003l) 
code have been carried out by Badnell, Gorczyca, and 
coworkers. These include calcula tions for th e vari- 
ous isoelectronic sequences H-Na ("Altun et al.', '2OO4 
■2005 2006; Badnell, 2006; Co lgan ct_al.,_ 2003 200^ 
boost iM itnik and Badnell, 2 0041 IZatsarTnnv et all I2OO: 
l2004ajh. .2005a.k .2006.) . In a complementary effort 
(20033) has used fac to calculate total DR rate coef- 
ficients for the H-like through Ne-like isoelectronic se- 
quences for the 7 elements Mg, Si, S, Ar, Ca, Fe and 
Ni. Together these remove the shortcomings of early cal- 
culations, namely that of LS coupling and the choice of 
energetically allowed channels for applications at coronal 
temperatures. For applications at temperatures charac- 
teristic of photoionized plasmas, large uncertainties re- 
main owing to the use of ab initio energy level structure 
intrinsic in this work. This has been shown by th e stor- 
age ring measurements of ISchippers eraP (|2004|) . who 
demonstrate the existence of strong resonances in the re- 
combination cross section for Mg®+ at electron energies 
20-70 meV. They point out that errors in the location of 
these resonances of only 100 meV can result in changes 
of a factor up to a factor ^3 in the total recombination 
rate coefficient at low temperatures. However, with these 
calculations and their experimental validation, the rate 
coefficients for DR at coronal temperatures for many ions 
of astrophysical interest have reached a level accuracy 
previously unattained. 



One remaining area of unce rtainty is the role of e xter- 
nal fields on DR. According to lBadnell et all lj2003D . this 
renders pointless efforts to compute field-free rate coeffi- 
cients to an accuracy of better than 20%. It has long been 
known that the high Rydberg states that frequently dom- 
inate t he DR process can be Stark-mixed by weak electric 
fields (iBurgess and Summ ers ,. 1969) . such as the plasma 
microfield (jJacobs. et all 119761(1 . and so increase the par- 
tial rate coefficients by factors of 2 or more over a wide 
range of n. Recently, the picture has been further compli- 
cated by the discovery that magnetic fields, when crossed 
with an electric field, strongly affect the e lectric field en- 
hancement, by reducing it in niost ca ses (jBartsch et al.l 
'1999' iRobicheaux and PindzolaLlT997h . This suppression 
of the electric field enhancement improves the applicabil- 
ity of field-free DR rates, but does not remove the un- 
certainty due to field effects in modeling of real plasmas. 
The importance of field effects, and the field strengths 
themselves, can be derived through collisional radiative 
modeling of the dynamic part of the plasma microfield 
jBadncU ct al., 2003). Previously it appeared that a 
reasonable approach would be to use the values of the 
plasma microfield for the electric field strength for use in 
the generation of field dependent data as input to plasma 
modeling, but the added sensitivity to the magnetic field 
makes this impractical. 
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FIG. 9 Left panel (a): Fe^^ "*" - Fe'^^^ rate coeffi cient for recombination. The thick sohd curve represents the experimentally 
derived rate coefficient from lSav in et al] l|2002b [). The thick e rror bars show the estima ted experime ntal uncertainty of 20%. 
Symbols show tokamak results of llsler. (>^ m^aj]j__ATnui^^ (1982') (filled c ircles) and IWang et al, (1988b) (open squares). 
Curves are DRrate coefficients calculated by Jacobs ot al. (T977a) as fitted by ShuU and van Steenbera ([l982 a b) (dotted curve), 
iBpszrnan ( 1987a) (short-dash ed curve), and of Dasgupta a nd Whit nev (1994) (long-dashed curve), and recommended PR rate 
coeffic ient of iMazzotta et al.l (1998) (dot-long-dashed curve), and recommended RR rate coefficient of lArnaud and R,avmondl 
lll992t) (thin solid curve). None of the experimental or theoretical DR rate coefficients include RR. The range of temperatures 
at which these processes are likely to occur in equilibrium are denoted for photoionized and coUisionally ionized pla smas. Right 
Panel ( b): As in (a) the thick solid curve represents the rate coefficient derived from the storage ring experiment of lSavin et alJ 
ll2002al) . and the thin solid curve shows the recommended RR rate. The dotted curve shows the MCBP calculations and the 
dashed curve shows the MCDF calculations of the DR rate coefficents calculated bv lSavin et alJ ll2002af) . Both calculations are 
for capture into states including capture states for all values of nl (i.e. Umax = oo) and both include DR via 1 — > 2, 2 ^ 2, and 
2 — > 3 c ore excitations. None of the experimental or theoretical DR rate coefficients in (a) or (b) include RR. From lSavin et all 
||2002eJ) . 



3. Collisional Excitation 

Background 

Collisional excitation or deexcitation, by electrons or 
protons, is closely associated with discrete diagnostics, 
since it is the dominant mechanism affecting the level 
populations associated with many diagnostic features. 
However, it also is key to the calculation of synthetic 
spectra and for cooling. For spectrum synthesis, compre- 
hensiveness is important in order to accurately calculate 
cooling rate coefficients and pseudo-continuum emission 
due to large arrays of blended weak lines. 

In tabulating collisional excitation data it is customary 
to w ork with the collision strength (jHebb and M enzd. 
1194(11) . rather than the cross section. The collision 
strength between a given pair of levels ilij is defined such 
that f2ij = riji, and the excitation cross section is given 

by 



where uji is the statistical weight of level i. The 
Maxwellian averaged collision strength is denoted Tij 
and preserves this symmetry. 

Although, calculations are of greater practical impor- 
tance to collisional excitation than experiments owing 
to the large number of transitions which may be of in- 
terest, measurement of plasma spectra can be used to 
infer collision rates. This makes use of assumptions 



about the excitation mechanism and with independent 
measurements of the gas density an d temperatri r e. Ex - 
amples include the measurements of iDatla et al 1 (ll97fiD . 
who deduced ionization rate coefficients for He-like B 
and C from the time histories of the lines emitted by 
these ions ions in a theta-pinch plasma and compared 
theoretical estimates using the semi-classical formula of 
Burgess with experir nental results, corrected for the ef- 
fects of metastables. iJohnston and Knnzd (ll97lD deter- 
mined electron collision-excitation rate coefficients exper- 
imentally for n=2, 3, and 4 levels of Be-like N^+, 0^+, 
Ne'^+, and Si^"^ using plasma produced in a theta pinch 
device. iTondello and McWhirteJ l)l97l[l measured exci- 
tation rate coefficients for eighteen transitions in Ne'^+. 
A review of m ethods and some results was presented by 
iKunzd (119721) . 

A general rev iew of exper i menta l techniques and re- 
sults is given bv lDunn et al ] (119951) . The use of merged 
beams for electron collisional excitation and for c harge 
transfer has been reviewed bv lPhaneuf et all l|l999j) . Re- 
cent measurements of 2s — 2p excitation in Li-like ions 
C'^ + and 0^+ have be e n carr ied ou t using merge d beams 
by iGreenwood et al. ' fl999l and 'Lo zano et all ()200ll) . 
and for bv ISrnTth et al. (,200.0_a) . Measurements of 
selected tr ansitions have provided c rucial checks on cal- 
culations (jBeiersdorfer et all l2002tl for spectra such as 
Fe^^+, and the EBIT apparatus has been used to bench- 
mark cross sections and density dependen t line ratios for 
N5+, Ar"+ and Fe2i+ (jChen et al.U2004l) . 

An early review of computational techniques for colli- 
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sional excitation is that o f'Belvl l|l966(l . These include the 
Born approximation, in addition to other computational 
methods, which require the consideration of many partial 
waves and so are slow when carried out on older comput- 
ers. For this reason the Bethe approximation was widely 
used in early work. This is based on the fact that at very 
high energy distant encounters are most important and 
the projectile remains outside the atom most of the time. 
Then the cross section can be expressed as a proportion- 
ality with the oscillator strength of the transition. The 
proportionality is expressed as a Gaunt factor, and val - 
ues of effective g were derived bv lVan Reeemorted l)l962|) . 
This allowed the use of radiative transition probabilities 
directly in modeling of coUisional plasmas. This formula 
was used, for example, in cal culation of structure and col- 
lision strengths for Fe^+ by ICzvzak and KrueeeJ l)l966j) 
and compared with and results from Hartree-Fock with 
exchange self-consistent wavefunctions. This approxima- 
tion is no longer in widespread use, since more accurate 
methods can now be conveniently applied to large num- 
bers of transitions. 

Coulomb Born Oppenheimer (CBO) calculations of 
collisional exc i tation in H-like ions were performed b y 
iBurgess etal\ (llDT Oh and i Golden an d Sampson' l)l97l|) . 
Screened hydrogenic calculations of collisional excitation 
for atoms with less t han 4 e lectrons were perfor med by 
ISampson and Parks! l)l974j) . iGolden et all l|l98l|) calcu- 
lated scaled collision stre ngths for hydrogen ic ions us- 
ing the CBO method, and lClark et al 

] 11982') calculated 

scaled collision strengths for excitation of highly charged 
ions. 

An early mention of the importance of CI for colli - 
sional excitation rate coefficients was by iLavzed l)l95ll) . 
I.TonesI l|l97(tl extended the Eissner and Nussbaumer CI 
code to include relativity, and illustrated that drastic 
changes in A values can occur depend ing on the treat- 
ment of CI. lErmolaev and .Tonej 1^22) explored the im- 
portance of CI in He-like ions and mixing between 
and "^P, especially for Z >2 using both Breit -Pauli and 
Z-expansion techniques. iNussha.umerl discussed 
the importance of CI for the ion C^+. iFlower ^ ("I 972^ es- 
timated collision strengths and showed that Li-like lines 
are useful temperature diagnostics i n solar corona if 2-2 
and 2-3 lines ca n both be measured. iBely-Diibaul l|l97,'^ 
and iFlower and Pincau dc s Foretal l)l973l) performed CI 
calculations for Fc^^+. The pitfalls of CI calculations 
which have deficient sets o f configurations and can lead 
to errors w ere examined by lNuaaba umed ljl97,'^) with the 
ion 

Fei2+ ("fI ower and Pincau des Forets, 1973] used as 
an example. .Loulergue and Nussbaumer (J972i) demon- 
strated the effects of blending on the 17.06 A line in a 
CI calculation of the structure and excitation rate co- 
efficients for Fe^^+. Jones (1974) calculated collision 
strengths for He-like ions of Si, Ca, Fe comparing the 
results of LS coupling and intermediate coupling using a 
DW method. The effects of relativity were examined by 
I Walked l)l974Jl . who calculated electron impact excitation 
of the n = I and n = 2 states of hydrogenic ions in 



the Coulomb-Born approximation using Dirac wavefunc- 
tions. 

A comparison of distorted wave and close coupling 
calculations of collision strengt hs for transitions in C^"*" 
excited by electron collisions fjFlower and Launavl Il972|) 
indicates that the distorted wave method is sufficiently 
accurate for many astrophysical and laboratory applica- 
tions. DW calculations for many ions have been car- 
ried out by Bhatia, Mason and coworkers, cited in Sec- 
tion IVII Calculations using similar techniques have 
been made for S^° + (Landi and Bhatia, 200aa) , Ca^+ 
(Land i and BhatiaLl2 003bl. Fei'^+ dCornille et allll992t) 
Fe^°+ l|Phillips et pli. 119 961. Ne2+ fLandi and Bhatia], 
I2nn5ch . Cfi^^+ llLandi a nd Bhatia, 2005b ), Ca^^ 
llLandi and Bhat ial l2005al). and A r"+ ()Eissner et aP 



l2005j) . iFawcett and MasonI l)l99l|) calculated collision 
strengths and oscillator strengths for Fe*+. Collisional 
excitation rate coefficients for H-like and He-like ions cal- 
culated using a relativistic DW met hod have been calcu- 
lated by ISampson et alJ l)l983(l and lZhang and SampsonI 
('1987'). This has also been applied to the F 
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^ isoelec- 

tronic sequences. Collections of references to cal- 
culations of co ll isiona l excitation have be e n pro - 
vided bv iKatd lll976l) iRavmond and SmithI lll977D 
i Mewe and Gronenschildl l|l98ll) . and 'Me we et alJ l|l98,'Ti 
[l^^Q)- Reviews of calculations of available electron ex- 
citation cross sections for many ions of interest were 
pre sented a s part of a 1994 conference on collisions 
bv iFossi ^nd L ant^ini (199 4). Cornille et alJ llT994a) 



iMasonI lll994l) . iPradha'nr jl9 941 . Dufto n and KingstonI 
lll994l) iBhatia et all lll994^■ i Lang and Summers^ (1994), 
i Katol lll994. ISainpson et all (19 94). Bcrrington (1994), 
iMcWhirted ]^ 9941) . iDubanil jl 994l) . and ICallawavi Il994l) . 

Recent Developments 

Close coupling calculations are available for 
an increasing fraction of transitions needed to 
model coronal plasmas. R-matrix calculations of 
collisional excitation have been carried out for 
hydrog enic ions of C I Aggarwal and Kingston . 
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interme- 



Ca jAggarwa^^n^Kingstonj, 

(|Aggarwal and KingstonI .1993^1 . R-matrix 
diate coupling calculatio ns have been made for He-lik e 
and Li-like Ar and Fe bv'Whit eford et all l|200ll 12002^. 
Calculation s for Be-like ions have been reviewed by 
iBerringt'onI (|1994|) . The breakdown of the isolated 
resonance approximation and the importance of ac- 
curate treatment of radiation damping in calculations 
of collisional excitation, in contr ast with the situation 
for DR, was pointed out by iBadnell and Pindzolal 
(1993). The point is that two or more Rydberg series 
of resonances will interfere only if they overlap because 
coupling through the background is weak, or absent. In 
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contrast, for electron- impact excitation there is strong 
coupling through the background. There is no dipole 
selection rule and consequently number of resonances 
increase greatly, so that close-coupling technique or a 
perturbative approach to overlapping resonances are 
needed. Such a situation is seldom met in photoion- 
ization and DR cross section because the number of 
accessible resonances is limited by dipole selection rules. 

A comparison of relativistic DW an d R-matrix cal- 
culations for Fe^^"*" was carried out bv iGul |2004*). In 
order to test for the importance of channel coupling 
and breadown of the isolated resonance approximation, 
the two calculations were compared using the same 
energy level structure, and good agreement between 
the two methods was found for most transitions, for 
both resonance structure and for the background cross 
section. This is in contrast to ^adncU and Pindzolal 
who 
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1993|) . who demonstrated factor ~2 errors in the 
dominant resonant contributions to the excitation of 
Mg-like ions. Thus uncertainties introduced by the 
DW + isolated resonance approximation can be sig- 
nificant. R-matrix calculation of excitation of many 
X-ray t ransitions have been carried out by the Iron 
project ijHummer et al.lll993|) . This work is continuing 
and the results are contained in the TIPTOPbase 
database ijCunto and MendozaLll992l). A sample of pub- 
licati ons include those for Fe^""*" |Butle £_aiid Z eippenL 
1200(1). Fei9+ fe utlcr a nd Zeippenl . l2001a^. Fe^^+ 
llButler and ZeippenL \2001^ . F e"+ llBinello et al, . 
1998albr) . the oxygen isos e quenc e llButler and Zeip penL 



1^^, Ca^+ llT.a,ridi et a,].L [M. Fei5+ iFissner et aU 
1999^, the CI isoseque nce (iPelan and Berringtonlll995|) . 
Fe^^+ fi ne structure llStorev et all Il996|), the B isose - 
quence l)Zhang et all 1199^ andfor Fe2+ l|ZhangL IT996|) . 

R-matrix calculations have been made for the fine struc- 
ture levels of Fe^^"*" and for Fe^"+ ( G upta and Taval 
Il998l Il999|) . Most of these do not include levels in 
the n = 4 manifold, but do account for relativistic 
effects using the Breit-Pauli Hamiltonian and inter- 
mediate coupling. Exceptions include the R-matrix 
calculations for He-like O bv iDelahave and PradhginI 
l|2002() and He-like Ne by iBautista et all l)2003(l . 
The mission of the Iron Project has been extended 
into the X-ray region by the UK RmaX Network 



I http: / / amdpp.phys.strath.ac.uk/UK_RmaX i . 

Close-coupling calculations which include n = 4 
and inter mediate coupling h a ve be en carried out for 
Fe^°+ bv iBadnell and GriffinI <)200l[l and for Fe^^+ by 
iBadncU et alJ l|200j) . and show general consistency with 
the Iron Project results. The difference between collision 
strengths calculated using the close-coupling approxima- 
tion from those calculated using DW is illustrated in Fig- 
urc^l taken from Eissncr ct al. (1999). This shows the 
effective collision strength for the 3s^ ^5*0 3s3p "^Pf 
transition in Fe^'*"'" as a function of temperature for close- 
coupling (circles) and DW calculations (squares). The 
difference at low temperatures is due to the cumulative 
effect of resonances, which are not accounted for by DW 




FIG. 10 The effective collision strength for the 3s^ ^SO 
3s3p^Pi transition in Fe^'*''' as a function of temperature for 
close-coupling (circl es) and DW calculations (squares), from 
lEissner et aiTjlQQgl) 




IcgT 

FIG. 11 Comparison of the collision strength for the 
2s^2p^(^P§) 2s^2p''(^D %) transition in Fe^°+ calc ulated 
using BPRM (solid curve; (Butler and ZeioDon", "200[f)) with 
calculated using MCDF by AgEarwal (1991) (dashed) as a 
function of temperature (from lButler and ZeiDDenl l|200fll ) 



method. At the temperature corresponding to the peak 
abundance of this ion in coronal equilibrium, log(T)~6.4, 
the two methods give very nearly the same result. The 
role of CI in coUisional calculations is illustrated in Fig- 
ure m This shows comparison of the effective colli- 
sion strength for the 2s'^2p'^{^P^) 2s'^2p^(^D^) tran- 
sition in Fe^^"*" as a function of temperature calculated 
using BPR M (Butler and Zcippcn, 2000) with those of 
lAggarwallfl991() which used MCDF. The latter, although 
it includes a more complete treatment of relativistic inter- 
actions, omits the n=3 configurations and so under pre- 
dicts the resonance structure at high temperature, where 
these resonances can be excited. 
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FIG. 12 Comparison of the collision strength for the 
3s^ ^Po 3s3p transition in A r^^ calculated usin g 
an R matrix calculation (solid curve) ijGriffin et all llflQ.^Tl . 
an isolated-resonance DW calculate (das hed curve), and ex - 
peri mental merged beam measurements iSmith et alll20o3) . 
from lSmith et al.l (l2003f) 



Measured line ratios and absolute rate coefficient are 
important for validating calculations. Absolute cross sec- 
tions have been measured for electron impact excitation 
of a few ions, for example 2s(^S-i /■?) — » 2ip(^ Pi .s /?) 
for energies near threshold bv lSavin et alJ l)l995|) (see also 
[janzeu et al. (1999) for further w ork and references t o 
othe r measurements), for Fe"'^^^ bv [ Brown et alJ l|2001b() . 
and iBeiersdorfer et alJ ll2002l |2004|) . and ~ Fe^°+ - 
Fe^3+ bv ictien et alJ l)2005j) . These serve as impor- 
tant tests for close- co upling calculations such as those 
oflCh en and PradhaiJ l)2002() . An example is shown in 
figure I12L which shows a comparison of the collision 
strength for the 3s^(^Po) 3s3p{'^D) tran sition in Ar^"*" 
calcu lated using an R matrix calculation <)Griffin et all 
an isolated-resonance DW calc ulations and experi- 
mental merged beam measurements (ISmith et all 120031) . 
These illustrate the importance of resonances, which are 
more prominent in the spin-forbidden transition shown 
here, and the difficulty in accurately calculating reso- 
nance strength. An added uncertainty in experimental 
measurements comes from the fact that experiments such 
as those shown in Figure E| must detect the scattered 
electrons in order to discriminate between various exci- 
tation channels, and the angular scattering distribution 
must be understood in order to correct for the electrons 
which miss th e detectors. Similar comparisons have been 
performed bv lDiuric et alJ l)2002j) . 

Bibliographies of the avail able data for c o Uision al ex- 
citation were p ublish e d by lltikawa et al.l l)l984l) and 
lltikawal l)l99ll Il996l |2002|) . all designed primarily 
for fu sion applications, and by iPra dhan and Gallagher 
|)1992() . Reviews organized by isosequence include 
the aforementioned conference proceedings from 1994 



(e.g. iBerringt'onl l)l994) etc.). Much recent work rel- 
evant to X-ray astronomy is contai ned in the CHIANT I 
database. These in cl ude papers b y iDere et alJ lll997l). 
iLandi et al.' ('1999'), 'Dere et al.' (2001"), 'Young et alj 
12003). Landiet al.. (.2004) .Landi and Phillips (2005^) , 
and iLandi et al.l ll2005tl. and papers describing t he 
Arcetri spectral code ^Landi and Landinilll998l 12002(1 . 



4. Radiative Transition Probabilities 

Background 

Traditional modeling of coronal plasmas assumed that 
radiative decay to the ground level is more rapid than 
excitation or decay to excited levels. As a consequence, 
all dipole-allowed line emission and cooling could be cal- 
culated without explicit reference to the rate coefficients 
for radiative decay. Recent plasma models, such as APEC 
l|Smith et al.l l200l|) . calculate all level populations ex- 
plicitly, and therefore require transition probabilities for 
transitions in all decay paths of excited levels. 

Radiative transition probabilities are conveniently pa- 
rameterized in terms of the oscillator strength, which has 
a value of unity by definition for a classical point charge. 
Non-relativistic quantum mechanica l oscillator strengths 
for hy drogen have been tabulated bv lMenzel and Pekerid 
l|l935l) . It is also customary to define the oscillator 
strength for hydrogen in terms of the semi-classical 
Kramers expression: 

^^^"'"'^ ^ ^ " ^) (^) ^'^"''"^ 

(30) 

along wi th a Gaunt factor q r (n/, n ). These are tabulated 
in, e. g. jBaker and Menzell l|l938|) . 

Classic tabulations of these in cl ude t he H-like os- 
cillator strengths of IWiese et al .' (1966'), along with 
the compilations of NIST l|Fuhr et al .. 1M9) . For 
hydrogen-lik e ions without rela t ivity exact calculation 
is possible l|Bethe and Salpeterl Il972() . Treatment of 
non- dipole allowed proce sses such as 2 photon decay 
((Shapiro andjre i^ Il959l) and magn etic dipole rate co- 
efficients~ ijParpia and JohnsonL [197^ requires treatment 
of the effects of QED and nuclear size. 

Examples of DW calculations include: tran- 
sition proba b ilities for lines from Al-like ions 
ijNussbaumeil Il977(l . transition probabilities within 
the 2s^ — 2s2p — 2p^ manifold in th e Be isosequence 
l|Muehlethaler and Nussbaumerl Il976tl. radiat iv e dat a 
for the Mg isosequence (|Christensen et all Il986fl . 
transition probabilities for ground configuration of 
S'^"'" |johnson ct al., 19Sfi), transition probabilities for 
transitions within the ground configuration of S^^ and 
other Si-like ions f Huand.ll985^ . and oscillator strengths 
for He-like allowed li nes and collision strengths among 
fine structure levels l|Zhang and Sampsonl Il987(l . It is 
worth noting that, although ground state fine structure 
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transitions often fall outside of what is typically con- 
sidered the X-ray band, the structure and transition 
probabilities are relevant to the X-ray spectrum from 
these ions seen in absorption, and they may also be 
useful for study of higher energy tran sitions which 
belong to the same isoelectronic sequence. iMartin et al.l 
|)1993|) calculated transition probabilities in the lithium 
sequence. DW calculations of oscillator strengths have 
be en carried out for 2s — 2p transitions Be, B, C, N, and 
O l)Fawcett et al.lll978|). and fo r An =0, 1 tr ansitions in 
the O. S. P llFawcettL ll986albO and C llFa^cctt. 1987) 
isosequences. 

ICowan et all l)l984|l presented theoretical calcula- 
tion of wavelengths and oscillator strengths for Fe^+. 
They pointed out the importance of spin-orbit and 
configuration-mixing effects in the calculation of wave- 
lengths and oscillator strengths for heavy-element mod- 
erately charged ions. This technique involves calcula- 
tion of the energy level structure and line spectrum for 
a variety of ions by collecting available experimental 
data from laboratory and solar measurements, adjust- 
ing in order to obtain smooth variation along isoelec- 
tronic sequence, and then using these levels to calcu- 
late the structure of model ions using a semi-empirical 
code ( Cowan. 19811) . This was done for Fe^^"*" and 
Pgi2+ 3p _ g^jj^ 3g _ 3p tra nsitions dBromage et all 
Il978a|) . Be-like and B-like iron l|Bromage et al.lll978b() . 
the — 2p3d transition array o f Fe^""*" and isoelectronic 
spectra l)Bromage and FawcettL Il977al) the 2s^2ip" — 
2£2p""Md Fe"+ and Fe^^+ Hues llBromaee et all 
1197^1), much of which was updated bv ICornihemd] 
l|l992(l . the 2p'^ — 2p^3d transition array in Fe^^+ and 
isoelectronic ions ( Bro magc and F awcfitt. 1977bl ). and 
the 3s^3)9" - 3s^3)9"-^3 d transitions in Fe^+ and Fe^^+ 
l|Bromage et al.lll977lj) . A compilation of coUisional ex- 
citation and radiative decay rate coefficients for lines of 
2s^2p'^, 252^*^+^, and 2p'^+^ configurations, and for the O, 
N, C, B, and Be isoelectron ic sequences of Ti, Cr, Fe, Ni 
was produced bv lFeldman e t al. (1983). A compilation of 
transition probabilities for Cr, Fe, Ni ions in the B, C, N, 
O, an d F isosequences was produced bv iFeldman et al.l 
lll98f]D . 

Recent Developments 

Close-coupling transition probabilities for many 
dipole-allowed transitions of interest have be e n cal- 
culated by the Opacity Project ijSeatonL Il987j) 
and are contained in the the TOPbase database 
l|Cunto and Mendoz a. 1992). These use LS coupling and 
are tabulated with theoretical wavelengths, and there- 
fore are not directly applicable to synthesis of spectra 
and do not include forbidden transitions. Their appli- 
cation to observations has been aided by the work of 
O^'uhr et al. (1999) who have matched these with experi- 
mental wavelengths where possible and redistributed the 
oscillator strength among fine structure levels. Other 
CI calculations of transition probabilities include oscil- 
lator strengths of F-like ions f or 18< Z <33 using CIVS 
l)Blackford and Hibbertl . li994|) . Iron Project calculations 



include intermediate coupling and relativistic effects in 
the Breit-Pauli approximation. These include extensive 
work on forbidden transitions within the ground configu- 
ration of complex ions, whose wavelengths generall y place 
them outside the scope of this review (Galavis et al.l 
An annotated bibliography of transition proba- 
bilities for allowed and forbidden transitions with some 
overlap to the X-ray and EUV band is pr esented by 
iBiemont and ZeippenI l)l996|) . lOuined 1)20001) calculated 
the wavelengths and oscillator strengths using hfr for 
the 3p4s, 3p4d, 3p5s and 3p5d transitions in Fe®"*"- 
Fe^^+ appearing in the soft X-ray region. Close cou- 
pling calculations of transition probabilities for fine struc- 
ture t r ansitions in S"^"*" hav e been calculated by iTavall 
(|1999'). 'St orev and ZeippenI ('2000) presented transition 
probabilities for transitions within the ground configura- 
tion of the carbon and oxygen iso-electronic sequences. 
Opacity Project line strengths have been incorporated 
into the line list for use in abso rption line studies by 
IVerner. Verner and Ferlandl ()l996l) . These used a compi- 
lation of experimental energy levels similar to that used 
by NIST and Opacity Project wavelengths derived using 
LS coupling rules for lines originating from ground-term 
multiplets in the spectral region 1-200 A . Experimental 
measurements of radiative lifetim es for excited l evels of 
many ions have been reviewed bv lTraberil (j2nn3)- 



B. Photoionized Plasmas 

In gases exposed to strong ionizing radiation, or in 
the absence of strong mechanical or non-radiative heat- 
ing, the ionization balance can be determined by the ef- 
fects of photoionization and recombination. In this case 
the gas temperature is determined by a balance between 
heating and cooling due to photon interactions. Heating 
processes include slowing down of fast photoelectrons, 
and Compton scattering, while cooling is due to emis- 
sion processes analogous to th ose in a coro n al pla sma. 
Early discussion of t his wa s bv iTarter etll\ (|1969() and 
iTarter and Salpeteil l)l969ft . Traditional application for 
photoio nization models was HII regions and planetary 
nebulae llFloweilll96ill983D . but it is now apparent that 
photoionization is dominant in many X-ray sources as 
well, such as active galaxies (e. g.. Figures H and 0J) . 

The condition for photoionization to dominate over 
collisions depends on the rate of mechanical heating, if 
any, and on the ratio of the ionizing flux to the gas den- 
sity, n, which is called the photoionization parameter. 
Various definitions are in use, those most widely quoted 
are U = N / (nc) where N is the number flux of ionizing 
photons (i.e., photons with energies greater than 1 Ry), 
and ^ = AttF/ti, where F is the energy flux of ionizing 
photons, typically in the 1-1000 Ry energy range. 

The choice of convention for ionization parameter defi- 
nition is arbitrary. Any accurate calculation of photoion- 
ization rate coefficients or heating must take into account 
the shape of the ionizing spectrum in detail; the ioniza- 
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tion parameter serves as a constant of proportionality for 
use in describing model results. Certain choices of ioniza- 
tion parameter are more directly related to the problem 
to be solved, which explains in part the origins of the 
differing conventions. Photoionization rate integrals can 
be written 



rate = / ts^aye) — 

J efh ^ 



(31) 



where eth is the ionization threshold energy for the bound 
level in question, is the ionizing flux in units erg 
cm~^ erg~^, and the photoionization cross section is 
a{e). This quantity scales with energy above threshold as 
(j{s) oc e^'^ where 7 ~ 3. So the photoionization rate and 
heating integrals are always dominated by the behavior 
at threshold if (ilog(Fe)/(ilog(e) < 7. In classical pho- 
toionized nebulae, such as HII regions or planetary neb- 
ulae, the behavior of the model depends most sensitively 
on the ionizing photon flux at 1 Ry. The total number 
of ionizing photons N is also weighted toward the flux 
at 1 Ry unless increases faster than cx e^. However, 
photons at 1 Ry have less influence on the ionization of 
ions typically observed in the X-ray band than do pho- 
tons at energies greater than ~0.25 keV. Study of X-ray 
photoionized plasmas differs from classical photoionized 
nebulae in that it is often possible to directly observe the 
photons responsible for the ionization of gas responsible 
for line emission or absorption, which is motivation for 
an ionization parameter which is more closely coupled 
to the continuum in the X- ray band tha n is N. This 
has led to the suggestion bv iNetzeil l)l996l) of the use of 
the parameter Ux = Nx / (nc) where Nx is the number 
flux of ionizing photons above 0.1 keV. The choice of ^ is 
motivated by the fact that it is more heavily weighted to- 
ward the X-ray band by using ionizing energy flux rather 
than number flux. Also, heating by Compton scattering 
is proportional to the photon energy flux. ^ has units 
implied (typically erg s cm^^), while U is dimensionless. 
The study of photoionized gases in which the pressure is 
prescribed rather than the density has led to the further 
definition of an ioniza tion parameter S = F/(cP), where 
P is the gas pressure (jKrolik et alJ . ll98H) . This quantity 
is dimensionless. It is important to note that no mat- 
ter what definition of ionization parameter is adopted, it 
serves only as a convenient scaling quantity. Real calcu- 
lations take into account the full frequency dependence of 
the radiation spectrum, and so do not depend ultimately 
on the convention which is used. 

Since the ionizing radiation spectrum may not, in gen- 
eral, be a pure thermal spectrum, photoionization is not 
conveniently parameterized by a rate coefficient depend- 
ing on temperature. Rather, it is customary to con- 
sider the full energy dependent cross section, and per- 
form the calculation of the ionization rate by integrat- 
ing this over the photon flux distribution. Photoion- 
ization cross sections are also important for calcula- 
tions of the inverse process, RR, via the Milne relation 



ijOsterbrock and FerlandLl2006fl . In the case of RR, cross 
sections to excited levels are important. In many cases of 
astrophysical interest the equilibrium populations of ex- 
cited levels are negligible, so that for the purpose of calcu- 
lating the ionization rates, it is only necessary to consider 
photoionization from the ground level or term. For this 
reason many traditional calculations were applied only 
to the ground term. Recent calculations have addressed 
rate coefficients for excited levels, which are increasingly 
being self-consistently included in level population cal- 
culations for most ions. An added feature of photoion- 
ization is that, owing to the non-thermal nature of the 
radiation field, ionization from inner electronic shells can 
be important for ions with 3 or more electrons. In con- 
trast, EII generally results in a smaller contribution from 
inner shell ionization. A thermal electron energy distri- 
bution has too few electrons at the high energies needed 
for inner shell ionization (at least 5-10 times the valence 
threshold), although excitation-autoionization is impor- 
tant for many isoelectronic sequences. 



1. Photoionization 

Background 

Analytic expressions for the photoionization cross sec- 
tion of hydrogen can be obtained in the same way as for 
bound-bound transition oscillator strengths if the upper 
level is allowed to have an imaginary principle quantum 
number, ik. The oscillator strength can then be written 
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gii{n,k) (32) 



where qrr(n,k) i s the bound- free Gaunt facto r 
(jKarzas and Latteil Il96lt iMenzel and PekerisL Il935l) . 

The cross section involves the density of states, and if 
this is calculated in the limit of large n then a reason- 
ably accurate approximation is 



(2.815 X 10^^)^4# 



(33) 



This is a hydrogenic approximation and is independent of 
angular momentum. The Gaunt factor can take angular 
momentum into account; the assymptot ic dependence of 
the c ross section on energy is oc e"'"'^/^ (jFano fc Coooeii 
119681) . The exact a nalytic non-relativist i c cros s section 
for Z=0,1 is given bv lBethe and Saloeted l)l972|) . 

A widely used simple param eterizat i on of p hotoioniza- 
tion cross sections is that of 'SeatorJ (^^^. This has 
been fitted to monoconfigurational Hartree-Fock results 
for the mo derately ionized ions i mport a nt to planetar y 
nebu lae bv iHenrv and Williams! lll968D. iHenrvl t97(t) 
Chapman and Henr^ ~ ifl97lh and IChaDman'^Liid Henrvl 
|)1972() . Photoioni zation cross s ections for He- like ions 
were calculated bv lBrownll|l97l[l . Photoionization cross 
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sections using Har t ree-Sl ater wavefunctions were calcu- 
lated by iBarfieldl lll979D . who derived a simple scal- 
ing behavior along isonuclear sequences which could be 
used to relate ionic cross se c tions to those for neutrals. 
iReilman and MansonI l)l978l Il979f) performed the first 
large scale calculations of photoionization cross sections 
using Hartree-Slater wavefunctions, resulting in tabula- 
tions of cross sections for all ions of the elements with Z 
up to 30, including inner shells. Cross sections for excited 
levels of many ions based on a hydroge nic approximation 
were calculated bv lClark et al~l l)l986^ . 

The behavior of photoionization cross sections near the 
thresholds in multi-electron ions is affected by resonance 
structure due to quasi-bound states. Accurate models 
for these require calculations using the close-coupling ap- 
proximation. Near inner shell edges these take the form 
of photo-excited core (PEC) resonances below the thresh- 
old, and the net effect is to fill in the cross section below 
threshold and smooth the inner shell edge. 
Recent Developments 

A large scale compilation of close co upling calcula - 
tions is included in the Opacity Project l)Seatonl Il987(l . 
with results contained in the TOPbase, which is in 
turn part of TIPTOPbase fCunto and Men dozaL Il992|) 
http://cdsweb.u-strasbg.fr/topbasc/home.html Also in- 
cluded as part of TIPTOPbase is the Iron Project 
database and the OP server, which can calculate Rosse- 
land mean opacities ( lylihalas, 1978) for various elemental 
mixtures using the Opacity Project cross sections. These 
are LS coupling calculations including channel couplings, 
carried out for ground states and excited states for all 
ions of elements with Z u p to 30. A widely used com - 
pilation was presented bv IVerner and YakovlevI l)l995(l . 
who calculated central field cross sections including in- 
ner shells based on the Dirac-Slater potential, for most 
astrophysically important ions. These were scaled to fit 
the behavior of the i?-matrix calculations near thresh- 
old, while also having the correct behavior in the high 
energy asymptotic limit. These are convenient to use, 
and provide an approximation to the resonance structure 
near threshold, and therefore can be used to calculate 
integral quantities such as photoionization rates. How- 
ever they do not contain the detailed resonance structure 
near threshold, and so should not be applied to calcu- 
lations of monochromatic opacity or synthetic spectra. 
Opacity Project opacities have been recently revised to 
include inner-shell contributions, using the AUTOSTRUC- 
TURE package, and these are available thr ough the Opac- 
ity Project on line database TIPTOPbase fBadnel l et all 
12005) . Figure Ha shows a com parison of BPRM calcula- 
tions for O''"'" with experiment jChampeaux et all . l2003|) 
in the vicinity of the threshold for ground state ioniza- 
tion, showing the extensive resonance structure and gen- 
eral consistency between the measured and calculated 
cross section. 

Even when the total ionization rate is so low that the 
gas is neutral, or if photoionization is not the domi- 
nant ionization mechanism, the spectrum of the radia- 
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FIG. 13 Comparison of BPRM calculation of the photoi on- 
ization of O*'*' with experiment ijChampeaux et al.L 120031) in 
the energy range 90-130 eV. Below 103 eV the cross section 
is due to ionization from the 2s2p excited configuration, and 
the 2s^ above, from lNahaj (|2004f) . 



tion field transmitted in the X-ray band is affected by 
the photoionization cross section. Thus a knowledge of 
the cross section is needed to interpret observations of 
X-ray absorption. In this case, spectroscopic accuracy 
can be important, since features in the cross section can 
be used to diagnose the conditions in the absorbing gas 
ijPaerels et al.U200Hl . Examples of such features include 
the line features due to K shell photoexcitation in oxy- 
gen and its ions, which is abundant in the interstellar 
medium. T he importance o f this process has been em - 
phasized bv iPradhanI (|200^ and 'Pradhan et al.' ('2003). 
An ear ly detection of interstellar oxygen K absorption 
was bv lSchattenburg and Canizared l)l986|) . 

The calculation of photoionization from inner shells is 
affected by the resonance structure associated with ex- 
citation of states with one rvp excited electron and a K 
shell vacancy. They decay predominantly by spectator 
Auger transitions, in which the electron does not par- 
ticipate, and have widths nearly independent of principle 
quantum number n. This leads to a series of resonances 
with constant width close to the threshold for photoion- 
ization, thereby smearing and lowering the location of the 
edge. This effect, called A uger damping, has a lso been 
observed in the laboratory ijParhat et al.L[l99^ . It was 
emphasized and pointed out earlier in calculations of in- 
ner shell ionizati on by Gorczvca and McLaug hlin ( 200(|) 
for oxygen (c. f.-jCorczvca and Robicheaua <fl999)l. and 
Ne llGorczvcal l2000() . Computation of these effects us- 
ing an i?-matrix technique relies on the use of an opti- 
cal potential to mimic the decays to spectator channels 
not explicitly included in the close coupling expansion. 
BPRM photoionization cross sections including Auger 
dampin g have been carried out for inner shells of all ions 
of iron feautista et al 2003. 2004: Kallma n et al. 200i 
iMendoza et alll2004lPalmeri et allTi5Ml2003allJ) . and 
for oxygen llGarcia e t al.'.'2005'). An example of these ef- 
fects is shown in figure ^1 which shows the photoioniza- 
tion cross sections of Fe^^+ and Fe^^+ near t he thr eshold 
for the K shell, taken from iBautista et all l)2004|) . The 
left panels show the resonance structure with the effects 
of damping by resonances with n > 2 included, while the 
right panels show the resonance structure with damping 
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FIG. 14 High-energy total photoabsorption cross sections of 
the ground level of: (a) Fe^^'"'" including radiative and Auger 
damping effects; (b) Fe^®+ excluding damping effects for res- 
onances with n > 2; (c) Fe^'^"'" including damping; (d) Fe^^""" 
excluding damping for resonances with n > 2. This demon- 
strates that when damping is included the resonance widths 
are c onstant for high n lea ding to a smearing of the K edge. 
From lBautista et alJ J2004|) . 



by resonances with n > 2 excluded. 

The vacancies created by inner shell photoioniza- 
tion are filled by fluorescence or Auger decays. The 
rates for these processes have been calculated by 
iJacobs and Rozsnv af flQSG") using configuration-average 
energy levels and by Kaastra and Mcwc ( 199^ based on 
isoelectroni c scaling from the calculations of iMcGuirj 
Ill969tll97(iri971alhl.ll972<) . The validity of these approx- 
ima tions has been discuss ed by iGorczvca et al.l l|2fl0,'^ 
and iGorczvca et "aP l)2006f) . Level resolved intermediate 
coupling calculations of K Auger and fluorescenc e rates 
have been calculated f or iro n bv iBautista et al 1 2003 , 
20041) ■ iMendoza et all l)2004|) . and iPalmeri et alJ l)2002 , 



2008, 



2. Radiative Recombination (RR) 

Radiative recombination (RR) is more important as a 
line emission mechanism in photoionized plasmas than 
in coronal plasmas and therefore receives more discus- 
sion in this context. Rate coefficients can be calculated 



using detailed balance argumen ts via the Milne relation 
ijOsterbrock and Ferlandll2006() . The rate coe f ficient can 
be expressed in the form defined by ISeatoiJ l)l959D for 

hydrogenic ions: 



S^{hvfe-^''/''^ar^{Z, hy)d(hy) (34) 

where /„ is the ionization potential and the photoioniza- 
tion cross section can be written 



cr„(Z, hv) 



2^a7rao n 



{l + n'e)-'gn{n,t) (35) 



where a is the fine structure constant, oq is the Bohr 
radius, and energy conservation requires that hu = 
hRcZ^(l/n^ + e). The Kramers-Gaunt factor can be ex- 
panded as a polynomial in u = n?e, and the rate coeffi- 
cient can be expressed as: 



a„{Z,T) 



Ai/2 

DZ x„5„(A) 



(36) 



where D = ^ (|)^ ^a^'cag 
A = hRcZ^/{kT) = 157890ZVT, x„ = \/r 



5.197 X lO^'^cm'^ s"^ and 



SniX) = 



gii{n,e)e 



du 



(37) 



and u 



n € 



ISeatonI lll959D uses the asymptotic expan- 
sion for gii{n,e) to evaluate S'„(A). Numerical tables 
for S'„(A) and for the total recombination rate coeffi- 
cients for as {Z, T) = E5^ia„(Z, T) and for the sp ectrum 
of reco mbining electrons have been published bv ISeatonl 
l|l959l) . When RZ'^jr? > kT then a„(l,T) oc n-^T-^l"^ 
and when RZ'^/n^ < kT then, for Z = 1, 



a„(l,T)cx n-^T-^' 



57890 



0.5772 



-H8.56 X IQ'^T^/^ - 2.3 x lO^^T^/^) (38) 

The first term in the brackets becomes less important at 
high temperatures, corresponding to the fact that excited 
levels contribute more at low temperatures than at high 
temperatures. 

Such calculations are typically tabulated as total re- 
combination rates, in which all the possible radiative 
transitions from the lowest continuum state (i.e., the 
ground state of the next highest ion stage) to the bound 
states of a given ion are summed. This requires a sum 
over all the rate coefficients into excited levels of an ion, 
and implicitly assumes that all recombinations to excited 
levels decay to the ground level. Such sums traditionally 
use ground state photoionization cross sections appropri- 
ate to the ion, together with hydrogenic photoionization 
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cross sections for the excited levels. iTarted l)l97lf) calcu- 
lated RR rate c oefficients for many ions of astrophysical 
interest, a s did ICouldl lll97(^. A widely used com pila- 
tion is by 'Aldrovandi and Peauignol] l)l973l Il976|) who 
calculated RR and DR rate coefficients for m any ions 
of ast rophysical interest, calculated using the iBurgessI 
11965") general formula for DR and the Milne relation. 
[Peauignot ct al. (1991J calculated total and effective RR 
coefficients for all ions with Z<10 f or important optica l 
and UV transitions of these ions. IWoods et alJ lll98lD 
calculated RR rate coefficients for iron , using the cross 
sections of lReilman and MansoiJ l)l979l) . iPradhaiJ l)l983|) 
calculated the RR for the ground states of Li-like ions. 
IVerner and Ferlandl l)l99(j) calculated rate coefficients for 
RR based on Op acity Project photoionization cross sec- 
tions. Qui ll2003cr) has provided total RR rate coefficients 
for the H-like through Ne-like isoelectronic sequences for 
the 7 e lements Mg, Si, S, Ar, Ca, Fe and Ni, and 'GJ 
has provided X-ray line emission rate coefficients 
for Fei^+ - Fe^^+ due to both RR cascades and DR. The 
need for tabulations of total rate coefficients or of effec- 
tive rate coefficients is reduced by advances in computer 
speed, which allow plasma modeling codes to calculate 
or make use of state-specific recombination rates. 



3. Dielectronic Recombination (DR) 

Dielectronic recombination (DR) in photoionized plas- 
mas is generally less important than in coronal plasmas, 
but it is not negligible. As shown in figure IHl the contri- 
butions of radiative and DR are comparable for Fe^^+ DR 
at temperatures characteristic of photoionized plasmas. 
For F e^^"*" - Fe^^+ DR dominates RR by a factor of 2 or 
more llSavin et al.1 . '2002a"b. •2003'). Theoretical rate coef- 
ficients are less certain at low temperatures, owing to the 
fact that the electron velocity distribution samples low 
energy states. Theory cannot calculate the resonance en- 
ergies for the relevan t DR re s onanc es at l ow energy with 
suffici ent precision. iNetzed ()2004j) and iKraemer et al.l 
l)2004j) provide examples of the sensitivity of astrophys- 
ical results to DR rates in photoionized plasmas, and 
also point out the consequences of errors in the rate co- 
efficients which are in widespread use for astrophysical 
modeling. 

The importance of low temperature DR was pointed 
out bv lStorevll)l98l|) . who showed that the 0^+ 2p'^ 
2s2p ^PoA2296 line observed from planetary nebulae 
is likely excited by DR. This was extended t o ion s 
of C, N, and O by iNussbaumer and StorevI l)l983(l . 
who calculated rates for these ions using LS cou- 
pling. The associated satellite spectra were calculated by 
[Nussbaumer and Storev ( 19 84). These were extended to 
M g, Al, Si by i Nussbaumer and StorevI l)l986|) and to Ne 
bv INussbaumer and StorevI l)l987|) . Calculations of low 
temperature DR are most affected by the uncertainties 
in the low energy resonance structure, and this in turn 
is affected by the assumption of LS coupling in many 



calculations. LS coupling can be applied when the fine 
structure splitting is small compared to the term separa- 
tion which is caused by the Coulomb interaction. 

Rece nt ex p erimen t al measurements u sing storage 
rings l|SavinL l2000l: ISavin et all l2002bl) and calcu- 
lations using autostructure Altun et all 2004 , 



lations usm g auto struc ture HAitun et ai.1. zvm , 
2005', '200^; Badnell", '2006': 'Colgan et al.', '200?, '2004 , 



2005; Mitnik and BadncU, 2004; Zat sarinnv et al., 2003, 
2004albl . l2005albl l2006|) of DR include intermediate cou 



pling structure calculations, and so are likely to be more 
accurate for low temperature DR than previous calcu- 
lations. None theless, as p ointed out bv lSchippers et afi 
l)2004j) and bylc] l|2003bfl . DR at low temperatures con- 
tinues to be a major source of uncertainty in calculations 
of ionization balance under photoionization conditions. 



4. Charge Transfer 

Background 

Charge transfer is of potential importance in photoion- 
ized plasmas owing to the fact that the charge transfer 
rate coefficient for an ion with neutral H or He can ex- 
ceed that for electro n recombination by a factor > lO'' 
l)Kingdon and Ferland . 1996). Photoionized plasmas can 
have a hydrogen neutral fraction as great as ~ 0.1 co- 
existing with several times ionized metals, depending on 
the shape of the ionizing spectrum. Another application 
for charge transfer rate coefficients was pointed out by 
ISavin et all J2004a . who calculated the rate coefficients 
for H+-I-H2 H-|-H^ using recently published theoreti- 
cal cross sections, and showed that uncertainties in these 
rate coefficients can affect the cooling and formation of 
primordial structures such as stars and galaxies. 

The importance of charge trans fer in the interstel- 
lar m ediu m was pointed o ut by iField and SteigmanI 
l)l97lj) and ISteigman et al.1 l)l97l(l who also calculated 
rate coefficients for C^+ -I- He ^ C"*" -I- He"*" l|SteigmanL 
Il975^ and discussed the influence of charge transfer 
on the O"*" and N+ ionization balance. The po- 
tential importance of charge transfer and, for the 
first time, charge transfer involving ions more than 
singly ionized on the structure of planetary nebula 
was pointed o ut bvlPeoui gnot. S tasinska and Aldrovandil 
(1978) and iPeouignotI (1980a'). This led to ef- 
forts to calculate cross sections and rates, includ- 
ing calculations of the charge transfer of N"*" with 
H, and of C+ and S+ with H and He usine 



calculation by 


'Butler and Dalgarnc| 


iButler and Dalgarno 




.1980b) calcu- 



lated the charge transfer of a variety of multiply 
charged ion s with H and He using a La ndau-Zener 
calculation. IWatson and ChristcnsciJ l)l979(l calculated 
charge transfer of C'^+ -l-H and N^+-|-H using close- 
coupling appro ximation. A review and compilation by 
[Kingdon and F erlandl l)l996|) provides a collection of rate 
coefficients in a form useful for modeling. 
Recent Developments 
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H llStancil et all 12001 . 
1)2003(1 . These show large (> 



The most accurate calculations at the low energies of 
interest to astrophysics (< 1 eV/amu) are techniques sim- 
ilar to molecular calculations, usin g the molecular-o rbital 
close-coupling (MOCC) technique l|Heil e t al.','l98 d). Re- 
cent MOCC calculations include 0^+ -|- H (,Wang et all 
l200.'^. Si4+ He IStancil et all n997^. + H 

llZvgelman ei~all 11 997). 8^+ 
and + He lWan g et al^ "~ 

10^) differences with Landau-Zener calculations in some 
cases. Recent measuremen ts include Si'^+_Si^J_j__Si^ 
+ He l|Tawara et all I2OOID. Ne^+ + R llReioub et all 
I2004D. Ne^+ + H llMroczkowski et all 120031) C'^+ + H 
("Bhe k et allll99'i) N^+ + H llPieksma et all 1199/1) , C+ 
+ H (ISta,Ticil et~ll99SalH^. 

A review of the charg e transfer process is given by 
iJanev and Winteil l|l985tl . Recent reviews of charge 
transfer data for photoionized plasmas are g iven in 
IStanciJ l(200ll) for theory and in iHaveneJ l)200l() for ex- 
periment. Although quantum mechanical MOCC calcu- 
lations provide the most r eliable results when no experi- 
mental data are available l|Rta,nci1I . Eo?)Tl) . still large dis- 
crepancies between theory and experiment at collision 
energies of less than a few hundreds eV/amu point to 
the difficulty of theory at this energy range an d the need 
for benchmark measurements l|Haveneill200l|) . 



C. Charge Transfer in Cometary and Planetary 
Atmospheres 

The discovery of X-ray emissio n from comets by the 
ROSAT X-ray astronomy satellite l|Lisse et allll996j) led 
to the appreciation of the importance of charge transfer 
of solar wind nuclei with neutrals in comets and gaseous 
planets. This process has a distinct spectral signature as- 
sociated with the cascade of the captured electron from 
the excited state where capture initially occurs. The 
most probable states are determined by the energetics 
o f the collision and the level structure of the reactants. 

iHaberli et al. extended this idea to look at in- 

dividual spectral lines produced from C, O, and Ne ions. 
The excited state of the ion following capture was as- 
sumed to have a principal quantum number n which 
is the nearest integer to the quantity g*^'^^, where q is 
the charge on the ion, followed by a cascading decay 
of Art — 1. This assumption leads to X-ray spectra 
which fits to the observed X-ray spectra in a subset of 
the cases. A recent review of t he physics of co metary 
X-ray emission was presented by ICravensI (|2002l) . EBIT 
si mulations of cometa r y cha rge exchange were performed 
bv iBeiersdorfer et all l)2003j) . Additional models for the 
cascade, level populations and X-ray line emission associ- 
at ed with cometary cha rge tra nsfer h ave been calculated 
bv iKharchenko and Dalgarnol l)2000l) . 

Measurements for various species have been carried 
out using beam techniques. T otal cross sect i ons fo r 
C, N and O are ta bulate d bv iPh aneuf et all l|l987(l . 
iHavener et all l|l989|) . and iHuo et al.i ((198^ For in- 



terpreting X-ray spectra, and for testing theoretical 
calculations, state-selective measurements are needed. 
This need, along with the relatively low collision en- 
ergy, make the X-ray astronomy data needs distinct 
from those for fusion plasmas. Techniques for mea- 
suring state-selected charge transfer include energy loss 
spectroscopy and photon emission spectroscopy. En- 
ergy loss spectropscopy provides greater counting rates, 
while photon emission spectroscopy which allows higher 
energy resolution. Measurements have been made us- 
ing energy loss spectrosc opy of 0^+ in H, H2 and He 
l)McLaughlin et a ll |1990 |), C^"*" ions in collisions with 
H-, and O^ llMcLaughlin et all 119921), Fe'^+ an d Fe^+ 
ions with H and He atom s (lMcLMieJilili-e t all ^993), 

with H and H2 fWil son et all Il990bl) . and 
with H, H2 an d He (.Wilson et al.Ul990a^ . Using sim- 
ilar techniques, iKimura et al.l l|l987(l have studied highly 
stripped Nc, O, N and C in collisions with H and H2, and 
Kambcr et al. ( 1996) have studied colliding with H2, 
He, Ne, and Ar. Beam measurements have been made 
using photon emission spectroscopy to obtain state se- 
lective cross sections of single electron cha rge transfer 
of C''+ on H and H^ jHoekstra et allll990^ . 0^+ on H 
and H2 llBeiiers et allll99fil). C/+ . . 0^+ with H, H2 
and He (iDiikkamp et all Il985al). C^+. ^+, and 

Ne^+ onto H^. He. Ar ijPiikkamp et all Il985b|) . Sim- 
ilar techniques ha ve been used to stud y C^"*" and N®+ 
with He and H2 llSuraud et all 
N'^+ with He and H9 l|Ciric et all 



1991), and and 



1985(1 . Double charge 



exchange can be an important loss process for highly 
charged ions incident on neutrals, but in many cases it 
leads to the production of autoionizing states, and there- 
fore does not contribute to X-ray emission. Cross sec- 
tions with molecular targets other than H2 have been 
measured for ions of C, N, O, and Ne collidi n g wit h 
He, H2, CO2, and H2O by ICreenwood et alJ 1I2OOII) . 
Some of these experimental data have been collected 
at the ORNL/UGA charge transfer database webs ite. 



http: / / cfadc.phy.ornl.gov/astro/ps/data/home.htmH in- 
cluding both total and state-specific cross sections and 
rates, along with fits to many of these. 

Computational techniques for charge exchange colli- 
sions i nclude the classical tr ajectory Monte-Carlo tech- 
nique jCornefius et alll2000j) . in which the motion of the 
projectile, target and electron are calculated by integrat- 
ing Hamilton's equations including only the Coulomb in- 
teractions. This technique is of limited use at energies 
below 1 keV/amu, where the structure of the projectile 
ion is important. Another method is the semi-classical 
impact parameter close coupling method, in which the 
wavefunction of the electrons is calculated from a set of 
orbitals appropriate to the projectile ion. This technique 
is mos t suitable to systems with a sma ll number of elec- 
trons l|Kumar and SahaL Il998t ISahaL Il995() . and so has 
not been applied to the study of many ions. The Landau 
-Zener model has also been widely applied, and involves 
the use of avoided level crossings in order to determine 
the cross sections for various processes. It turns out that 
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especially for energies below 200 eV/amu there are con- 
siderable differences between experiment and theory both 
in relative and state-selective cross sections. 



VIII. DISCUSSION 

Basic needs and physical processes important to 
atomic data for X-ray astronomy have been known since 
the advent of solar X-ray astronomy. This has led to 
steady production of rate coefficients and cross sections 
for use in modeling and interpreting observed astrophysi- 
cal X-ray spectra, and with considerable overlap with the 
needs for fusion plasmas. Prior to the launch of Chandra 
and XMM — Newton a great deal of work had been 
done toward the goal of developing an atomic database 
for coronal plasmas. Much of this was devoted to the 
study of discrete diagnostics, since the spatial and tem- 
poral structure of the X-ray emission limits the accuracy 
of global modeling for the Sun. 

The launch of instruments capable of observing spec- 
tra from extrasolar objects with comparable resolution 
and good counting statistics has added motivation for 
data production and it has changed the emphasis some- 
what from that of studies of solar X-rays. Extrasolar 
objects often cannot be spatially resolved and the sen- 
sitivity to temporal variability is limited. Arguments 
based on other knowledge of their properties motivates 
attempts to construct global models based on a single set 
of physical conditions, such as the assumption of a single 
temperature or a cooling flow for the coronal emission 
from a cluster of galaxies. For some objects a simple dis- 
tribution of conditions is assumed, such as a differential 
emission measure analysis performed for a cool star. In 
addition, the spectra of distant objects often do not have 
sufficient statistical accuracy to allow the application of 
detailed discrete diagnostics. This creates a greater need 
for global modeling, in which the ionization balance and 
spectrum can be calculated based on simple assumptions 
about the conditions, such as temperature and density. 
In addition extrasolar data have broadened the range of 
physical processes of interest to include charge exchange, 
photoionized plasmas, inner shell processes, and opacities 
associated with interstellar gas. 

Significant advances in the calculation and measure- 
ment of atomic cross sections and rate coefficients needed 
for X-ray astronomy have occurred in parallel with the 
launch of the new observatories. Notable among these are 
the energetic application of experimental apparatus such 
as the EBIT and storage rings, and the improvements 
in computer technology and campaigns to calculate large 
quantities of data such as the Iron Project and Opacity 
Project. At the same time, long term efforts to compile 
reliable data from more traditional laboratory sources 
have continued to yield results, and the advent of free 
on line databases, such as the NIST database, has aided 
in their use. Other important databases are those devel- 
oped primarily for the fusion energy program at ORNL 



('http: / /www-cfadc. phy.ornl.gov/) and ALADDIN and 
AMBDIS at the IAEA (http://www-amdis.iaea.org/). 
As a result, there now exist accurate experimental 
datasets for line wavelengths and cross sections or rate 
coefficients for some key processes, along with computa- 
tions using the most accurate known algorithms for many 
quantities. Laboratory measurements have been made 
for DR of many ions, and comparison with calculations 
allows benchmarking of computational techniques. Com- 
prehensive calculations of state-specific DR, which have 
been benchmarked against the measurements, are avail- 
able for many ions of interest. Measurements of EII exist 
for most ions of astrophysical interest, although verifica- 
tion is needed for many of these. Calculations, including 
close-coupling calculations with adequate treatment of 
intermediate coupling, relativistic effects, CI and radia- 
tion damping, have been made for radiative transition 
probabilities and electron impact collision strengths of 
many ions of interest. Measurements of absolute elec- 
tron impact excitation cross sections have proven to be 
crucial in benchmarking these calculations. Beam mea- 
surements have yielded state-selective cross sections for 
charge exchange for many ions of interest to solar system 
X-ray studies. 

In addition to the databases discussed above, a great 
deal of useful data has been collected and made pub- 
licly available as part of the databa ses a ssociated with 
analys is packag es such as CHIANT I l|Land i and Bhati^ 
l2005a^. APEC jSmith et al. . 20 0j|), and mekal/spex 
ijKaastra. Mewe and^ieuwenhuiizenl Il996j) . It is impor- 
tant to point out that the ultimate source of all databases 
and compilations is extensive computational and exper- 
imental work, and that the credit for this work is of- 
ten neglected when the compilation or database is used. 
Astrophysicists and modelers should, whenever possible, 
attempt to cite original sources even when using the com- 
pilation or database as a guide or repository for atomic 
data. 

As a result of accuracy demanded by the new instru- 
ments for X-ray astronomy some of the atomic calcu- 
lations which were in widespread use for modeling as- 
trophysical plasmas are no longer adequate for applica- 
tion to many observations. The approximations necessi- 
tated by early computers or analytic work provide valu- 
able insight and constraints on the more recent work, and 
therefore have been crucial in the development of modern 
tools. However, there are few remaining processes or ions 
for which rate coefficients calculated using simple approx- 
imations are all that is available. These include the use 
of the Born approximation, or CBO, in coUisional ion- 
ization cross sections, the Bethe approximation for coUi- 
sional excitation, ab initio wavelengths for strong lines in 
the X-ray band, the Burgess general formula for DR, hy- 
drogenic or central field cross sections for photoionization 
or radiative recombination, and Landau-Zener rate coeffi- 
cients for charge exchange. An added consequence of the 
capabilities of new computers is that it is no longer nec- 
essary to use total recombination rates when calculating 
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ioniza tion balance. Rather, coUisional-radiative models 
(eg. (jSummers et all l2006j) ') can be quickly calculated 
which take into account state-specific recombination and 
ionization rates, and so are applicable to a wide range of 
gas densities and radiation environments. 

Areas where there are still critical needs include the 
accumulation of energy level structures and transition 
wavelengths which are of sufficient accuracy. These 
are needed for applications including line identification, 
where the observations have an accuracy of 10~^ in many 
cases, and for calculations of DR at low temperature. Ex- 
periments have the most promise for useful work in this 
area, but new theoretical techniques are needed for this 
challenge as well. There are few experimental measure- 
ments of inner shell photoabsorption for ions of interest 
to astronomy. Similarly, many lines in observed spectra 
such as that of NGC 3783 remain unidentified, and these 
may be associated with inner shell transitions not typi- 
cally observed from coronal plasmas. X-ray processes in- 
volving molecules, including detailed spectra associated 
with charge transfer, and inner shell opacities, have not 
been studied for many species. Experimental work is 
needed to continue the campaign to measure all the DR 
and CI rate coefficients needed to benchmark the cal- 
culation of coronal ionization balance in both electron 
ionization and photoionization driven plasmas, and in 
cometary and planetary atmospheres. 

The most probable choices for future X-ray instrumen- 
tation will likely be only a partial continuation of the 
trends of the recent past. Both technological challenges 
and astronomical priorities suggest that X-ray astronomy 
satellites following Chandra and XMM — Newton will 
have significantly greater sensitivity, so that the number 
of objects which can be observed spectroscopically, the 
statistical quality of the spectra and the ability to study 
time variability will all be greatly enhanced. However, 
it is less likely that either the spectral resolution or the 
spectral bandpass will be improved in such instruments. 
Thus, spectral resolution comparable to that available 
from the best optical or UV instruments will not be at- 
tained. In view of this it seems likely that the interpre- 
tation of data from such future missions will continue to 
rely on tools similar to those currently in use, tools which 
attempt to calculate ionization, excitation and synthe- 
size the spectrum over a wide range of wavelengths. New 
targets for these observations will likely include galax- 
ies, clusters of galaxies, and intergalactic gas in order 
to study the formation of structure, nucleosynthesis, and 
cosmological parameters. More sensitive instruments will 
also study spectra of fainter nearby objects such as stars 
and supernova remnants, and allow improved statistics 
and more detailed study of time variability of the bright- 
est objects. This is in contrast to what might be expected 
if the trend toward both greater spectral resolution and 
greater sensitivity were to occur, which might lead to 
more extensive application of discrete diagnostics and a 
reduced reliance on global modeling. If so, the atomic 
data needed for the foreseeable future will not be greatly 



altered, in the sense that the precision required will be 
comparable to the best currently available, and physical 
processes of interest will span those covered in this re- 
view. Spectra with improved statistical accuracy, or time 
resolution, are likely to reveal physical effects which have 
not been incorporated into the available atomic database. 
These may include departures from ionization equilib- 
rium, non-stationary processes, optical depth effects, the 
effects of magnetic fields, and interactions with energetic 
particles and non-Maxwellian electron energy distribu- 
tions. This implies the need for cross sections at ener- 
gies both significantly above and below the traditional 
range of energies prescribed by coronal equilibrium. In 
addition, it is likely that X-ray astronomers will further 
explore molecular, or solid matter in astrophysics using 
charge transfer and also absorption features near edges 
or inner shell lines. 
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TABLE II Glossary 



Acronym 



Meaning 



AMA Angular Momentum Average 

APEC Atomic Physics Emission Code 

BP Breit-Pauli 

BPRM Breit-Pauli R-Matrix Code Package 

CB Coulomb Born Approximation 

CBE Coulomb Born with Exchange 

CBO Coulomb Born Oppenheimer 

CCD Charge Coupled Device 

CI Configuration Interaction 

DI Direct Ionization 

DR Diclectronic Recombination 

DW Distorted Wave Aprroximation 

EA Excitation-Autoionization 

EBIS Electron Beam Ion Source 

EBIT Electron Beam Ion Trap 

ECIP Exchange Classical Impact Parameter 

EII Electron Impact Ionization 

FAC Flexible Atomic Code 
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